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ABSTRACT 

The  emission  spectra  of  40  amp-2.5  atm  xenon  and  krypton  arc 
discharges  are  calculated.  The  spectral  overlap  of  these  spectra  and 
the  measured  emission  spectra  of  tungsten  and  potassium-mercury  arc 
discharges  with  the  measured  absorption  spectra  of  Nd:YAG  are  used  to 
determine  relative  system  efficiencies  which  agree  well  with  measured 
values.  It  is  found  that  the  K-Hg  lamp  is  four  times  more  efficient 
in  pumping  Nd:YAG  than  the  studied  W,  Kr,  and  Xe  lamps.  The  spectral 
overlap  calculation  for  a  typical  ruby  oscillator-amplifier  at  the  peak 
of  the  flash  lamp  pulse  gives  a  pump  efficiency  which  is  much  higher  than 
measured  system  efficiencies.  This  discrepancy  may  be  caused  by  the  time 
dependent  behavior  of  the  flash  lamp  pulse  and  the  oscillator  pulse. 

A  new  technique  is  developed  to  obtain  electronic  radial  wave  functions  from 
known  radial  Integrals.  In  a  test  case,  the  Cu  4s  wave  function  is  determined 
within  +  8%.  Group  theory  is  applied  to  the  calculation  of  radial  integrals 
involving  hydrogen  wave  functions  and  some  selection  rules  are  explained. 


i 


TABLE  OF  CONTENTS 


PAGE 

1.  INTRODUCTION .  1 

Nd.'YAG  LASERS  5 

2.  LAMP  EMISSION  SPECTRA .  5 

2.1  Theoretical  Results  for  Xenon  and  Krypton .  7 

2.1.1  Arc  Composition . . .  8 

2.1.2  Transport  Properties.... .  9 

2.1.3  Line  Radiation .  11 

2.1.4  Continuum  Radiation .  13 

2.1.5  Temperature  Profile.. .  15 

2.1.6  Arc  Emission .  16 

2.2  Comparison  of  Theoretical  and  Experimental  Results .  1? 

2.2.1  Xenon- .  18 

2.2.2  Krypton .  20 

2.3  Measured  Spectra  for  W  and  K-Hg  Lamps .  21 

2.4  Discussion. .  22 

3.  LASER  ABSORPTION  SPECTRA .  26 

4.  SYSTEM  EVALUATION. .  27 

4.1  Lamp  Efficiency- . 27 

4.2  Slab  Efficiency*. .  28 

4.3  Pump  Efficiency. . 29 

4.4  Activator  Efficiency . . .  52 

4.5  Discussion . 33 

5.  SUMMARY .  34 

iii 


TABLE  OF  CONTENTS  (Continued) 


* 


PAGE 

RUBY  LASERS  37 

6.  INTRODUCTION. .  38 

7.  LAMP  EMISSION  SPECTRUM. .  39 

8.  RUBY  ABSORPTION  SPECTRA .  40 

9 .  SYSTEM  EVALUATION. .  41 

9.1  Lamp  Efficiency. .  41 

9.2  Slab  Efficiency .  42 

9.3  Pump  Efficiency .  43 

9.4  Activator  Efficiency .  43 

10.  SUMMARY .  45 

RADIAL  WAVE  FUNCTIONS  47 

11.  PARAMETER  FITTING  FROM  SLATER  INTEGRALS .  48 

11.1  PF  Constraints .  49 

11.2  PF  Basis  Sets .  51 

11.3  Numerical  Procedures .  53 

11.4  Screening  Constants .  57 

11. 5  Application  to  the  Cu  4s  Electron. .  59 

12.  GROUP  PROPERTIES .  63 

13 .  SUMMARY . . . . .  .rrr. . . . . .  1 1 . .  .T  65 

APPENDIX  1  -  A  Comparison  of  Lamps  for  Use  in  High  Continuous 

Power  Nd:YAG  Lasers .  68 

APPENDIX  2  -  Rod  to  Slab  Transformation .  72 

APPENDIX  3  -  Newton' 8  Method  for  Parametric  Fitting  of  Radial 

Wave  Functions  in  Integrals .  76 

iv 


TABLE  OF  CONTENTS  (Continued) 


PAGE 

APPENDIX  4  -  Comparison  of  Direct  Search  and  Newton’s  Method .  81 

APPENDIX  5  -  Group  Theoretical  Properties  of  Radial  Wave 

Functions .  84 

APPENDIX  6  -  Group  Properties  of  Hydrogenic  Wave  Functions .  92 

REFERENCES .  101 

TABLES .  105 


FIGURES, 


128 


May  31,  1970 


OPTICAL  PUMPING  OF  LASERS 

L.  H.  Taylor,  I.  Liberman, 

R.  Liebermann,  A.  Federowicz,  L.  Armstrong* 

1.  INTRODUCTION 

The  optically-pumped  solid-state  laser  provides  the  highest 
peak  power  from  coherent  sources.  However,  one  of  the  main  energy  losses 
in  optically  pumped  lasers  is  in  the  energy  transfer  between  the  optical 
pump  and  the  laser  medium.  This  report  s'lmmarizes  the  progress  toward 
measuring  and  optimizing  this  energy  transfer. 

The  energy  transfer  can  best  be  optimized  by  maximizing  the 
spectral  overlap  between  the  emission  spectrum  of  the  optical  pump  and 
the  absorption  spectrum  of  the  lasing  medium.  Since  the  lasing  medium 
is  optimized  for  its  lasing  properties,  it  is  unlikely  that  its  absorp¬ 
tion  spectra  can  be  significantly  altered  without  reducing  its  desirable 
lasing  properties.  On  the  other  hand,  the  single  job  of  the  optical 
pump  is  to  convert  electrical  energy  to  optical  energy  in  such  a  manner 
that  the  lasing  medium  absorption  of  the  optical  energy  is  maximized. 

It  then  follows  that  the  best  approach  is  to  change  the  optical  pump 
until  the  spectral  overlap  is  maximized  for  a  particular  losing  medium. 
The  most  common  optical  pump  is  a  wall-stabilized  arc  discharge,  and 
therefore  this  type  of  discharge  is  the  subject  of  this  report. 

The  work  in  this  contract  is  divided  into  two  Phases.  In  the 
first  Phase  the  emission  spectra  of  xenon  and  krypton  arc  discharges 
are  calculated  and  measured.  This  permits  an  evaluation  of  the  theo¬ 
retical  results.  The  measured  absorption  spectra  of  Nd:YAG  and  ruby 
laser  slabs  are  then  converted  to  digital  form  and  the  pump  efficiencies 
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of  different  lamps  for  these  laser  materials  are  calculated.  In  the 
second  Phase  an  atomic  radial  wave  function  is  expanded  in  a  power 
series  with  fitting  parameters,  and  the  fitting  parameters  are  then 
adjusted  to  reproduce  the  Slater  integrals  determined  from  spectroscopic 
measurements.  Thus  the  radial  wave  function  is  "experimentally"  deter¬ 
mined  and  is  independent  of  the  atomic  shell  structure  complexities. 
Moreover,  the  group  theoretical  properties  of  hydrogenic  radial  wave 
functions  are  investigated  to  establish  group  theoretical  methods  for  the 
evaluation  of  radial  integrals. 

These  two  Phases  complement  each  other.  The  pump  efficiency 
calculation  will  be  possible  for  both  experimental  and  theoretical 
spectra.  However,  at  present  the  emission  spectra  of  arc  discharges  can 
only  be  calculated  for  the  noble  gases  and  the  alkali  metal  gases.  When 
the  second  Phase  is  successfully  developed,  the  emission  spectra  will  be 
calculable  for  any  monatomic  gas  provided  only  that  the  energy  levels  of 
the  emitting  atoms  are  known  (this  condition  already  exists  for  the 
first  Phase). 

The  Phase  One  work  in  this  report  is  divided  into  two  parts. 

The  first  part  investigates  the  narrow  absorption  line  Nd:YAG  laser 
whereas  the  second  part  investigates  the  broad  absorption  band  ruby  laser. 
In  each  part  the  lamp  emission  spectra,  the  slab  laser  absorption  spectra, 
and  the  useful  overlap  of  these  spectra  are  measured  or  calculated.  In 
the  first  part  the  emission  spectra  of  a  xenon  arc  discharge  and  of  a 
krypton  arc  discharge  are  calculated  and  compared  to  experimental  measure¬ 
ments.  The  spectral  overlaps  of  these  calculated  spectra  and  of  the 
tungsten  and  potassium-mercury  measured  emission  spectra  with  the  measured 
0.7%  Nd:YAG  absorption  spectra  are  then  computed.  The  calculations  thus 
evaluate  the  lamp  emission  efficiencies,  the  Nd:YAG  slab  fluorescence 
efficiencies  when  pumped  by  the  various  lamps,  and  the  total  system 
efficiencies.  The  latter  two  efficiencies  are  evaluated  as  a  function 
of  the  slab  thickness.  Since  most  laser  materials  are  in  the  form  of  a 
cylindrical  rod,  a  relationship  between  slab  thickness  and  rod  diameter  is 
given . 


In  the  second  part  the  spectral  overlap  for  a  typical  ruby 
laser  is  calculated  from  the  measured  spectrum  at  the  peik  radiance  of 
a  xenon  flash  lamp  and  the  measured  absorption  spectra  (there  are  two 
spectra  because  ruby  is  anisotropic).  The  efficiencies  are  calculated 
from  these  spectral  overlaps.  The  difference  between  a  four-level  laser 
such  as  Nd:YAG  and  a  three-level  laser  such  as  ruby  are  explicitly 
considered. 

Phase  Two  of  this  report  is  concerned  with  the  evaluation  of 
radial  wave  functions.  The  energy  levels  of  any  free  atom  or  ion  are 
experimentally  determined  by  its  emission  spectrum  and  theoretically 
determined  by  radial  integrals  such  as  the  spin-orbit  coupling  constant 
and  the  Slater  integrals.  The  normal  procedure  is  to  parameterize  these 
radial  integrals  and  fit  them  to  the  experimentally  determined  energy 
levels.  The  parameter  fitting  (PF)  technique  developed  and  tested  in 
Phase  Two  adjusts  the  parameters  in  a  power  series  expansion  of  the  radial 
wave  function  until  “he  experimentally  determined  radial  integrals  are 
reproduced.  The  resulting  power  series  representation  of  the  radial  wave 
function  can  then  be  used  to  calculate  the  transition  probabilities  which 
are  required  for  the  prediction  of  emission  spectra  in  l&aps. 

A  more  powerful  alternative  to  the  PF  technique  may  be  group 
theoretical  techniques  for  the  evaluation  of  radial  integrals.  The  group 
properties  of  hydrogen  radial  wave  functions  are  examined  since  they  are 
presumably  the  simplest  wave  functions.  Although  this  investigation  has 
been  short  and  is  not  completed,  it  has  already  produced  significant 
results  which  explain  some  apparent  "selection  rules"  of  radial  integrals. 

To  accomodate  the  reader  who  may  not  be  interested  in  all  of 
this  report,  the  material  has  been  divided  into  three  large  parts: 

Nd :YAG  LASERS,  RUBY  LASERS,  and  RADIAL  WAVE  FUNCTIONS.  Although  the 
three  parts  are  interrelated,  they  are  written  in  a  self-contained  manner 
so  that  the  reader  may  read  them  individually  or  in  any  order. 
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2.  LAMP  EMISSION  SPECTRA 

The  necessary  input  information  to  obtain  the  laser  pumping 
efficiency  of  a  lamp  is  its  spectral  emlttance  and  the  excitation 
spectrum  of  the  lasing  medium.  The  lasing  media  of  NdjY^Al^O^CNdtYAG) 
and  CrrAl^O^Cruby)  will  be  considered  in  later  chapters,  whereas  the 
spectral  emlttance  of  several  lamps  is  considered  in  this  chapter. 

The  lamp  spectra  are  difficult  to  obtain  experimentally. 

Most  continuous  gas  discharges  and  low  to  medium  energy  pulsed  discharges 
have  considerable  line  spectra  which  are  difficult  and  tedious  to  measure 
experimentally  and  then  convert  to  computer  data  cards.  It  would 
simplify  the  work  considerably  if  the  spectra  could  be  calculated  rather 
than  measured  experimentally.  To  do  this  for  an  arbitrary  gas  discharge 
under  arbitrary  operating  conditions  Is  not  possible  at  the  present 
state  of  the  art.  In  fact,  the  major  objective  of  this  contract  is  to 

make  some  headway  in  calculating  the  properties  of  one  class  of  atoms. 

12  3 

However,  because  of  many  years  effort  *  *  we  believe  we  can  calculate 
to  a  reasonable  degree  of  precision  the  atomic  and  ionic  spectra  of  the 
noble  gases  and  the  alkali  metals.  Fortunately,  these  atoms  create 
discharges  which  have  proved  to  be  of  interest  for  punping  lasers. 
Therefore  we  believe  we  can  measure  and  optimize  spectral  overlap 
between  pump  and  rod  quicker  and  better  by  using  calculated  rather  than 
experimental  spectra. 

The  work  done  in  1969  was  to  verify  that  the  calculated 
spectra  compare  favorably  with  the  experimental.  Therefore  careful 
controlled  measurements  of  the  spectral  radiance  of  a  xenon  and  a 
krypton  continuous  arc  were  made.  Computer  calculations  of  these  arcs 
were  performed,  given  only  the  experimental  values  of  lamp  dimensions, 
current  and  pressure.  To  save  time  and  money  only  the  measured  spectral 
radiances  were  used  for  a  tungsten  lamp  and  a  potassium-mercury  lamp. 
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2.1  Theoretical  Results  for  Xenon  and  Krypton 


The  gas  discharge  consists  of  a  cylindrical  arc  confined  by  a 
cylindrical  tube.  In  this  wall-stabilized  arc,  electrical  current 
carries  power  into  the  positive  column  of  the  arc.  Radiation  and  thermal 
conduction  to  the  wall  take  power  out  of  the  positive  column.  Convective 
power  transport  is  negligible  in  wall-stabilized  arcs  and  will  not  be 
considered.  An  arc  discharge  has  a  central  plasma  core  at  the  tempera¬ 
ture  of  the  fully  ionized  plasma  and  a  cold  sheath  adjacent  to  the  wall 
at  a  temperature  determined  by  the  rate  at  which  power  flows  from  the 
wall  to  the  outside  world.  (A  fully  ionized  plasma  is  a  misnomer  which 
means  that  the  electron-electron  and  electron-ion  interactions  are  the 
dominant  particle-particle  interactions.  It  does  not  mean  that  the  gas 
is  completely  ionized.) 

The  emission  from  a  cylindrical  wall-stabilized  gaseous  arc 
discharge  is  calculated  with  a  minimal  amount  of  input  information. 

The  requisite  macroscopic  parameters  are:  (1)  the  length  and  diameter 
of  the  discharge  tube,  (2)  the  electrical  current  into  the  tube,  (3)  the 
tube  wall  temperature,  and  (4)  the  ambient  pressure  in  the  discharge. 
These  parameters  are  easily  obtained  for  dc  discharges.  The  requisite 
microscopic  parameters  are:  (1)  the  amounts  of  each  substance  (e.g., 
water  vapor,  water,  xenon,  etc.)  present  before  the  discharge  is  struck, 
(2)  the  heat  capacities  as  a  function  of  the  temperature  T,  (3)  the 
heats  of  formation  AH(T)  for  each  species  present  in  the  tube  before 
and  during  the  discharge,  (4)  the  lower  lying  energy  levels  and  their 
SLJ  classification  for  each  species  in  the  discharge,  and  (5)  the 
momentum-transfer  cross  sections  for  the  species  in  the  discharge. 

These  parameters  are  difficult  to  obtain  and  are  approximated  whenever 
our  knowledge  is  insufficient. 

Given  the  necessary  input  data,  the  arc  discharge  properties 
are  calculated  as  a  function  of  temperature  in  the  following  order: 

(1)  the  arc  composition,  (2)  the  electrical  and  thermal  conductivities 
of  the  arc  gases,  (3)  the  line  radiation,  (4)  the  continum  radiation,  and 
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(5)  the  line  shifts  and  broadening.  The  temperature  profile  and  the 
resultant  arc  discharge  emission  can  then  be  calculated.  The  calcula¬ 
tion  of  the  line  radiation  by  the  present  computer  programs  restricts 
these  calculations  to  arc  discharges  with  only  rare  and  alkali  metal 
gases.  The  development  of  Phase  Two  of  this  contract  should  extend  the 
applicability  of  these  calculations  to  any  atomic  gases  in  arc 
discharges. 

It  is  easily  seen  that  a  little  input  information  generates  a 
lot  of  output  information.  The  computer  programs  needed  in  the  calcula¬ 
tions  were  developed  by  Westinghouse  on  internal  funds  and  on  previous 
government  contracts.  They  were  updated  and  corrected  whenever  necessary 
during  the  arc  discharge  calculations  which  are  summarized  in  the 
remainder  of  Section  2.1. 

2.1.1  Arc  Composition 

The  equilibrium  composition  of  an  arc  at  a  given  temperature 
and  pressure  is  usually  calculated  by  assuming  local  thermodynamic 

4 

equilibrium  (LTE)  and  then  using  the  Saha  Equations.  These  equations 

need  partition  functions  and  ionization  potentials  for  the  various 

species  in  the  discharge.  However,  in  a  plasma,  charges  of  one  sign 

are,  on  an  average,  surrounded  by  an  excess  charge  of  the  opposite  sign. 

This  lowers  the  ionization  potential  and  Debye  has  derived  the  equation 

4 

which  gives  the  magnitude  of  this  lowering.  In  practice  this  Debye 
lowering  is  introduced  into  the  Saha  Equations  and,  after  an  initial 
guess  of  this  lowering,  the  Saha  Equations  are  iterated  until  a  self- 
consistent  Debye  lowering  is  obtained.  This  method  is  adequate  for 
simple  gases  such  as  xenon  or  krypton  but  totally  inadequate  for 
complicated  discharges. 

R.  Liebermann  has  developed  a  more  complete  approach  which  is 
applicable  to  any  system  (it  has  been  applied  to  the  SF^  discharge  in 
which  each  SF^  molecule  is  decomposed  into  22  different  species).  The 
Gibbs  free  energy  is  calculated  from  given  heat  capacities  and  heats  of 
formation.  The  Gibbs  free  energy  is  then  minimized,^  under  mass  balance 
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and  charge  conservation  constraints,  via  the  geometric  programming 
techniques  of  Duff in  et  al.  The  Debye  lowering  is  again  introduced  as 
an  iteration  parameter  and  the  simultaneous  equations  iterated  until  a 
self-consistent  Debye  lowering  is  obtained. 

For  the  simple  pure  xencn  or  krypton  discharge,  either 
approach  can  be  used  but  because  the  Llebermann  program  is  more  general, 
it  was  used.  Since  the  heat  capacity  should  vary  slowly  with  the  Debye 
lowering,  it  was  held  constant  during  the  iterations  on  the  Debye 
lowering.  The  heats  of  formation  were  used  as  a  function  of  temperature. 
To  simplify  the  calculation  it  was  assumed  that  the  discharge  was 
optically  thick  below  2000  X  and  optically  thin  above  2000  X.  Thus  the 
ultraviolet  radiation  effects  were  completely  ignored.  The  results  of 
the  calculations  for  a  pressure  of  2.5  atm  are  given  in  Figs.  1  and  2. 

The  temperature  range  was  between  1000°K  (the  assumed  wall  temperature) 
and  18,000°K.  We  expected  the  highest  temperature  in  the  discharges  to 
lie  between  6,000°K  and  10,000°K,  but  another  calculation  on  a  Na-Br 
discharge  led  us  to  believe  that  a  temperature  as  high  as  15,000°K  might 
be  possible.  A  comparison  of  Figs.  1  and  2  shows  that  the  higher  first 
and  second  ionization  potentials  of  Kr  require  higher  temperatures  for 
ionization.  In  fact  the  concentration  of  ions  is  not  significant  below 
about  10,000°K  for  Xe  and  below  about  11,0008K  for  Kr. 

2.1.2  Transport  Properties 

The  electrical  and  thermal  conductivities  are  calculated  by 
the  first  Chapman-Enskog  approximations  as  formulated  by  Yos  and  by 

g 

Brokaw.  The  electron-neutral  momentum-transfer  cross  sections  for 

9 

xenon  have  been  measured  by  Frost  and  Phelps  as  a  function  of  tempera¬ 
ture.  To  simplify  the  xenon  calculations  a  value  of  13.17  X^  was  taken 
for  all  temperatures.  This  was  the  value  measured  by  Frost  and  Phelps 
at  10,000°K.  In  the  highly  ionized  regime  Sptizer's  formula  was  used. 
The  neutral-neutral  cross  sections  were  estimated  by  Brokaw' s  method. ^ 
The  neutral-ion  cross  sections  are  unknown  but  based  on  our  previous 
experience  with  other  discharges,  a  value  of  2.5  times  the  neutral- 
neutral  cross  sections  seemed  reasonable  and  was  taken.  Since  the  ions 
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should  constitute  a  small  part  of  the  discharge,  this  approximation  Is 

not  critical  and  may  even  be  Insignificant.  It  also  seemed  reasonable 

(2  2) 

to  Identically  equate  the  unknown  Q  *  averaged  cross  sections  to  the 
averaged  cross  sections. ^  The  xenon  conductivity  results  are 
shown  In  Figs.  3  and  4. 

Included  In  Figs.  3  and  4  are  the  results  of  Devoto  which 
12  13 

were  published  in  1967  and  1969.  When  the  present  calculations 
were  done  only  the  1967  Devoto  results  were  known,  and  they  were  thought 
to  be  In  error.  The  1969  Devoto  results  confirm  this  and  are  quite 
close  to  the  present  results;  the  differences  are  easily  explained.  At 
low  temperatures  Devoto  used  a  temperature-dependent  electron- neutral 
cross  section  whereas  our  value  was  temperature-independent.  At  high 
temperatures  Devoto  used  an  ion-neutral  cross  section  which  was  about 
10  times  the  neutral -neutral  value  whereas  we  used  a  multiplicative 
factor  of  2.5. 

The  electrical  and  thermal  conductivities  were  more  accurately 

calculated  for  the  krypton  discharge.  The  temperature  dependent 

electron- neutral  momentum-transfer  cross  section  was  used  and,  as  given 
13 

by  Devoto,  is  shown  in  Fig.  5.  As  a  check  on  our  cross  sections  and 

calculations,  the  electrical  conductivity  was  calculated  at  1.0  atmo- 

13 

spheres  and  compared  with  Devoto' s  results  (Devoto  did  the  calculation 
at  1.0  atmospheres, but  not  at  2.5  atmospheres).  Our  results,  labeled 
Liebermann,  are  compared  to  Devoto' s  results  in  Fig.  6.  At  high  tempera¬ 
tures  the  electron-ion  momentum-transfer  cross  sections  are  significant. 
The  curves  differ  in  these  regions  because  Devoto  used  different  values 

for  these  cross  sections.  At  low  temperatures  the  fourth-order  calcula- 
14 

tion  by  Devoto  deviates  from  our  first-order  results  whereas  between 

7,000  and  8,000°K  the  first-  and  fourth-order  theories  agree.  Devoto 

14 

found  a  similar  relationship  for  an  argon  arc  discharge. 

The  conductivity  calculations  for  krypton  were  completed  by 
equating  Q^2’2^  to  Q^1’1^  for  the  neutral-neutral  and  the  electron- 
neutral  interactions.  The  neutral-neutral  cross  sections  were  again 


calculated  by  Brokaw's  method. ^  The  Q^1’1^  neutral-ion  cross  sections 

were,  estimated  as  twice  the  neutral-neutral  cross  sections,  whereas  the 

(2  2) 

Q  ’  neutral-icn  cross  sections  were  estimated  as  ten  times  the  neutral- 
neutral  cross  sections.  The  conductivity  results  are  given  in  Figs.  3 
and  4.  These  curves  show  a  surprisingly  small  difference  between  the 
xenon  and  krypton  results,  particularly  for  the  electrical  conductivity. 

In  fact,  although  small,  the  different  cross  section  approximations  used 
by  Devoto  give  larger  differences  than  changing  the  gas  from  xenon  to 
krypton. 

2.1.3  Line  Radiation 

Much  of  the  discharge  radiation  originates  from  bound-bound 
transitions,  l.e.,  line  radiation.  Although  much  of  the  science  is 
known  for  these  transitions,  the  problem  is  most  assuredly  nontrivial. 

Two  simplifying  assumptions  which  are  usually  made  is  that  there  is  no 
configuration  interaction  and  that  magnetic  dipole  radiation  is  insig¬ 
nificant  compared  to  the  electric  dlpcle  transitions.  The  first 
approximation  should  be  very  good  for  simple  gases  such  as  xenon  but  is 
very  bad  for  certain  gases,  e.g.,  mercury.  The  second  approximation 
should  be  excellent  except  for  the  -  £n  transitions  in  the  rare 
earths.  Both  approximations  are  made  iTor  xenon  and  krypton. 

The  wave  finction,  which  determines  the  electric  dipole 

transition  probability  rate,  is  determined  from  a  SLJ  listing  of  the 

energy  levels.  For  xenon  I  only  the  lowest  five  configurations  are 

completely  known;^^  these  known  levels  are  listed  in  Table  I.  For 

krypton  1  the  lowest  nine  configurations  are  completely  known.  By 

3  1 

setting  the  unknown  Pq  energy  level  equal  to  the  energy  level 
within  each  ns  configuration  for  n  between  eight  and  twelve,  five  high 
energy  ns  configurations  can  also  be  included.  The  energy  levels  for 
these  fourteen  configurations  are  listed  in  Table  II. 

With  this  energy  level  information  the  theoretical  expressions 
for  the  energy  levels  can  be  adjusted  to  give  a  least  squares  fit  to  the 
experimental  levels  by  using  the  radial  integrals  as  fitting  parameters. 
The  resultant  intermediate  coupling  angular  wave  functions  for  xenon  are 
given  in  Table  III.  For  example,  in  the  5p^  5d  configurations: 
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|1P1  >  -  0.717 1  3D1>  -  0.2981^)  +  0.6301^)  (1) 

where  |  >  Is  the  intermediate  coupling  and  the  |)  Is  the  SL  coupling 
wave  function.  The  Intermediate  coupling  angular  wave  functions  for 
krypton  are  given  In  Table  IV  (Tables  III  and  IV  do  not  list  the  trivial 
one-component  eigenvectors).  These  last  two  tables  dramatically  show  the 
well  known  collapse  of  SL  coupling  in  the  rare  gases. 

The  Westinghouse  Intermediate  coupling  program  Is  only  valid 
for  the  rare  gases,  and  moreover,  was  found  to  give  incorrect  answers. 
Previous  to  this  contract  the  very  versatile  electric  dipole  transition 
computer  program  written  by  Cowan1^  had  been  obtained  and  converted  from 
the  IBM  STRETCH  Computer  to  the  CDC  6600  Computer.  A  comparison  of  the 
two  computer  results  for  xenon  located  the  Westinghouse  program  error 
which  was  then  circumvented.  (If  the  Cowan  computer  program  were  to  be 
Incorporated  into  the  Westinghouse  intermediate  coupling  program,  any 
atom  with  two  open  shells  or  less  could  be  handled.)  The  resultant 
radial  integral  values  for  xenon  and  krypton  are  listed  in  Table  V.  The 
n'd  spin-orbit  coupling  parameters  must  be  positive.  The  two  small 
negative  values  are  therefore  unrealistic  and  are  a  result  of  cumulative 
errors  or  configuration  interaction.  In  any  event  these  values  are  too 
small  to  significantly  affect  the  results. 

The  radial  part  of  the  wave  functions  must  also  he  known  to 

2 

calculate  the  value  of  I  IrU  >1  which  is  needed  to  calculate  the 

transition  probability  rate.  These  radial  wave  functions  are  very 

difficult  to  calculate.  For  simple  atoms  (less  than  two  equivalent 

electrons)  the  radial  wave  functions  are  calculated  from  the  Bates  and 

18 

Damgaard  Coulomb  approximation  or  the  Thomas-Ferml  approximation, 
whichever  is  most  appropriate.  Both  of  these  methods  require  one  electron 
energies  and  consequently  break  down  whenever  more  than  one  equivalent 
electron  is  present.  The  development  of  Phase  2  of  this  contract  will 
remove  this  restriction  and  also  make  this  part  of  the  calculation  con¬ 
sistent  with  the  intermediate  coupling  determined  radial  integrals. 
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The  Thomas-Fermi  method  was  exclusively  used  for  xenon  and 

for  the  five  lowest  configurations  of  krypton.  The  more  Important  values 

of  the  resultant  transition  probability  rates  and  transition  wavelengths 

are  listed  in  Table  VI.  As  a  check  of  our  results,  the  thirty  5s-5p 

transition  probabilities  for  krypton  were  compared  to  the  results  of 
19 

Murphy.'  The  agreement  is  within  8%  which  is  very  satisfactory.  The 
results  of  Table  VI  show  that  both  atoms  emit  considerable  energy  in 
the  ultraviolet  although  xenon  radiates  about  twice  as  much  energy  in 
tills  wavelength  region  as  does  krypton.  Consequently,  the  ultraviolet 
optical  thickness  assumption  made  in  the  arc  discharge  calculations 
will  be  a  worse  approximation  for  xenon  than  for  krypton.  In  the 
infrared  region  xenon  essentially  has  no  radiation  whereas  krypton  does 
have  a  small  amount  of  radiation.  The  total  line  radiation  as  a 
function  of  temperature  has  been  calculated  and  is  given  in  the  next 
subsection  where  it  :.s  compared  to  the  continuum  radiation. 

2.1.4  Continuun  Radiation 

The  free-bound  and  free-free  electronic  transitions  yield 
continuum  radiation,  out  the  theory  for  these  transitions  is  rather 
crude.  The  spectral  absorptivity  for  these  transitions  in  given  by 

K(A,T)  -  1  °i  <*>  e*P  <-E±/kT)/Zo  (2) 

where  Zq  Is  the  internal  partition  function,  g^  is  the  statistical  weight 

of  the  initial  state,  is  the  energy,  and  is  the  photoelectric  cross 

section  for  the  ith  species.  The  exponential  factor  is  simply  the 

20 

Boltzmann  factor.  Unsdld  showed  for  hydrogen  how  the  summation  could 

2L 

be  changed  to  an  integration.  SchlUter  has  generalized  Unsbld's  work 
to  give 

2  6 

K  *  3,2Tr  \k4  z~  *3  exP  ("X/kT)  [exp(hc/XkT)-l]yC  (3) 

3/J  he  o 

where  x  Is  the  ionization  energy,  y  is  the  statistical  weight  of  the  core 
state,  and  S  is  a  temperature  dependent  factor  with  0(1).  The  latter 


* 
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factor  Is  very  difficult  to  determine.  The  values  given  in  SchlUter 's 

paper  for  xenon  are  reproduced  in  Fig.  7.  The  Schltiter  curves  should  be 

22 

the  most  accurate.  Blberman  et  al.  used  the  hydrogenic  approximation 
23 

whereas  Yankov  used  a  modified  hydrogenic  approximation.  For  short 
wavelengths  £  Increases  exponentially  with  increasing  wavelength.  In 

O 

the  Intermediate  range  between  5000  and  6500  A,  the  £  factor  is  extremely 
difficult  to  determine.  In  this  region  the  absorption  edges  appear,  but 
they  are  shifted  to  longer  wavelengths  by  plasma  interactions  and  there 
is  no  satisfactory  theory  for  these  shifts.  Hence,  the  £  factor  is 

O 

extremely  crude  between  5000  and  6500  A.  The  absorption  edges  should  be 
sharper  than  those  shown  by  SchlUter,  and  might  appear  as  Llebermann's 
curve  shows  in  Fig.  7.  The  £  factor  for  krypton  is  given  in  Fig.  8. 

The  evaluation  of  Eq.  (3)  in  the  literature  has  led  to  some 

confusion.  The  value  of  y  for  xenon  should  technically  be  four,  since 

2  2 

the  core  state  is  split  into  the  2^/2  and  Pl/2  8tates  and  tlie  J  “  3/2 
state  is  the  ground  state.  However,  a  few  people  such  as  SchlUter  use 
a  value  of  six.  SchlUter  has  shown  that  the  value  of  y  is  not  Important 
if  one  is  consistent  because  the  £  factor  contains  a  y  *  factor.  However, 
we  used  the  partition  function  for  y  (which  is  approximately  four)  but 
used  SchlUter 's  £  which  contained  the  value  of  y  ■  6.  This  confusion 
lead  to  a  few  Incorrect  early  results  but  was  eventually  corrected. 

The  continuum  radiation  is  calculated  from  the  effective 
spectral  absorptivity 


K'  -  K[1  -  exp (-hc/XkT) ]  (4) 

which  is  the  difference  between  the  true  absorption  and  Induced  emission. 
Equation  (4)  is  also  valid  for  the  line  absorptivities .  The  radiation  is 
then  calculated  as  the  radiated  power  per  unit  volume,  as  is  given  by  the 
spectral  emission  coefficient 

e(X,T)  -  BXK'  (5) 
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where  is  th-1  Planck  black  bod;  Intensity  which  is  given  by 

=  2hc/\^[exp (hc/XkT)-l]  (6) 

The  line  and  continuum  radiation  results  are  given  in  Figs.  9 
and  10.*  The  former  figure  shows  that  the  line  radiation  is  more  important 
in  xenon  than  in  krypton  and  is  always  significant  for  both  arc  discharges 
within  the  specified  temperature  range.  The  latter  figure  shows  that 
below  12,000°K  the  xenon  arc  discharge  emits  more  radiation  than  the 
krypton  arc  discharge.  The  radiation  peaks  in  Fig.  10  are  a  result  of 
two  competing  factors;  as  the  temperature  is  increased  each  atom  increases 
its  radiation  but,  at  constant  pressure,  the  atomic  density  decreases. 

This  situation  is  in  contrast  to  the  normal  experimental  situation  where 
the  atomic  density  is  held  constant  and  the  pressure  is  allowed  to  change. 

2.1.5  Temperature  Profile 
24 

Lowke,  at  Westinghouse,  has  developed  a  very  powerful  computer 
program  for  the  determination  of  temperature  profiles.  Given  the  tube 
radius,  the  wall  tenperature  and  the  central  temperature  of  the  arc,  the 
program  uses  a  relaxation  method  to  determine  a  profile,  the  current  and 
the  electric  field  strength.  Of  cours.*,  the  electrical  und  thermal  con¬ 
ductivities  and  the  radiation  must  also  be  given.  The  electric  field 
strength  is  guessed  and  then  iterated  until  a  self-consistent  solution 
is  obtained.  The  discharge  is  assumed  to  be  optically  thin. 

The  electrical  characteristics  of  the  xenon  and  krypton  arc 

discharges  are  shown  In  Fig.  11.  Two  characteristics  are  shown  for 

12 

xenon;  one  was  determined  from  Devoto's  conductivities  and  one  was 
determined  from  our  calculated  conductivities  (see  Figs.  3  and  4).  The 
measured  values  are  also  given.  At  the  measured  current  of  forty  amperes, 
the  electric  field  strength  calculated  from  our  conductivities  is  4%  too 


The  xenon  continuum  results  in  this  report  differ  from  those  given  in 
Technical  Report  No.  1.  Subsequent  to  the  printing  of  the  latter  report 
a  keypunching  error  in  the  continuun  computer  program  was  discovered 
and  corrected. 
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low  for  xenon  and  8%  too  high  for  krypton.  This  agreement  between  the 

measured  and  calculated  values  of  the  electric  field  strengths  is  very 

25 

good,  particularly  since,  as  Mottschmann  explains,  the  ultraviolet 
radiation  which  was  ignored  should  lower  the  calculated  value. 

A  qualitative  dependence  of  the  electrical  characteristic  upon 
the  thermal  and  electrical  conductivities,  and  thus  on  the  cross  sections, 
is  obtained  by  comparing  the  two  xenon  characteristics.  Figures  3  and  4 
show  that  Devoto's  conductivities  are  higher  than  ours  and  Fig.  11  shows 
that  as  a  result  of  Ohm's  law  these  higher  conductivities  yield  a  lower 
electric  field  strength  (40%  too  low  at  40  amp)  which  in  turn  leads  to 
a  lower  central  core  temperature  of  the  arc  discharge. 

The  calculated  temperature  profiles  for  a  central  core 
temperature  of  9200°K  (xenon)  and  9775#K  (krypton)  are  given  in  Fig.  12. 

They  are  quite  flat  over  most  of  the  tube  volume.  Every  other  point  is 
calculated  and  the  intervening  points  are  interpolated.  Linear  inter¬ 
polation  explains  the  small  dip  at  about  2  mm.  The  wall  temperature  was 
fixed  at  1000°K  for  the  calculation,  although  the  lamp  was  water  cooled. 

The  steep  slope  near  the  wall  shows  that  a  crude  approximation  of  wall 
temperature  is  sufficiently  accurate. 

2.1.6  Arc  Emission 

The  wavelength  shifts  and  broadening  of  the  line  radiation  is 

4 

first  calculated.  Griem's  electron  impact  mechanism  is  evaluated  as 
well  as  the  van  der  Waals  interaction  between  different  species. 

Resonance  broadening  is  calculated  for  all  transitions  in  which  either 
the  lower  or  upper  state  has  an  allowed  electric  dipole  transition  to 
the  ground  state.  Whenever  ions  are  present  the  charges  cause  a 
polarization  effect  which  also  shifts  the  lines  and  which  is  also 
calculated.  The  broadening  varied  from  zero  to  two  or  three  angstroms, 
and  a  typical  wavelength  shift  was  one  angstrom. 

These  broadenings  and  shifts  are  calculated  for  each  line  at 
six  different  temperatures  (see  Fig.  12).  The  absorptivities  of  these 
lines  are  calculated  and  extended  by  interpolation.  The  line  radiation 

}' 
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is  then  integrates  along  a  diameter  to  give  an  "effective"  line  which 
is  the  actual  emission  observed  by  the  experimentalist.  The  integration 
is  also  performed  for  the  continuum  radiation.  The  xenon  and  krypton 
results  are  shown  with  the  experimental  results  in  the  next  section. 

2 . 2  Comparts on  of  Theore tical  and  Experimental  Results 

The  lamps  used  for  spectral  measurements  were  built  at  the 

Hcsuat'ch  Laboratories.  They  have  been  already  tested  as  laser  pumps, 

and  spectroscopic  measurements  were  previously  made  under  slightly 

26 

Jifferent  operating  conditions.  An  important  parameter  in  all  calcu¬ 
lations  is  atom  density  or  pressure.  In  sealed  lamps,  the  fill  pressure 
is  known,  but  the  operating  pressure  changes  with  power  and  can  only  be 
calculated  if  the  inactive  or  dead  volume  of  the  lamp  is  known.  To 
eliminate  this  possible  source  of  error,  the  lamps  were  not  sealed  but 
were  connected  to  a  large  reservoir  of  fill  gas  kept  at  2.5  atm.  Other 
measured  parameters  vere  arc  diameter  «’  0.56  cm,  arc  length  ■  7.3  cm, 
current  «*  40  A,  and  voltage  ■  82  V  for  xenon  and  69  V  for  krypton. 

Another  lamp  of  14.1  cm  arc  length  was  operated  under  ideitical  conditions 
in  order  to  determine  the  electrode  loss  which  then  was  calculated  to  be 
12  V.  Therefore,  the  electric  field  in  the  positive  column  is  9.6  V/cm 

for  xenon  and  7.8  V/cm  for  krypton,  and  the  electrical  input  power 

2 

densities  (power  per  unit  surface  area)  are  226  and  192  watts/cm 
respectively. 

Spectral  measurements  were  performed  using  a  Jarrell-Ash  1/2  m 
Ebert  monochromator,  for  measurement  of  the  line  spectra  the  curved 

O 

slits  were  40  y  wide  corresponding  to  a  dispersion  of  0.64  A.  In 
measuring  the  continuum  the  slits  were  opened  to  100  y  in  order  to 
increase  the  detected  signal  of  the  relatively  weak  continuum  radiance. 

The  detector  used  was  a  silicon  p-i-n  phetodiode.  The  system  was  cali¬ 
brated  using  a  calibrated  tungsten  spectral  radiance  standard  and  auto- 

O 

mat ically  corrected  for  nonlinearity  between  5000  -  10,000  A  by  using  a 
variable  attenuator. 
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2.2.1  Xenon 

The  experimental  results  are  shown  and  compared  with  the 
computer  calculations  in  Fig.  13.  Figure  13a  shown  the  experimentally 

O 

obtained  radiance  spectrum  of  xenon  between  7500  -  10,000  A.  The 
computerized  equivalent  spectrum  is  shown  in  Fig.  13b.  The  two  spectra 
compare  well  in  the  location  of  the  lines.  Three  lines  around  9500  X 
which  appear  in  the  experimental  spectrum  do  not  appear  in  the  calculated 
spectrum.  The  reason  for  this  discrepancy  is  not  certain,  but  it  appears 
that  these  lines  originate  from  high  lying  configurations  of  xenon  which 
were  net.  included  in  the  calculation. 

Other  differences  in  data  include  the  relative  amplitudes  of 
the  lines.  By  examining  the  experimental  data,  it  appears  that  all  the 
lines  have  about  the  same  width.  The  narrower  lines  have  greater  intensity 
relative  to  the  measured  spectrum.  The  wider  the  line  the  closer  is  the 
agreement.  This  discrepancy  could  be  a  result  of  inadequate  spectral 
resolution  of  the  monochromator. 

When  taking  the  experimental  spectrum  we  also  measured  in 

O 

detail  the  line  shape  of  the  most  intense  line  occurring  at  about  8232  A. 
This  measurement  is  compared  with  the  calculated  results  at  two  different 
core  temperatures  in  Fig.  14.  The  ordinate  and  absissa  are  the  same  for 
both  curves.  A  comparison  of  the  experimental  line  and  the  8450°K 

®  o 

computed  line  shows  that  the  experimental  half  width  of  3.66  A  is  1.0  A 
greater  than  the  calculated  half  width.  The  dispersion  limited  resolu- 

O 

tion  of  the  monochromator  is  0.64  A  but  the  aberration  in  the  instrument 

O 

could  decrease  the  true  resolution  to  about  1  A  although  the  manufacturer 

O 

places  the  aberration  resolution  at  0.1  A.  Regardless  of  the  resolution, 
the  energy  in  the  line,  which  is  proportional  to  the  area  under  the  line, 
i.e.,  the  product  of  radiance  and  angstroms,  remains  constant.  Comparing 
the  energies  of  the  calculated  and  experimental  results,  we  find  that  the 
calculated  energy  is  15%  smaller.  This  difference  is  surprisingly  small 
since  neither  spectral  curve  was  expected  to  be  accurate  to  within  20%. 
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A  c'-npirison  of  the  experimental  line  and  the  9200°K  computed 

O 

line  shows  that  the  experimental  half  width  is  1.48  A  smaller  than  the 
computed  half  width,  and  that  the  calculated  energy  is  125%  larger. 

These  results  indicate  that  a  central  core  temperature  of  8450  °K  is 
slightly  too  low  whereas  9200°K  is  considerably  too  high.  These  line 
shape  measurements  indicate  a  central  core  temperature  of  8530  +  1508K. 

The  calculated  curves  also  show  some  fine  structure  not 
present  in  the  experimental  spectrum  because  of  limited  resolution.  The 
notch  in  the  short  wavelength  side  of  the  peak  is  real  and  due  to  self 
absorption  of  the  hot  core  by  the  cool  gas  near  the  wall  which  has  its 
peak  absorptivity  at  shorter  wavelengths  than  the  core.  There  is  also 

O 

a  difference  in  the  location  of  the  peak  radiance  of  about  1.3  A. 
However,  neither  value  is  sufficiently  accurate  to  determine  which  curve 
is  in  error. 

In  addition  to  calculating  t:he  spectral  radiance,  the  computer 

also  calculates  the  peak  absorptivity,  and  the  emittance.  These  ealeu- 

0 

lations  for  the  8232  A  line  are  shown  in  Figs.  15  and  16.  The  emittance 

is  the  information  needed  to  calculate  the  spectral  overlap.  From  the 

spectral  absorptivity  we  can  determine  the  degree  of  optical  thickness 

0 

of  the  line.  In  thin  way  we  have  deteimined  that  the  8252  A  line  does 
not  have  the  strongest  core  radiance,  but  it  is  less  absorbed  by  the 
cooler  gas  so  that  its  emitted  radiance  is  greatest.  However,  its 
radiation  is  also  appreciably  absorbed;  so  that  the  line  radiation  from 
a  relatively  small  diameter  low-current  discharge  is  already  becoming 
optically  thick.  In  fact  at  9200°K  the  discharge  is  80%  optically  thick 
and  below  6000°K  the  discharge  is  100%  optically  thin  (based  on  a 
maximum  spectral  absorptivity  within  the  specified  wavelength  region  of 
10  cm  *  as  being  100%  optically  thick). 

The  experimentally  measured  continuum  spectrum  is  shown  in 
Fig.  17  and  a  calculated  continuum  spectrum  in  Fig.  18.  Because  of  a 
programming  error  the  ordinate  was  corrected  as  shown  in  the  curve. 

O 

Beyond  7000  A,  comparison  becomes  difficult  because  of  the  contribution 
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of  Che  wings  of  Che  lines  Co  Che  conCinuum.  A  quick  glance  aC  Che  Cwo 
curves  shows  ChaC  Che  agreemenC  is  poor  since  Che  shapes  of  Che  curves 

O 

below  6000  A  are  differenC.  (A  keypunching  error  rendered  impossible 

any  absoluCe  comparison  of  conCinuum  radiaCion  ampliCudes  for  Chis  figure 

alchough  Che  line  radiaCion  resulCs  are  correcC.)  We  also  Cried  Che  less 

23 

sophisCicaCed  approach  of  Yankov.  The  resulC  appeared  Co  be  jusC 
slighcly  beCCer  buC  Che  shape  of  Che  curve  sCill  disagreed  wich  Che 
measured  curve. 

As  explained  in  Che  secCion  describing  Che  meChod  of  calculaCing 

O 

Che  conCinuum,  Che  bound-free  absorpCion  edges  occur  in  Che  5000  -  6500  A 
region.  Therefore,  Che  shape  of  Che  conCinuum  depends  scrongly  on 
correccly  chosing  SchlliCer's  £  funcCion.  SchlliCer  could  noc  himself 
gee  good  agreemenC  wich  experimenc  for  xenon.  We  considered  Che  possi- 
bilicy  of  ionic  species  causing  Che  difficulCy,  buC  our  calculaCed  resulCs 
show  chac  Chis  is  noC  a  significanC  facCor.  We  also  modified  SchlliCer 's 
£  funcCion  as  shown  by  Che  doCCed  line  in  Fig.  7  Co  accounC  for  CemperaCure 
shifCs  and  obCained  Che  specCrum  shown  in  Fig.  19.  The  change  in  speccrum 
is  considerable,  and  Che  shape  above  5100  X  is  in  much  beCCer  agreemenC 
wich  experimenc.  A  comparison  of  Che  measured  and  calculaCed  conCinuum 

O 

ampliCudes  shows  fair  agreemenC  below  6000  A  buC  Che  calculaCed  ampliCude 

O 

is  more  Chan  a  facCor  of  Cwo  Coo  low  above  6000  A.  Our  purpose  is  noC 
Co  find  a  beCCer  SchlliCer  facCor  by  fleeing  iC  Co  Che  measured  speccrum 
buC  racher  Co  show  chac  a  much  beCCer  knowledge  of  Che  SchlliCer  facCor 
is  necessary  if  beCCer  conCinua  are  Co  be  calculaCed.  In  Che  remaining 
calculaCions  Che  calculaCed  conCinuum  of  Fig.  19  is  used. 

2.2.2  KrypCon 

The  experimenCal  resulCs  are  shown  and  compared  wich  Che 
compuCer  calculaCions  in  Fig.  20.  Figure  20a  shows  Che  experimenCally 

O 

obCained  specCral  radiance  of  krypCon  beCween  7500  -  10,000  A.  The 
compuCerized  equivalenC  speccrum  is  shown  in  Fig.  20b.  A  comparison 
of  Che  line  posiCions  shows  excellenC  agreemenC  beCween  Che  Cwo  specCra. 

As  explained  in  Che  previous  subsecClon  Che  ampliCudes  of  Che  Cwo  specCra 
do  noC  compare  as  well  because  Che  compuCer  resoluCion  is  finer  Chan  Che 
experimenCal  resoluCion. 


20 


The  experimentally  measured  7,606  A  line  and  the  same  line 
computed  with  a  core  temperature  of  9,775°K  are  compared  in  Fig.  21. 

O 

The  computed  line  exhibits  a  -5  A  shift  relative  to  the  measured  line. 

15  6 

Since  the  unshi:ted  line  should  occur  at  7,604  A,  the  experimental 

O 

apparatus  is  in  error  by  at  least  2  A.  The  remaining  discrepancy  is 
small  but  may  indicate  a  slightly  too  high  core  temperature.  A  com¬ 
parison  of  the  energy  under  each  line  shows  that  the  computed  line  has 
21%  more  energy.  This  very  good  agreement  is  well  within  the  experi¬ 
mental  and  computed  accuracies  but  again  Indicates  that  the  computed 
core  temperature  may  be  just  slightly  too  high. 

As  for  the  xenon  discharge  the  emlttance  and  spectral 
absorptivlties  were  also  calculated.  The  krypton  discharge  is  23% 
optically  thick  at  9775°K  and  100%  optically  thin  below  7,000°K. 

The  experimentally  measured  continuum  spectrum  is  shown  in 
Fig.  22  while  the  calculated  continuum  spectrum  is  shovm  in  Fig.  23. 

As  explained  in  the  previous  subsection  the  continuum  calculation  is 
crude.  The  Schlllter  factor  shown  in  Fig.  8  was  used  although  it  again 

O 

lead  to  the  wrong  curve  shape  below  6,000  A.  The  amplitudes  agree 

O 

fairly  well  below  6.000  A,  but  the  calculated  amplitude  is  about  a 

O 

factor  of  two,  too  low  above  6,000  A.  These  figures  again  show  that  a 
better  knowledge  of  the  SchlUter  factor  is  necessary  if  better  con¬ 
tinuums  are  to  be  calculated. 

2 . 3  Measured  Spectra  for  W  and  K-Hg  Lamps 

The  spectral  radiance  of  the  commonly  used  tungsten  lamp  was 

obtained  by  multiplying  the  black  body  radiance  at  3,200°K  by  the 
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measured  emissivity  at  3,200°K.  The  spectral  emittance  was  obtained 
by  multiplying  by  n  steradians,  This  multiplication  factor  is  correct 
for  an  optically  thick  infinite  plane  emitter  radiating  onto  an  infinite 
plane  absorber  of  finite  thickness,  or  for  an  infinitely  long  cylinder 
of  finite  radius  radiating  onto  an  infinitely  long  concentric  annulus 
where  the  annulus  thickness  is  much  smaller  than  the  annulus  radius . 
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The  input  power  to  the  tungsten  lamp  was  calculated  as  eoT  where  e 
is  the  average  emissivity  at  T  =  3,200°K  and  a  is  the  Stefan-Boltzmann 
constant . 

The  potassium-mercury  spectral  emittance  could  be  theoretically 
calculated  as  in  the  last  section  hut  to  save  time  and  money  the 
experimental  spectral  radiance,  shown  in  Fig.  24,  was  used  since  it 
was  already  available.  The  spectral  emittance  was  obtained  by 
multiplying  the  spectral  radiance  by  three.  The  factor  of  three  was 
obtained  by  averaging  the  factors  calculated  for  the  xenon  arc  dis¬ 
charge  which  should  be  about  as  optically  thick  as  the  potassium-mercury 
discharge.  The  input  power  density  is  given  by  the  100  volt  voltage 
drop  multiplied  by  the  5  amp  current  divided  by  the  2.36  cm  circum¬ 
ference  times  the  3.5  cm  arc  length  .  This  answer  is  then  reduced  by 
15%  to  exclude  the  end  losses. 

2.4  Discussion 

Because  of  the  narrowness  of  the  lines,  it  takes  about  80 

0 

minutes  to  make  one  monochromator  scan  from  5,000  -  10,000  A,  but  the 
major  time  consumption  is  spent  in  building  the  lamp, set-up  time, 
erroneous  starts,  etc.  The  entire  procedure  would  have  to  be  repeated 
for  every  parameter  change  such  as  current  density,  pressure,  and  wall 
diameter.  In  addition,  after  obtaining  the  data  it  still  has  to  be 
converted  to  computer  data  cards.  This  spectrum  provides  the  spectral 
radiance  through  the  center  of  the  arc.  To  get  the  spectral  emittance, 
i.e.,  the  radiation  emitted  at  the  wall  from  all  portions  of  the  arc, 
one  would  have  to  obtain  the  radiance  from  various  rays  through  the 
arc  and  calculate  the  emittance.  The  above  discussion  points  out  the 
difficulty  of  doing  the  spectral  overlap  calculation  from  experimental 
data.  Furthermore,  if  a  pulsed  arc  were  studied  the  experimental  degree 
of  difficulty  would  be  further  increased. 

The  object  of  comparing  the  Xe  and  Kr  calculated  spectra  with 
the  experimental  spectra  was  to  determine  how  close  the  calculations 
match  the  experiments.  The  results  indicate  that  the  line  spectra 
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calculations  are  appreciably  better  than  the  experimental  data.  There 
is  no  absolute  way  of  knowing  which  spectrum  is  correct,  but  it  is 
certain  that  the  experimental  data  is  limited  in  resolution  to  at  least 

O 

0.7!)  A.  Since  the  measured  lines  are  only  a  few  angstroms  wide,  a 

O 

resolution  of  0.7f>  A  will  cause  noticeable  error.  However,  the  measured 
energy  within  a  line  is  independent  of  the  resolution.  The  agreement 
for  krypton  to  within  21%  between  calculation  and  experiment  is  better 
than  the  accuracy  of  each  method  and  therefore  impossible  to  determine 
which  method  is  nearer  reality.  The  xenon  line  profile  at  different 
core  temperatures  shows  how  accurate  line  profiles  can  be  compared  to 
determine  the  actual  core  temperature.  We  therefore  conclude  that  it 
is  very  desirable  to  calculate  the  line  spectrum.  It  is  stressed  that 
the  measurements  and  calculations  were  carried  out  completely 
independently.  The  calculated  line  spectra  were  never  altered  to  fit 
the  experimental  spectra.  If  the  experimental  measurements  were  never 
made,  the  identical  calculated  spectra  would  have  been  presented. 

The  calculated  continua  spectra  were  only  about  a  factor  of 
two  in  agreement  with  experiment.  The  difference  was  net  only  in 
magnitude,  but  in  wavelength  dependence.  Since  the  experimental  radiance 
is  estimated  to  have  an  absolute  certainty  of  ±20%  and  a  relative 
accuracy  with  wavelength  of  better  than  5%,  it  is  believed  that  the 
calculated  continua  are  in  error.  We  believe  that  the  problem  is 
associated  with  the  £  factor  of  SchlUter  although  attempts  at  modifying 
it,  after  making  spectral  comparisons,  were  not  entirely  successful 
even  though  it  was  apparent  that  this  factor  can  exert  a  strong  influence. 
As  shown  in  Figs.  7  and  8,  the  SchlUter  factor  appears  to  be  too  small 

O 

above  6,000  A. 

The  true  value  for  the  continuum  radiation  from  the  xenon 
discharge  can  be  crudely  approximated  by  forming  an  energy  balance  of 
the  arc: 
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where  r  is  Che  discharge  Cube  radius.  The  lefC  hand  side  is  Che 
elecCrical  inpuC  power  per  cm  and  Che  righc  hand  side  is  Che  oucput 
power  per  cm  as  deCermined  by  Che  Chermal  conducCion  and  Che  line  and 
conCinuum  radiacion,  given  by  R^  and  R£,  respecCively  (see  Fig.  10). 

By  assuming  a  nearly  isochermal,  opCically  Chin  discharge,  Che  lefC  hand 
side  becomes 


2  2 

ttE  or  (8) 

The  CemperaCure  gradienC  dT/dr  can  now  be  evaluaCed  by  using  Che  cal- 
culaCed  value  of  9.2  volCs/cm  for  E  and  using  Figs.  3,  4,  and  10  Co 
obcain  Che  ocher  unknown  values  aC  9,200°K.  This  value  of  Che  Cempera¬ 
Cure  gradienC  is  now  assumed  consCanC  and  Che  conducCiviCies  are  used 
for  8,530°K  (Che  core  CemperaCure  deCermined  by  a  comparison  of  Che 
calrulaCed  line  profiles  Co  Che  measured  profile) .  The  value  of  R^  aC 

8,530°K  is  deCermined  from  Figs.  9  and  10,  and  Che  value  of  R  is 

C 

deCermined  from  Eqs.  7  and  8.  In  ocher  words,  Che  elecCrical  inpuC 
power  is  held  consCanC  buC  Che  CemperaCure  of  Che  discharge  is  lowered 
from  9,200°K  Co  8,530°K.  This  CemperaCure  decrease  lowers  Che  con¬ 
ducCion  and  line  radiacion  losses  and  consequenCly  Che  Crue  conCinuum 
radlaClon  loss  musC  Increase  Co  exacCly  cancel  Chese  decreased  losses. 

The  resulC  is  ChaC  Che  conCinuum  radiacion  aC  8,530°K  is  approximacely 
-3  -1 

99.4  waCCs  cm  sCer  which  is  nine  Cimes  as  much  conCinuum  radiacion 

as  given  by  Figs.  9  and  10,  alchough  Che  CoCal  radiacion  has  only  been 

-3  -1 

reduced  from  165  Co  155  waCCs  cm  sCer  .  Thus  Che  conCinuum  radiacion 
aC  8,530°K  is  64%  of  Che  CoCal  radiacion.  This  compares  favorably  wich 
Che  55%  calculaCed  for  Che  krypCon  discharge  and  leads  us  Co  believe 
ChaC  Che  conCinuum  radiacion  calculaCed  for  Che  krypCon  discharge  is 
abouC  righc. 

In  Che  deCerminaCion  of  Che  CemperaCure  profiles  and  elecCrical 
characCerisCics,  Che  discharge  is  assumed  Co  be  opCically  Chin.  The 
effecC  of  Chis  assumpCion  is  Co  increase  Che  volCage  gradienC  value  aC 
a  given  currenC.  This  increase  can  be  seen  in  Che  xenon  discharge  by 
assuming  an  isochermal  CemperaCure  profile.  Then 
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where  E  and  I  are  given  in  Fig.  11  and  a  is  the  electrical  conductivity 
given  in  Fig.  3.  The  constant  K  is  determined  to  be  6.44  by  setting 
I  equal  to  40  amp  and  using  the  E  and  a  values  appropriate  for  9,200°K. 

If  the  electrical  conductivity  is  now  changed  to  the  value  for  8,5306K 
(the  core  temperature  determined  by  a  comparison  of  the  calculated  line 
profiles  to  the  measured  profile),  the  voltage  gradient  is  determined 
from  Eq.  (7)  to  be  12.27  volts/cm. 

On  the  other  hand,  the  krypton  discharge  is  23%  optically 
thick  and  the  calculated  voltage  gradient  is  8%  too  high.  Since  the 
xenon  discharge  is  80%  optically  thick  its  voltage  gradient  should  be 
(80/23)  8%  too  high.  In  other  words  the  calculated  voltage  gradient 
should  be  12.27  volts/cm. 

The  above  two  calculations  give  the  same  result  and  show 
that  the  difference  between  the  calculated  and  measured  voltage  gradients 
increases  in  direct  proportion  to  the  optical  thickness  of  the  discharge. 
However,  the  calculated  core  temperature  is  independent  of  the  optical 
thickness.  These  results  also  show  that  the  -4%  agreement  between  the 
calculated  and  measured  voltage  gradients  for  xenon  is  fortuitous  and 
should  be  27.8%. 

The  purpose  of  the  previous  discussion  is  to  understand  the 
present  limitations  of  the  calculations.  Thus  the  errors  have  been 
discussed  and  estimated.  It  is  consequently  very  easy  to  be  too 
concerned  about  the  limitations  and  overlook  the  successes.  We  believe 
that  the  theoretical  results  are  so  good  for  the  purpose  of  evaluating 
relative  pump  efficiencies  that  the  original  calculations  of  Section  2.1 
are  used  without  modification.  In  fact,  despite  the  errors,  we  believe 
that  the  best  spectral  overlap  calculations  will  be  done  with  calculated 
rather  than  experimental  spectra  because  the  calculated  spectra  are 
better  resolved.  This  increased  resolution  is  particularly  important 
for  line  emitting  arc  discharges  pumping  a  line  absorbing  laser  rod 
such  as  the  rare  earth  activated  lasers. 
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3.  LASER  ABSORPTION  SPECTRA 


The  necessary  input  information  to  obtain  the  pumping 
efficiency  of  a  lamp  is  its  spectral  emittance  and  the  excitation 
spectrum  of  the  lasing  medium.  The  lasing  media  we  wish  to  consider 
are  Nd^^Al^O^  (Nd;YAG)  and  CriAljO^  (Ruby).  The  absorption  spectra 
of  these  crystals  which  closely  follow  the  excitation  spectra  in  the 
region  of  interest  can  be  obtained  by  experimental  measurement.  This 
data  must  then  be  converted  from  a  chart  recorder  to  computer  data 
cards.  This  operation  is  easily  done  for  ruby  since  the  absorption 
spectrum  varies  slowly  with  wavelength  and  will  be  included  in  a 
later  chapter. 

The  Nd:YAG  absorption  spectrum  consists  of  a  large  number  of 
narrow  lines  and  requires  considerably  more  data  points  to  describe  it 
compared  to  ruby.  However,  prior  to  this  contract,  we  have  developed  a 
program  which  requires  knowing  only  a  few  data  points  per  line  to 
calculate  and  describe  the  spectrum.  The  general  interpolation 
technique  used  is  called  GOOP,  and  for  determining  the  YAG  spectrum 
it  is  called  YAGOOP.  The  YAGOOP  program  is  complete  and  available  for 
use.^  The  general  technique  of  YAGOOP  is  to  fit  a  series  of  Lorenzian 
profiles  to  the  measured  spectrum.  Each  Lorenzian  profile  is  completely 
specified  by  its  center  frequency,  peak  value  and  width  at  half  the 
peak  value.  Prior  to  this  contract  the  Nd.’YAG  absorption  spectrum 
shown  in  Fig.  24  was  fitted  by  YAGOOP.  Fifty-one  Lorenzian  profiles 
are  needed  to  reproduce  the  measured  spectrum.  In  other  words,  153 
specified  variables  give  a  complete  mathematical  description  of  the 
absorption  spectrum.  The  details  of  this  fit  are  given  in  Ref.  1  and 
will  not  be  described  here. 
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4.  SYSTEM  EVALUATION 


How  well  does  a  particular  lamp  pump  a  given  laser  rod?  This 
is  a  question  of  considerable  interest  which  we  will  attempt  to 
answer  in  this  chapter.  In  particular  the  efficiencies  of  the  potassium- 
mercury,  tungsten,  krypton,  and  xenon  lamps  specified  in  Chapter  2  will 
be  evaluated  for  pumping  slab  lasers  of  different  thicknesses. 

The  system  which  is  considered  is  a  simplication  of  actual 
laser  systems.  The  lamp  is  considered  to  be  an  infinitely  long  cylinder 
whose  light  is  focused  by  an  idealized  optical  system  onto  an  infinite 
laser  slab  of  thickness  t.  The  lamp  light  is  transmitted  to  the  slab 
through  a  perfectly  transparent  medium  and  perpendicularly  strikes  the 
slab  with  the  same  emittance  at  which  it  left  the  lamp.  The  light  which 
enters  the  slab  is  then  absorbed  or  transmitted.  The  absorbed  light  pumps 
the  laser  ion  whereas  the  transmitted  light  is  lost  to  the  system.  Although 
this  simple  system  is  obviously  different  than  real  laser  system  it  will 
be  shown  that  it  produces  excellent  quantitative  results. 

4. 1  Lamp  Efficiency 

The  lamp  efficiency  is  defined  here  as  the  percentage  of 

O 

electrical  power  into  the  lamp  which  is  emitted  in  the  3,000  to  10,000  A 

wavelength  region  which  includes  all  the  useful  absorption  bands  of 

Nd : YAG.  The  input  power  per  unit  surface  area  for  the  four  lamps  under 

2 

consideration  was  given  in  Chapter  2.  The  output  power  per  cm  is 
obtained  by  numerically  integrating  the  spectral  emittance  over  the 

O 

5,000  to  10,000  A  region.  This  integration  has  an  accuracy  of  better 
than  5%  and  shows  that  the  lamp  efficiencies  are  67%,  36%,  35%,  and 
14%  for  the  potassium-mercury,  tungsten,  xenon,  and  krypton  lamps, 
respectively.  Note  that  the  xenon  lamp  is  2.5  times  more  efficient 
than  the  krypton  lamp.  These  lamp  efficiencies  do  not  include  electrode 
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losses  which  are  normally  about  15%  of  the  input  power  to  a  3"  lamp 
and  can  be  included  in  the  above  efficiencies  by  a  multiplication 
factor  of  0.87. 


The  efficiency  of  the  simple  system  described  in  the  beginning 
of  this  chapter  is  particularly  interesting  if  the  slab  is  0.7%  Nd:YAG. 
This  laser  material  is  a  line  absorber  and  is  widely  used  in  laser 
systems.  The  former  property  makes  it  a  difficult  test  of  the  lamp 
emission  calculations. 

4.2  Slab  Efficiency 

The  slab  efficiency  is  defined  here  as  the  percentage  of  the 
2 

power  density  (watts/cm  )  incident  on  the  slab  which  is  emitted  at  the 
fluorescing  wavelength  of  1.064  y.  In  other  words,  the  slab  efficiency 
measures  the  effectiveness  of  the  slab  in  absorbing  incident  energy 
and  transforming  this  absorbed  energy  to  useful  emission.  The  power 
density  incident  on  the  slab  was  calculated  in  the  last  subsection. 

The  power  density  emitted  at  1.064  y  is  determined  by  the  spectral 
overlap: 
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Ijl(v)[l  -  exp (-a(v)t) ] (vq/ v)Q(v)dv 


(10) 


where  1^  is  the  spectral  emittance  of  the  lamp,  a  is  the  spectral 
absorptivity  of  the  slab  and  Q  is  the  quantum  efficiency  of  the  laser 
material.  The  vq  is  the  fluorescent  frequency  corresponding  to  1.064  y; 

and  are  the  frequencies  corresponding  to  wavelength  of  1.0  and 
0.5  y,  the  wavelength  region  of  the  lamp  emission.  The  quantum 
efficiency  is  assumed  to  be  unity.  The  spectral  emittance  was  calculated 
in  Chapter  2,  and  the  slab  absorptivity  was  obtained  as  described  in 
Chapter  3.  Thus  everything  in  Eq .  (10)  is  specified  and  the  integrand 
can  be  numerically  integrated.  Before  doing  the  integration,  however, 
it  is  interesting  to  study  the  integrand,  I,  for  various  values  of  slab 
thickness. 
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The  spectral  overlap  integrand  for  the  tungsten  lamp  is  shown 
in  Fig.  25.  As  the  thickness  is  increased  the  broad  emission  bands 
from  the  tungsten  are  increasing  absorbed  by  the  unsaturated  absorption 
lines.  the  absorption  lines  become  saturated  ,the  spectral  distribution 

of  the  tungsten  lamp  becomes  prominent  and  is  very  evident  in  Fig.  25f. 

It  would  be  too  bulky  to  include  all  the  spectral  overlap 
integrand  plots.  Therefore,  only  the  results  for  a  1.2  cm  slab  thickness 
are  shown  in  Figs.  26,  27  and  28  for  the  remaining  tiiree  lamps.  These 
plots  clearly  show  that  most  of  the  spectral  overlap  contribution  i« 

O 

from  above  7,000  A.  The  line  emission  of  these  three  lamps  gives  a 
more  jagged  appearance  to  the  plots  than  present  in  the  tungsten  plots. 
Xenon  obviously  has  the  worst  spectral  overlap. 

The  spectral  overlap  integrand  is  numercially  integrated  to 
an  accuracy  of  at  least  5%,  and  the  slab  efficiency  is  obtained  by 
dividing  by  the  incident  power  density.  The  results  are  shown  in 
Fig.  29.  The  krypton  lamp  clearly  has  the  best  spectral  match  with  a 
Nd :YAG  slab. 

A.  3  Pump  Efficiency 

The  pump  efficiency  is  defined  here  as  the  percentages  of  the 
electrical  power  density  into  the  lamp  (neglecting  electrode  losses) 
which  is  emitted  by  the  slab  at  the  fluorescing  wavelength  of  1.064  u. 

Thus  the  pump  efficiency  measures  the  effectiveness  of  the  slab  laser 
system.  Since  the  slab  properties  remain  constant  as  tie  optical  pump 
is  changed,  it  is  really  a  measure  of  the  optical  pump  effectiveness 
and  explains  the  tern  "pump"  efficiency. 

The  pump  efficiency  is  the  product  of  the  lamp  efficiency 
and  the  slab  efficiency  and  is  the  final  goal  of  these  calculations 
since  it  is  the  best  method  of  evaluating  slab  laser  systems.  All  the 
required  numbers  have  been  calculated  in  the  previous  subsections  and 
the  pump  efficiencies  are  plotted  in  Fig.  30.  The  most  efficient  pump 
by  far  is  the  potassium-mercury  lamp.  This  is  a  result  of  being  the 
most  efficient  lamp  of  the  four  and  hating  a  very  good  spectral  overlap, 
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although  the  spectral  overlap  of  the  tungsten  lamp  is  considerably  better 
than  that  of  the  xenon  lamp.  The  krypton  lamp  is  quite  inefficient  but 
the  spectral  overlap  is  so  large  that  the  pump  efficiency  is  very  close 
to  that  of  the  tungsten  lamp. 

A  comparison  of  the  krypton  and  tungsten  pump  efficiencies 
is  rather  interesting.  The  tungsten  pump  is  more  efficient  for  slabs 
thinner  than  1.2  cm  or  thicker  than  2.4  cm.  The  krypton  pump  is  very 
slightly  more  efficient  in  the  intermediate  range  of  thicknesses. 

This  behavior  demonstrates  how  the  absorption  saturation  of  individual 
lines  can  change  the  relative  efficiencies  of  two  pumps  with  pump 
efficiencies  which  are  very  close  to  each  other. 

These  results  can  be  compared  with  experimental  laser  measure¬ 
ments.  Appendix  1  shows  our  published  results  on  comparing  tungsten, 
krypton,  and  xenon  lamps  for  pumping  a  5  mm  diameter,  50  mm  long  Nd:YAG 
laser  rod.  Figure  31  includes  these  results  and  also  shows  some  more 
recent  work  on  potassium-mercury  lamps.  These  lamps  appear  to  be  very 
promising  as  Nd:YAG  pumps.  Currently  they  must  operate  in  a  vacuum 
to  prevent  atmospheric  attack  of  the  Nb  end  caps  which  operate  at 
temperatures  in  excess  of  700°C.  The  quartz  vacuum  jacket  presently 
used  degrades  the  coupling  efficiency  to  the  laser  rod  and  the  vacuum 
operation  considerably  decreases  the  input  power  to  the  lamp  since  the 
sapphire  walls  of  the  lamp  are  cooled  only  by  radiation.  By  eliminating 
the  vacuum  jacket  and  protecting  the  end  cap*.  through  coatings  or  other 
means  we  have  predicted  previous  to  this  contract  the  experimental 
results  designated  as  air-cooled  K-Hg  lamp. 

Although  the  slope  efficiencies  are  easily  found  from  Fig.  31, 
the  relationship  between  the  rod  diameter  and  the  slab  thickness  must 
still  be  determined  in  order  to  compare  the  measured  and  calculated 
efficiencies.  An  approximation  can  be  made  by  equating  the  slab 
thickness  to  the  effective  chord  length  in  the  rod,  c: 

t  ■  cD  (11a) 
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where  D  is  the  ratio  of  the  activator  doping  in  the  rod  to  the  doping 
in  the  slab.  The  effective  chord  length  in  a  cylindrical  rod  of  diameter 
d  is  shown  in  Appendix  2  to  be  given  by 

i.  -d-yfa  +  o-27)d  <Ub) 

r 

This  calculation  does  not  include  the  cavity  coupling  efficiency  and 
assumes  a  Lambertian  distribution  of  the  light  intensity  in  the  rod.  For 
the  1.3%  doped,  0.5  x  5  cm  rod  used  for  the  experimental  measurements 
of  Fig.  31,  Eqs.  (11a)  and  (lib)  give  an  equivalent  slab  thickness  of 
0.84  cm  (see  Appendix  1).  The  experimental  slope  efficiencies  are 
included  in  Fig.  30  at  this  slab  thickness.  The  agreement  between  the 
calculated  and  measured  relative  pump  efficiencies  is  excellent.  This 
agreement  is  rather  surprising  when  all  the  approximations  are  considered. 

The  comparison  for  the  tungsten  pump  only  tests  the  spectral 
overlap  calculation  and  the  rod  to  slab  thickness  transformation  since 
a  measured  lamp  spectrum  was  used.  The  agreement  at  the  sharply  rising 
curve  for  the  potassium-mercury  pump  is  good  considering  the  crude 
"measured"  value.  This  comparison  is  another  test  for  the  spectral 
overlap  and  rod  to  slab  transformation  calculations. 

The  calculated  xenon  pump  efficiency  is  too  low  and  must  be 
attributed  to  the  low  value  calculated  for  the  contiuum  radiation.  The 

calculated  value  can  be  adjusted  by  reducing  the  total  lamp  radiation 

-3  -1 

from  165  to  155  watts  cm  ster  and  by  assuming  that  64%  of  the  radiation 
is  continuum  (see  Section  2.3).  The  lamp  efficiency  is  thus  reduced  from 
35.2%  to  33.1%.  Since  the  radiation  contributing  most  to  the  spectral 
overlap  is  the  continuum  radiation,  the  slab  efficiency  should  be  close 
to  that  for  tungsten.  The  adjusted  xenon  pump  efficiency  is  2.5%  which 
is  in  good  agreement  with  the  measured  laser  efficiency  of  3.1%. 
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The  real  test  for  the  theoretical  calculations  is  given  by 
a  comparison  of  the  calculated  and  measured  results  for  the  krypton 
pump.  The  continuum  calculation  is  good  and  no  adjustments  need  to 
made.  The  measured  value  of  3.7%  is  30%  higher  than  the  calculated 
value  of  2.6%.  This  agreement  is  excellent  in  view  of  the  lengthy 
ab  initio  calculation.  This  is  the  first  time  such  a  calculation  of 
this  type  has  succeeded.' 

4.  4  Activator  Efficiency 

Although  the  electrical  input  energy  and  the  pump  efficiency 

is  all  that  is  normally  needed,  another  efficiency  is  interesting  to 

know.  This  additional  efficiency  will  be  called  here  the  activator 

efficiency  and  is  the  percentage  of  time  spent  fluorescing  by  an 

-2 

average  activator  ion.  If  the  slab  emits  R  watts  cm  ,  each  activator 

3 

emits  R/tn  watts  where  n  is  the  number  of  activator  ions  per  cm  .  Each 

-4 

emitted  photon  has  an  energy  of  hc/AQ  where  is  1.064  x  10  cm  for 
Nd:YAG.  Thus  the  number  of  photons  emitted  per  second  by  the  activator 
ion  is  given  by 


N  -  (R/tn) (A  /he)  (12) 

o 

The  fraction  of  time  spent  by  the  activator  ion  in  the  fluorescing 
process  is  simply  N  times  the  lifetime  of  the  excited  activator  ion, 
x.  The  activator  efficiency  is  thus 


E  -  Nt  (13) 

a 
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The  lifetime  has  been  measured  at  230  ysec.  Equation  (13)  can  then 
be  evaluated  and  the  results  are  plotted  in  Fig.  32  (Note  that  the 
input  power  densities  have  been  normalized).  The  low  values  show  that 
the  use  of  the  Nd  ion  is  surprisingly  inefficient  because  an  activator 
never  spends  more  than  about  0.06%  of  its  time  in  fluorescence. 


32 


4. 5  Discussion 

The  important  results  in  this  chapter  are  summarized  in  Table 
VII.  The  agreement  between  the  relative  experimental  laser  efficiencies 
and  the  relative  calculated  pump  efficiencies  is  excellent!  Ranked  in 
order  of  decreasing  lamp  efficiency,  the  four  lamps  are  potassium-mercury, 
tungsten,  xenon  and  krypton.  The  ranking  for  the  greatest  slab  efficiency 
is  krypton,  potassium-mercury,  tungsten,  and  xenon  whereas  the  ranking 
for  the  largest  pump  efficiency  is  potassium-mercury,  tungsten,  krypton, 
and  xenon.  These  rankings  clearly  indicate  the  interplay  of  the  lamp 
and  slab  efficiencies  to  give  the  pump  efficiency. 

The  pump  efficiencies  as  a  function  of  thickness  are  important 
but  it  would  be  more  convenient  if  pumps  could  be  ranked  by  a  single 
number.  One  possible  system  would  be  to  use  the  activator  efficiency 
at  zero  thickness.  This  efficiency  is  finite  and  is  independent  of  the 
slab  thickness.  It  is  given  in  Table  VII  in  the  Eq  column.  This  ranking 
system  which  is  independent  of  slab  thickness  would  rank  the  pumps  in 
descending  order  as  the  potassium-mercury,  tungsten,  krypton,  and  xenon 
pumps . 
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5.  SUMMARY 

The  emission  spectra  of  a  xenon  and  a  krypton  arc  discharge 
have  been  measured  on  this  contract,  and  the  absorption  spectra  of 
Nd:YAG  and  ruby  were  measured  previous  to  this  contract.  The  xenon  and 
krypton  arc  discharge  spectra  have  been  calculated  from  theory.  The 
line  spectra  (which  comprise  about  50%  of  the  total  emission)  agree  so 
well  with  the  measured  spectra  that  one  cannot  determine  which  is  the 
more  accurate,  although  the  calculated  line  profiles  appear  to  be  more 
accurate  than  the  measured  line  profiles  (due  to  the  low  resolving  power 
of  the  spectrometer).  The  calculated  continua  only  agree  with  the 
measured  continua  to  within  a  factor  of  two.  In  order  to  save  time  and 
money,  available  measured  emission  spectra  of  a  tungsten  and  a  potassium- 
mercury  lamp  have  been  used. 

The  overlap  between  the  emission  spectra  of  these  four  lamps 
and  the  Nd:YAG  absorption  spectrum  has  been  calculated  as  a  function  of 
the  slab  laser  thickness.  Combined  with  the  calculated  lamp  efficiencies, 
the  spectral  overlap  yields  the  system  efficiency  for  an  optically  pumped 
slab  laser.  A  cylindrical  rod  laser  in  a  cavity  was  approximately  trans¬ 
formed  to  an  equivalent  slab  thickness  of  a  fixed  thickness  and  thus 
actual  laser  system  efficiencies  were  calculated. 

The  agreement  between  the  theoretical  and  measured  relative 
laser  efficiencies  for  the  four  lamps  is  excellent.  This  agreement  is 
remarkable  because  the  calculations  for  the  krypton  and  xenon  pumps 
only  used:  (1)  lamp  dimensions,  (2)  lamp  gas,  (3)  lamp  pressure  and 
current,  (A)  laser  rod  absorption  spectra,  and  (5)  laser  rod 
dimensions  and  activator  doping  as  the  necessary  laser  system  variables 
which  have  to  be  experimentally  measured.  Therefore,  for  the 
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first  time  it  has  been  shown  that  relative  laser  system  efficiencies  can  be 
calculated  from  very  basic  laser  system  parameters.  The  results  show 
that  the  potassium-mercury  lamp  is  four  times  more  efficient  in  pumping 
Nd:YAG  for  all  rod  diameters  of  interest  than  the  tungsten  and  krypton 
lamps  which  are  slightly  better  than  the  xenon  lamp.  The  calculations 
clearly  demonstrate  how  the  theory  can  be  applied  to  evaluate  simple 
optical  pumps  for  a  given  laser  rod,  thereby  bypassing  costly  and 
lengthy  experimental  measurements. 
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6.  INTRODUCTION 

29  30 

The  ruby  laser  was  the  first  successful  laser  ’  and  Is  still 
not  completely  understood.  Since  it  is  a  three-level  laser^  more  than 
50%  of  the  ground  state  chromium  ions  must  be  excited  before  lasing  can 
occur.  This  requirement  necessitates  a  very  powerful  flash  lamp  which 
can  invert  the  population  before  the  excited  ions  can  return  to  the 
ground  state  via  spontaneous  emission  (fluorescence). 

The  required  50%  inversion  automatically  decreases  the  laser 
system  efficiency.  However,  theoretical  calculations  of  this  and  other 
losses  in  ruby  amplifiers  predict  laser  efficiencies  which  are  usually 
an  order  of  magnitude  higher  than  the  measured  values.  This  discrepancy 
is  still  unexplained  although  many  possible  explanations  have  been  pro¬ 
posed. 

One  proposed  explanation  is  that  the  flash  lamp  is  not  efficiently 
pumping  the  ruby.  This  special  report  investigates  this  problem  by  cal¬ 
culating  the  pumping  efficiency  for  fluorescence.  It  is  found  that  the 
pumping  efficiency  for  fluorescence  is  quite  high  at  the  peak  of  the 
flash  lamp  pulse.  In  fact,  it  is  found  that  a  1.27  cm  diameter  ruby 

amplifier  should  have  an  overall  efficiency  of  6%  based  on  pumping  con- 

2 

siderations  alone.  Depopulation  of  the  E  energy  level  by  excited  state 
absorption  would  decrease  the  efficiency  by  a  small  amount  but  excessive 
bleaching  of  the  ruby  rod  or  premature  arrival  of  the  ruby  oscillator 
pulse  appear  to  be  the  best  explanations  for  the  measured  low  amplifier 
efficiencies. 
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7.  LAMP  EMISSION  SPECTRUM 


The  necessary  input  information  to  obtain  the  pumping 
efficiency  of  a  lamp  for  a  given  laser  material  is  the  spectral  emittance 
of  the  lamp  and  the  excitation  spectrum  of  the  lasing  medium.  The 
excitation  spectrum  of  ruby  will  be  considered  in  the  next  chapter,  where¬ 
as  the  spectral  emittance  of  a  xenon  flash  lamp  is  considered  in  this 
chapter. 

The  spectral  radiance  of  a  typical  xenon  flash  lamp  has  been 
3 1 

measured  by  Church  et  al  and  is  shown  in  Figure  33.  The  12300  °K  black 
body  radiance  is  also  shown  for  comparison.  The  xenon  radiation  is 
clearly  continuum  radiation.  A  comparison  of  the  calculated  radiance 
and  the  measured  radiance  again  illustrates  the  inadequacy  of  the 
continuum  radiation  theory.  For  the  following  calculations  the  measured 
radiance  is  always  used. 

Although  Figure  33  shows  the  peak  radiance  from  the  flash  lamp, 
the  light  pulse  is  actually  very  time  dependent.  This  is  clearly  seen 
in  Figure  34.  The  calculation  in  this  report  neglects  this  time  dependence 
and  simply  uses  the  peak  radiance  values.  This  is  a  good  assumption 
since  by  employing  pulse  forming  techniques  nearly  rectangular  pulses 
may  be  obtained. 

Efficiency  calculations  require  the  electrical  power  into  the 

lamp.  Tills  xenon  lamp,  at  its  peak  radiance,  operated  at  a  current 

2 

density  of  4480  amps/cm  with  an  electrical  conductivity  of  63  mho/cm. 

4  2 

The  electrical  input  power  is  thus  1.0  x  10  watts  per  cm  of  surface 
area.  This  calculation  assumes  that  the  voltage  drop  at  the  electrodes 
is  negligible. 

Although  the  radiance  is  measured,  the  emittance  is  required 
for  the  efficiency  calculations.  The  emittance  is  obtained  to  a  good 
approximation  by  multiplying  the  radiance  by  tt  steradians,  the  value  for 
a  perfect  Lambertian  radiator. 


8.  RUBY  ABSORPTION  SPECTRA 


The  absorption  spectra  of  ruby  closely  follow  the  excitation 

spectra  in  the  region  of  interest  and  are  easily  measured  with  a  Cary 

32 

spectrophotometer.  The  results  of  Maiman  are  shown  in  Figure  35.  Since 

ruby  is  anisotropic,  two  absorption  spectra  are  shown,  one  for  light 

propagating  parallel  and  one  for  light  propagating  perpendicular  to  the 

c-axis  of  the  ruby  crystal.  The  chromium  concentration  was  determined^ 

by  chemical  analysis  to  be  0.0515  weight  percent  of  Cr,0  :A190„  which 

19  ^  ^  ^ 

corresponds  to  a  chromium  ion  density  of  1.62  x  10  ions/cm3.  The 

absorption  data  for  the  R  lines  are  shown  in  Figure  35  but  were  not  used 

33,34 

in  the  present  calculations.  Excited  state  absorption  was  also  not 
considered;  however,  this  is  expected  to  cause  no  major  error. 

A  comparison  of  the  two  spectra  in  Figure  35  indicates  that 
light  propagating  parallel  to  the  c-axis  should  pump  the  ruby  crystal 
more  effectively  than  light  propagating  perpendicular  to  the  c-axis. 

We  will  show  later  that  this  is  indeed  the  case,  although  the  difference 
is  not  as  large  as  one  might  expect  from  the  absorption  spectra. 

The  absorption  coefficient  given  in  Figure  35  is  for  small 
amplitude  pumping  when  virtually  all  the  absorbing  ions  are  in  the  ground 
state.  However,  under  hard  pumping  at  least  half  of  the  ions  are  in  the 
excited  state  and  unavailable  for  excitation.  In  this  case,  the  absorption 
coefficients  are  half  the  values  shown.  During  long-pulse  laser  operation, 
stimulated  emission  automatically  keeps  the  ground  state  about  one-half 
filled  and  therefore  the  efficiencies  are  calculated  from  absorption 
coefficients  which  are  one-half  the  values  given  in  Figure  35. 
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9.  SYSTEM  EVALUATION 


How  well  does  the  xenon  lamp  pump  the  ruby?  This  Is  the  basic 
question  which  will  be  examined  in  this  chapter.  However,  only  the 
fluorescence  pumping  efficiency  is  calculated;  the  details  of  the  lasing 
process  which  could  modify  this  efficiency  are  not  considered. 

The  model  used  is  a  simplification  of  the  actual  laser  system. 

It  is  thought  to  be  a  good  approximation  to  the  "first  pass"  in  a 
focusing  elliptical  or  spherical  pumping  cavity.  The  lamp  is  considered 
to  be  an  infinitely  long  cylinder  whose  radiation  impinges  onto  an 
infinitely  long  ruby  slab  of  thickness  t.  The  lamp  radiation  is  trans¬ 
mitted  to  the  slab  through  a  perfectly  transparent  medium  and  perpendicularly 
strikes  the  ruby  slab  with  the  same  emittance  at  which  it  left  the  lamp. 

The  light  which  enters  the  slab  is  then  absorbed  or  transmitted.  The 
absorbed  light  pumps  the  chromium  ions  whereas  the  transmitted  light  is 
lost  to  the  system.  Although  this  simple  model  is  obviously  different 
than  the  real  laser  system,  It  produced  excellent  results  for  the 
pumping  of  Nd:YAG  and  should  be  equally  applicable  to  the  pumping  of 
ruby. 

9 . 1  Lamp  Efficiency 

The  lamp  efficiency  is  defined  here  as  the  percentage  of 

0 

electrical  power  into  the  lamp  which  is  emitted  in  the  3000  to  7000  A 
wavelength  region  which  includes  all  the  useful  absorption  bands  of  ruby. 

The  input  power  per  unit  surface  area  for  the  xenon  lamp  was  calculated 
in  Chapter  7,  The  output  power  per  unit  surface  area  is  obtained  by 

O 

numerically  integrating  the  spectral  emittance  over  the  3000  to  7000  A 
region.  This  integration  has  an  accuracy  of  better  than  5%  and  shows 
that  the  xenon  lamp  efficiency  is  36%.  This  calculation  assumes  that 
the  electrode  losses  are  negligible. 
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9.  2  Slab  Efficiency 

The  slab  efficiency  is  defined  here  as  the  percentage  of  the 
2 

power  density  (watts/cm  )  incident  on  the  slab  which  is  emitted  as  ruby 
fluorescence.  In  other  words,  the  slab  efficiency  measures  the  effective¬ 
ness  of  the  slab  in  absorbing  incident  energy  and  transforming  this 
absorbed  energy  into  useful  emission.  The  power  density  incident  on 
the  slab  was  calculated  in  the  last  section.  The  power  density  emitted 

O 

at  6943  A  (assuming  that  all  the  ruby  fluorescence  takes  place  at  this 
wavelength)  is  determined  by  the  spectral  overlap: 

v2 

S  -  /  I^v)  l1  "  exP  (_a(v)t)](vo/v)  Q(v)  dv  (10) 

V1 

where  1^  is  the  spectral  emittance  of  the  lamp,  a  is  the  spectral 
absorptivity  of  the  ruby  slab  and  Q  is  the  quantum  efficiency  of  the 
ruby.  As  explained  in  Chapter  8 ,  the  absorptivity  used  in  Equation  (10) 
is  one-half  of  the  values  shown  in  Figure  35.  The  v  is  the  fluorescent 

O  ° 

frequency  corresponding  to  6943  A;  v.  and  v„  are  the  frequencies 

1  4  o 

corresponding  to  wavelengths  of  7000  and  3000  A,  the  effective  wavelength 

region  of  lamp  emission.  The  quantum  efficiency  has  been  measured  by 
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Maiman  and  is  shown  in  Figure  36  (although  the  results  of  Maiman  ’  are 

rather  old,  they  are  in  a  more  convenient  form  to  use  and  do  not 

33  34 

significantly  differ  from  more  recent  results  ’  ).  The  power  efficiency 

curve  shown  in  this  figure  is  the  vq/v  ratio  times  the  quantum  efficiency. 
Thus  everything  in  the  integrand  is  specified. 

The  spectral  overlap  integrand  is  numerically  integrated  to 
an  accuracy  of  at  least  5%,  and  the  slab  efficiency  is  obtained  by  dividing 
by  the  incident  power  density.  The  results  are  shown  in  Figure  37.  As 
expected  from  the  absorption  data,  light  propagating  parallel  to  the  ruby 
c-axis  gives  better  slab  efficiencies. 
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9 . 3  Pump  Efficiency 


The  pump  efficiency  is  defined  here  as  the  percentage  of  the 

electrical  power  density  into  the  lamp  (neglecting  electrode  losses) 

0 

which  is  emitted  by  the  slab  at  the  fluorescing  wavelength  of  6943  A 
(i.e.,  the  product  of  lamp  and  slab  efficiency).  Thus  the  pump  efficiency 
measures  the  fluorescence  effectiveness  of  the  slab  laser  system.  Since 
the  slab  properties  remain  constant  as  the  optical  pump  is  changed,  it 
is  really  a  measure  of  the  optical  pump  effectiveness  and  explains 
the  term  "pump"  efficiency. 

All  the  required  numbers  have  been  calculated  in  the  previous 

sections  and  the  pump  efficiencies  are  plotted  in  Figure  38.  It  was 

shown  in  Section  4.3  that,  to  a  good  approximation,  the  diameter  of  a 

laser  rod  may  be  equated  to  about  the  same  value  of  slab  thickness. 

Consequently,  the  abscissas  of  Figures  5  and  5  may  be  considered  as 

laser  rod  diameters.  For  a  rod  with  the  same  diameter  as  the  xenon  lamp 

(1.27  cm)  the  pumping  efficiency  is  about  6%.  It  is  then  apparent  that 

the  pumping  efficiency  of  ruby  is  quite  high  and  does  not  explain  the 

35 

measured  low  ruby  laser  efficiencies  . 

9 .4  Activator  Efficiency 

Although  the  electrical  input  energy  and  the  pump  efficiency 
are  all  that  are  normally  needed,  another  efficiency  is  needed  for  lasers 
which  are  pumped  hard.  This  additional  efficiency  will  be  called  here 
the  activator  efficiency  and  is  the  percentage  of  time  spent  in  the 
excited  state  by  an  average  chromium  ion.  If  the  ruby  slab  emits  R  watts  per 
unit  surface  area,  each  chromium  ion  emits  R/tn  watts  where  n  is  the 

3 

number  of  activator  ions  per  cm  and  t  is  the  slab  thickness.  Each 

0 

emitted  photon  has  an  energy  of  hc/AQ  where  is  taken  as  6943  A.  Thus 
the  number  of  photons  emitted  per  second  by  each  chromium  ion  is 
approximately  given  by 

N  ■  (R/tn)  (AQ/hc)  (12) 
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Equation  (12)  is  slightly  different  for  a  laser  rod.  If  we 
take  the  rod  diameter  d  to  be  the  same  as  the  slab  thickness  (this  is 

2 

probably  accurate  to  within  10%),  each  chromium  ion  emits  4(irdi)R/ (ird  £)n 
watts,  where  £.  is  the  length  of  the  rod.  In  other  words,  each  chromium 
ion  emits  4R/dn  watts  and  Equation  (12)  becomes 

N  =  (4R/dn)  (A  /he)  (12a) 

o 

for  a  ruby  rod.  This  expression  assumes  that  no  pump  light  is  absorbed 
by  the  ends  of  the  rod. 

The  fraction  of  time  spent  by  the  activator  ion  in  the  excited 
state  is  simply  N  times  the  lifetime  of  the  excited  chromium  ion,  t. 

The  activator  efficiency  is  thus 

E  =  Nt  (13) 

a 

32 

The  measured  lifetime  at  low  pumping  levels  is  3.0  msec  at  300°K.  At 
this  lifetime  the  calculated  activator  efficiency  varies  with  slab  thick¬ 
ness  from  150%  to  1300%!  Since  the  activator  efficiency  can  never  exceed 
100%,  it  is  clear  that  the  activator  efficiency  cannot  be  calculated  as 
described.  The  problem  is  the  excited^s tate  lifetime.  At  low  pumping 
levels  the  measured  values  are  appropriate,  but  at  high  pumping  levels 
many  excited  ions  will  go  to  the  ground  state  via  stimulated  emission. 

If  stimulated  emission  is  prevented,  the  power  absorbed  and  reemitted 
by  the  ruby  decreases  due  to  bleaching.  By  assuming  that  stimulated 
emission  takes  place^the  number  of  stimulated  emission  transitions  is 
sufficient  to  lower  the  average  excited-state  lifetime.  This  lowered 
lifetime  due  to  stimulated  emission  can  be  calculated  from  Equation  (13) 
by  setting  the  activator  efficiency  at  the  desired  value  of  50%.  The 
results  are  plotted  in  Figure  39,  which  shovs  that  the  lifetime  of  a  1.27 
cm  diameter  ruby  rod  is  lowered  from  30  msec  at  low  pumping  levels  to 
about  120p  sec  at  the  high  pumping  levels  of  the  xenon  flash  lamp  used 
for  these  calculations. 


44 


10 .  SUMMARY 


The  calculated  pump  efficiency  of  about  6%  fora  1.27  cm  diameter 

ruby  rod  is  considerably  higher  than  measured  laser  efficiencies,  particularly 

35 

in  oscillator-amplifier  combinations.  One  possible  cause  for  this  difference 

2 

is  depopulation  of  the  E  state  (see  Figure  40)  via  transitions  to  the  charge 
36  37  0 

transfer  band.  1  This  transition  occurs  at  3600  A  and  decreases  the  pump 
efficiency.  However,  the  xenon  flash  lamp  radiance  shown  in  Figure  33 
decreases  in  this  wavelength  region,  and  can  be  filtered  out  without  changing 
the  pump  efficiency  appreciably.  Depopulation  should  consequently  have  a 
measurable  but  small  effect  on  the  pump  efficiency. 

Another  possible  source  of  the  discrepancy  between  calculated 
pump  efficiencies  and  laser  efficiencies  is  the  time  dependence  of  the 
pump  pulse.  Figure  34  shows  the  time  history  of  one  xenon  line  and  Figure 
shows  the  energy  level  diagram  of  ruby  with  measured  lifetimes.  As  the 
pulse  radiance  increases,  the  pump  spectrum  changes  and  becomes  a  more 
efficient  pump.  Thus,  the  initial  part  of  the  pump  pulse  decreases  the 
total  pump  efficiency.  This  effect  is  also  present  during  the  final 
stages  of  the  pulse  when  the  pump  radiance  is  decreasing.  This  loss  can 
be  minimized  by  pulse  forming  networks. 

The  above  arguments  hold  for  oscillators  as  well  as  amplifiers. 
Amplifiers  have  an  additional  mechanism  for  loss  of  efficiency  which  we 
believe  could  be  very  large.  Above  it  was  assumed  that  stimulated  emission 
kept  the  ground  state  population  at  50%.  In  oscillators  this  occurs 
automatically,  but  in  amplifiers  this  condition  depends  on  the  proper 
timing  of  the  input  laser  energy.  If  the  input  laser  power  arrives  too 
soon, it  will  be  absorbed  by  the  amplifier  rod  rather  than  amplified.  If 
it  arrives  too  late^the  pump  pulse  will  have  bleached  the  rod  to  the 
point  where  the  pump  radiation  is  no  longer  being  absorbed  by  the  rod. 

This  problem  can  only  be  accurately  analyzed  by  a  time-dependent  study. 


>* 
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A  time-dependent  analysis  of  lamp  emission,  excited  populations,  absorption 
coefficients,  quantum  efficiencies,  and  pump  efficiencies  would  give  a 
better  understanding  of  ruby  amplifiers  and  perhaps  lead  to  increased 
efficiencies . 
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RADIAL  WAVE  FUNCTIONS 


47 


11.  PARAMETER  FITTING  FROM  SLATER  INTEGRALS 


« 


The  objective  of  this  work  phase  is  to  determine  atomic  radial 
wave  functions  by  fitting  a  power  series  expansion  with  parametric  co¬ 
efficients  of  the  radial  wave  function  to  experimentally  determined  Slater 
integrals  of  the  radial  wave  functions.  These  wave  functions  can  then 
be  used  to  calculate  the  transition  probabilities  required  in  the  arc 
discharge  emission  calculations.  The  parameter  fitting  (PF)  technique 
has  the  advantage  of  being  independent  of  the  complex  electronic  shell 
structure.  This  independence  is  in  contrast  to  the  conventional  methods 
such  as  the  Thomas-Fermi  and  Bates  and  Damgaard  coulomb  approximations 
which  only  work  for  simple  shells  of  one  or  two  electrons.  We  con¬ 
sequently  expect  PF  to  succeed  where  all  other  methods  fail. 

To  initiate  PF  the  constraints  must  be  enumerated  which  a  radial 

wave  function  must  satisfy,  and  reasonable  forms  for  the  PF  power  series 

expansion  with  fitting  parameters  must  be  formulated.  The  constraints 

and  expansions  are  given  in  the  next  two  sections.  The  selection  of  the 

numerical  procedure  to  adjust  the  fitting  parameters  to  satisfy  the 

constraints  and  reproduce  the  wxperimentally  determined  Slater  integrals 

is  a  critical  selection  which  forms  the  heart  of  the  PF  technique.  Two 

iterative  numerical  procedures  are  examined  in  Section  11.3.  The  initiation 

of  both  procedures  requires  an  approximation  of  the  fitting  parameters. 

This  approximation  is  best  achieved  with  the  aid  of  screening  constants 

which  are  discussed  in  Section  11.4.  Finally,  the  PF  technique  must  be 

proven  in  a  case  where  the  radial  wave  function  is  known.  The  Hartree- 
38 

Fock  calculations  of  the  Cu  4s  Slater  Integrals  and  radial  wave  function 
provide  an  excellent  test  case. 
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11.1  PF  CONSTRAINTS 


An  atomic  wave  function  Is  usually  given  by 

|n£m)  =  r_1Pn)i(r)  Y£m<e,*)  (14) 

where  P  (r)  Is  the  normalized  radial  wave  function  to  which  PF  will  be 

nn 

applied  and  is  given  below  in  atomic  units.  The  radial  wave  function 
must  satisfy  the  following  general  requirements. 

38 

As  r  approaches  zero  Mann  has  shown  that 

Pn«<r>  ^  +  1(1  ■  TTT  r>  (15) 

where  Z  is  the  nuclear  charge.  On  the  other  hand,  as  r  approaches 

39 

infinity  Handy  et  al.  have  shown  that  the  asymptotic  form  of  the  radial 
wave  function  is 

P„o<r>rT*  f<*>  exp  <-<-E  )1/2r)  (16) 

n*.  r  -*»  n»a» 

where  f  (r)  is  a  function  of  r  which  does  not  have  as  critical  dependence 
on  r  as  the  exponential.  The  E  ,  ,  (whose  value  is  negative)  is  the 

3g  n  “ 

Hartree-Fock  binding  energy  of  the  least  bound  electron  of  the  atom, 
l.e.,  not  necessarily  of  the  electron  under  consideration.  The 
constraints  in  Eq.  (15)  and  (16)  result  from  the  Hartree-Fock  Equation 
and  may  not  be  physically  correct.  However,  because  the  Hartree-Fock 
formulation  is  so  extensively  used  and  generally  gives  good  results, 
these  constraints  will  be  accepted  as  valid  for  PF. 

Three  other  common  constraints  are  the  number  of  nodes  and 
orthonormality.  The  number  of  nodes,  excluding  end  values,  is  given  by 
n  -  i  -  1.  This  condition  will  be  used  although  certain  wave  functions 
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are  more  accurately  obtained  when  this  condition  Is  dropped.  The 
orthonormality  condition: 


K  Vr)  p»'t(r> 


nn 


is  one  of  mathematical  convenience. 


(17) 


Of  the  above  constraints,  only  the  asymptotic  form  given  in 

Eq.  (16)  poses  any  problem  in  principle.  The  f (r)  is  unknown  as  is  the 
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binding  energy  (except  for  the  ground  states  of  atoms).  Consequently, 
the  implementation  of  this  constraint  requires  an  approximate  approach. 

The  f(r)  will  be  guessed  based  on  our  experience.  The  choice 
of  f(r)  should  not  be  critical  because  (a)  it  is  only  pertinent  in  the 
extreme  outer  portion  of  the  wave  function,  (b)  its  r  dependence  is 
less  than  that  of  the  exponential,  and  (c)  any  error  in  f(r)  will  be 
largely  ameliorated  by  the  fitting  of  the  adjustable  parameters  to 
reproduce  the  measured  Slater  integrals. 

The  determination  of  the  binding  energy  is  more  important 
since  it  appears  in  the  exponential.  Fortunately,  the  binding  energy 
is  valuable  for  other  purposes  and  approximations  for  its  value  have 
already  been  devised.  The  finding  energy  is  given  by 

E„J  -  -(z  -  <’„e,/n  <18> 

where  n  is  the  principal  quantum  number  and  o  is  a  screening  constant. 
Thus  Z  -  a  ^  is  the  effective  nuclear  charge  as  seen  by  an  electron  in 

the  (n£)th  orbit.  Methods  for  estimating  the  screening  constant  have 

40  41  42 

been  developed  '  and  tables  of  their  values  have  appeared. 
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11.2  PF  BASIS  SETS 


The  choice  of  a  basis  set  of  functions  into  which  the  radial 
wave  function  can  be  expanded  is  a  rather  arbitrary  choice.  Except  for 
trivial  problems  the  true  analytic  form  (if  one  exists)  of  the  radial 
wave  function  is  unknown.  Approximate  analytic  forms  are  usually 
obtained  by  assuming  an  analytic  form  with  variational  parameters  and 
then  varying  these  parameters  until  the  energy  is  minimized.  These 
variational  techniques  are  very  dependent  upon  a  good  choice  for  the 
analytic  form.  On  the  other  hand,  the  least-squares  fit  of  an  analytic 
function  with  fitting  parameters  to  a  numerical  function  is  not  very 
dependent  upon  the  choice  of  the  analytic  form  (this  fact  has  often 
allowed  experimentalists  to  fit  their  data  to  wrong  theories,  and 
vice  versa).  This  advantage  of  PF  will  tend  to  be  offset  by  the  large 
smoothing  of  the  integrand  e.g.,  the  radial  wave  functions,  by  the 
double  Integrals  in  the  Slater  integrals.  In  other  words,  integrals 
tend  to  be  independent  of  the  detailed  integrand  structure  and  this 
independence  allows  poor  radial  wave  functions  to  give  good  Slater 
integrals.  That  is  why  constraints  on  the  analytic  radial  wave  function 
are  imperative. 

Since  the  choice  of  the  PF  basis  set  is  somewhat  arbitrary, 

it  would  be  best  to  investigate  several  sets  so  as  to  obtain  the  best 
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one.  The  five  best  known  sets  are  hydrogenic  orbitals,  Slater 

44  45  46  46 

orbitals,  Boys  orbitals,  Bessel  orbitals  and  Gaussian  orbitals. 

Economic  realities  dictate  that  only  one  basis  set  be  investigated  (at 

least  until  PF  becomes  a  proven  and  useful  tool).  The  Slater  orbitals 

are  chosen  because  they  satisfy  constraints  (15)  and  (16),  are  very 

flexible,  and  are  the  most  commonly  used  basis  set.  The  PF  expansion  is 

they  written  as 
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n  -  1 


w<r) 


AF^Ux)  +  £ 
"l  i 


R  P(0^h  (19) 

t  +  2BjPnj  (bjx) 


+  C  P^o)  i (cx) 
n,n  -  1 


where,  for  Slater  orbitals 


(2a/(2 1  +  2)1)  (2ar)*  +  ^xp  (-ar)  (20) 


This  PF  expansion  has  A,  B  ^ ,  C  a,  b  ^  ,  and  c  for  the  fitting  parameters. 
As  stated  the  fitting  parameters  always  occur  in  pairs,  even  though 
there  frequently  are  an  odd  number  of  Slater  integrals  which  can  be 
used  for  fitting.  Whenever  this  situation  exists  the  value  of  c  will 
be  approximated  by  Eq.  (18).  Since  only  Slater  orbitals  are  considered 
here,  these  fitting  parameters  will  be  called  the  Slater  parameters. 
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11.3  NUMERICAL  PROCEDURES 


The  numerical  evaluation  of  integrals  such  as  the  Slater  in¬ 
tegrals  and  those  in  Eq .  (17)  is  a  mundane,  but  important  numerical 
problem  since  these  integrals  are  evaluated  thousands  of  times  in  a  long 

PF  iteration.  The  most  commonly  used  technique  is  to  use  a  variable 

47  48 

step  size  Simpson's  rule.  Herman  and  Skillman  and  Mann  tabulate 
their  wave  functions  in  a  step  size  which  starts  from  0.01  for  x  between 
0.00  and  0.10  and  is  doubled  for  succeeding  ranges  until  1.28  is  the 
largest  step  size  used  for  large  values  of  x.  It  is  important  to  note 
that  the  variable  which  is  usually  tabulated,  and  used  as  the  variable 
of  integration  is  not  r,  but  x,  which  is  a  Thomas-Fermi  normalized 
radius.  The  relationship  between  the  two  variables  is  given  by 


(9tt2/  128  Z)1/3  x  -  Kx 


(21) 


where  Z  is  the  nuclear  charge. 

Although  Simpson's  rule  is  the  most  convenient  integration 

method  to  use,  it  is  neither  the  fastest  nor  the  most  accurate.  The 

trapezoidal  rule  is  faster  and  only  slightly  less  accurate.  However, 

the  best  integration  procedure  for  integrals  which  contain  radial  wave 

49 

functions  is  probably  Laguerre  integration  of  the  Gaussian  type: 


i  £<x)  dx  \  f  <\> 


(22) 


th 


The  abscissa  x^  are  fixed  at  the  zeros  of  the  n  order  Laguerre 

polynomial  and  the  weights  w^  are  determined  from  the  values  of  the 
ttl 

(n+1)  order  Laguerre  polynomial  at  the  x,  values.  The  weights  and 

K  50 

abscissa  are  tabulated  in  the  literature  for  many  values  of  n.  It  is 
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easily  shown  that  when  f(x)  contains  radial  wave  functions,  the  Laguerre 
integration  rapidly  converges  to  the  correct  value  as  n  Increases  without 
limit. 

A  comparison  of  these  numerical  integration  methods  is  given 

in  Tables  VIII  and  IX  where  the  integrals  are  between  the  Cu  4s  electron 

and  inner  shell  electrons.  The  true  values  of  the  integrals  have  been 
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evaluated  by  Mann.  Simpson's  rule  with  100  points  is  the  most  accurate 
method,  but  extremely  slow.  The  trapezoidal  rule  is  slightly  faster  and 
slightly  less  accurate.  The  15-polnt  Laguerre  integration  is  about  flftv 
times  faster,  but  a  little  too  Inaccurate  for  most  work.  The  48-Doint 
Laguerre  integration  is  five  times  faster  and  as  accurate  as  Simpson's 
rule  integration. 

It  is  clear  that  Laguerre  integration  is  the  best  method  for 
PF,  but  the  application  of  Laguerre  integration  requires  more  facts.  The 
Hartree-Fock  radial  wave  functions  which  appear  in  the  f(x)  of  Eq.  (22) 
are  evaluated  at  the  x^  values  by  six-point  Lagrange  interpolation.  The 
upper  limit  of  the  integral  cannot  practically  be  infinite,  but  must  be 
some  value  of  x  beyond  which  the  Integrand  must  be  negligible.  For  the 
15-point  integration  the  upper  limit  is  simply  taken  to  be  x^,  i.e., 
48.03.  For  the  48-polnt  Integration  the  upper  limit  is  the  last  tabulated 
value  of  the  radial  wave  function  which  is  larger  than  0.0001.  For  the 
Cu  4s  radial  wave  function  this  upper  limit  is  66.46. 

The  selection  of  the  numerical  integration  technique  now 
allows  the  iterative  PF  method  to  be  selected.  A  conventional  problem 
in  physics  is  to  fit  via  least  squares  parameterized  functions  to  a  given 
curve.  The  PF  problem  is  different  in  that  the  parameterized  functions 
are  first  Integrated  and  the  integrals  are  fitted  to  known  values.  This 
difference  poses  a  problem  as  to  which  iterative  procedure  should  be  used. 
A  Newton's  method  for  systems  of  simultaneous  non-linear  equations  and 
a  Direct  Search  method  have  been  tried  in  which  the  sum  of  the  squares 
of  the  differences  between  the  known  and  calculated  radial  Integrals  is 
minimized. 
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The  Newton's  method  for  solving  a  set  of  simultaneous  non¬ 
linear  equations  has  been  slightly  modified  for  PF.  This  method  should 
quickly  converge  and  has  the  advantage  over  other  methods  in  that  if  It 
does  not  work,  the  method  specifies  the  reason  for  the  impasse  and,  hopefully, 
corrective  action  may  be  taken.  The  details  of  the  procedure  are  given 

in  Appendix  3  where  f  is  the  difference  between  the  n^  calculated  and 

n  * 

experimentally  determined  Slater  integral,  and  the  x^  are  the  n  fitting 
parameters  of  the  basis  set. 

A  constant  danger  in  all  non-linear  iterations  is  that  the 
procedure  will  settle  into  a  rather  high  local  minimum  and  be  unable  to 
adequately  search  parameter  space  for  lower  minima  (in  non-linear  problems 
many  minima  usually  exist  and  the  lowest  minimum  is  usually  not  found). 
Unfortunately,  this  danger  was  actually  reality  for  the  Newton's  method. 

It  was  all  too  clear  from  the  beginning  that  the  Newton's  method  was 
inadequate  for  PF.  The  cause  of  the  problem  was  traced  to  the  almost 
linear  dependency  of  some  of  the  Slater  parameters.  A  detailed  analysis 
of  the  problem  is  given  in  Appendix  4. 

The  Direct  Search  method  should  work  in  all  cases  because  it 
is  essentially  an  organized  trial  and  error  technique,  and,  indeed,  it 
does  work  very  well,  if  a  bit  slow.  The  Direct  Search  program  has  been 
internally  developed  by  Westinghouse  and  basically  works  in  the  following 
manner.  From  the  initial  set  of  Slater  parameters,  the  sum  of  the  square 
of  the  differences  between  the  known  and  calculated  radial  integrals  is 
calculated.  The  base  set  of  Slater  parameters  is  equated  to  the  initial 
set  and  a  set  of  exploratory  moves  is  made,  adjusting  individual  Slater 
parameters  one  at  a  time.  The  sum  of  the  squares  is  computed  after  each 
move.  A  particular  move  is  retained  if  it  produces  a  lower  sum  of 
squares,  but  is  otherwise  rejected.  If  a  set  of  exploratory  moves  lowers 
the  sum  of  squares,  the  adjusted  Slater  parameters  replace  the  base 
Slater  parameters  and  the  procedure  is  repeated.  If  a  set  of  exploratory 
moves  does  not  lower  the  sum  of  squares,  the  base  set  of  Slater  parameters 
is  retained  and  a  new  set  of  exploratory  moves  tried.  Direct  Search  con¬ 
tains  a  set  of  rules  for  determining  the  step  size  of  the  exploratory  moves 


for  individual  Slater  parameters.  These  step  sizes  depend  on  the  success 
or  failure  of  previous  exploratory  moves. 
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11.4  SCREENING  CONSTANTS 


Tiie  usage  of  PF  requires  the  approximation  of  the  Slater 
parameters.  The  PF  expansion  is  given  by  Eq.  (19),  but  Is  written  In 
a  more  compact  form  as 


(x) 
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-  t. 

m*l 


A  P  (0) 
m  m 


(ax) 
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(23) 


where  the  P^°^  are  the  Slater  orbitals.  This  expansion  Is  more  sensitive 
m 

to  the  a  values  than  to  the  A  values  since  the  a  parameters  appear  in 
m  mm 

exponentials.  This  expansion  Is  also  lmpllcltely  dependent  on  the  choice 
of  the  M  Slater  orbitals,  i.e.,  on  the  choice  of  the  nt  for  each  orbital 
as  given  In  Eq.  (20)  (Although  the  Slater  orbitals  are  Independent  of  n, 
the  nil  notation  Is  retained  to  Indicate  which  atomic  crbltal  the  Slater 
orbital  Is  approximating).  Since  each  Slater  orbital  is  nodeless,  the 
selection  of  the  M  orbitals  will  be  made  such  that  the  sum  of  orbitals 
reproduces  the  nodes  of  the  P  . (x)  as  closely  as  possible. 


The  selection  of  the  M  orbitals  is  mainly  based  on  experience 
and  is  explored  in  the  next  section.  The  approximation  of  the  A^  is  not 
critical  and  crude  methods  (such  as  guesses)  suffice.  However,  the 
approximation  of  the  a^  is  an  important  function  and  is  based  on  the  con¬ 
cept  of  screening  constants.  From  Eqs.  (16),  (18),  and  (21)  we  see  that 


a»  ■ K  i (z  - 


(24) 


gives  a  method  for  approximating  these  parameters. 

40 

Slater  has  given  some  rules  for  the  estimation  of  screening 

constants  which  are  convenient  to  use,  but  very  inaccurate.  A  more 

41 

accurate  set  of  rules  is  given  by  Burns  who  fit  a  single  Slater  orbital 
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to  more  than  90  Hartree-Fock  wave funcc ions .  The  results  are  given  In 

Table  X.  If  n'  is  less  than  n-2,  a  ,  is  1  and  If  n'  is  larger  than 

n  * 

n  +  2,  0  ,  ,  is  0.  The  screening  constant  in  Eq.  (24)  is  obtained  bv 
n 

summing  the  screening  constants  for  each  n'i'  electron  (excluding  the 
n£  th  electron)  in  the  atom.  Thus, 


0 


ni, 


(25) 


For  example,  the  screening  constant  for  the  Cu  4s  electron  is  given  by 

o4a  -  10+2(0. 90)+6(0.75)+10(0. 50)  -  21.3  (26) 

In  other  words,  the  4s  electron  sees  an  effective  nuclear  charge  of  7.7 
sines  Z  for  copper  is  29. 


An  alternate  method  of  obtaining  screening  constants  is  to 

38 

use  Eq.  (18)  by  putting  Mann’s  values  of  the  binding  energy  on  the  left 
hand  side  and  solving  for  the  screening  constant.  This  can  be  done  for 
the  first  103  atoms  and  the  resulting  screening  constants  averaged.  The 
result  is  given  in  Table  XI.  In  this  case  the  screening  constant  for  the 
Cu  4s  electron  is  given  by 


o4a  -  10+2(0. 875)+6(0.776)+10(0. 918)  -  25.586  (27) 

or  an  effective  nuclear  charge  of  3.414.  This  value  is  considerably  lower 
than  the  value  calculated  from  Bums'  rules.  This  is  consistent  with  the 
different  rule  derivations  since  Bums  fit  the  entire  wavefunction  whereas 
the  Mann's  rules  are  obtained  from  Eq.  (16)  which  is  only  valid  for  the 
outer  portion  of  the  wavefunction. 
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11.5  APPLICATION  TO  THE  Cu  4s  ELECTRON 


The  final  proof  of  any  new  numerical  technique  is  in  the  ap¬ 
plication  of  the  technique  to  a  good  test  case.  A  good  test  case  for 
PF  is  the  Cu  4s  electronic  radial  wave  function.  This  wave  function  is 
sufficiently  compressed  that  15-point  Laguerre  integration  can  be  used, 
thereby  reducing  computer  expenses.  On  the  other  hand  it  is  sufficiently 
complicated,  e.g. ,  it  has  three  nodes,  to  tax  any  radial  wave  function 
generating  technique.  It  may  at  first  appear  that  a  s  wave  function  is 
too  simple  and  that  a  p  or  d  wave  function  should  be  selected.  However, 
these  latter  wave  functions  are  too  extended  for  15-point  Laguerre  in¬ 
tegration  or  have  too  few  nodes.  Furthermore,  although  there  is  a  con¬ 
siderable  difference  between  different  l  value  angular  wave  functions, 
there  is  no  explicit  difference  between  different  l  value  radial  wave 
functions . 

Values  of  some  <rn>  integrals  and  Slater  integrals  between 
the  4s  electron  and  some  inner  shell  electrons  are  given  in  Table  XII. 

The  F°  Slater  integral,  the  normalization  integral,  and  the  three  ortho¬ 
gonality  integrals  with  the  3s,  2s,  and  Is  electrons  gives  values  of 
five  integrals  which  the  PF  4s  wave  function  must  reproduce.  Thus,  the 

PF  expansion  will  have  five  fitting  parameters.  The  3s,  2s,  and  Is  wave 

48 

functions  needed  for  the  orthogonality  integrals  will  be  those  of  Mann, 
as  will  be  the  4s  wave  function  which  will  be  considered  as  the  "true" 
wave  function  with  which  the  PF  wave  function  will  be  compared. 

There  are  many  different  ways  of  selecting  the  initial  PF 
expansion.  Seven  different  methods  were  tried;  the  first  four  expansions 
have  three  Slater  orbitals  with  the  last  exponent  being  fixed,  and  the 
last  three  expansions  have  five  Slater  orbitals  with  all  the  exponents 
being  fixed.  Thus,  the  first  four  expansions  can  have  at  most  two  nodes 
while  trying  to  reproduce  a  three  node  wave  function.  This  may  lead 
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to  some  inaccuracy  near  the  origin.  However,  these  same  expansions  have 
two  fitting  parameters  in  the  exponents  which  would  appenr  to  give  them 
a  greater  flexibility  than  the  five-term  expansions.  These  seven  ex¬ 
pansions  were  tested  by  comparing,  after  100  seconds  of  computing  time 
on  a  Burroughs  B5S00  Computer,  the  standard  deviation  between  the  five 
PF  integrals  and  the  five  Martree-Fock  integrals.  The  results  are 
summarized  in  Tables  XIII  and  XIV.  These  expansions  will  now  be  described 
in  detail. 

Set  1:  Switched  Input.  This  set  was  actually  set  3,  but  an 
accidental  shuffling  of  data  cards  exchanged  the  a ^  and  parameters. 

The  result  was  that  the  first  and  second  Slater  orbitals  almost  cancelled. 
In  other  words,  this  set  represents  an  extremely  bad  initial  expansion 
and  should  consequently  yield  relatively  bad  results.  However,  as 
Table  XIII  shows,  the  PF  procedure  was  able  to  recover  from  this  extremely 
bad  set  (note  the  initial  standard  deviation  of  8.983!)  and  produce  a 
respectable  final  standard  deviation  of  0.069!  Although  this  final 
standard  deviation  was  the  highest  found,  it  is  not  bad  and  proves  the 
great  stability  of  PF. 

Set  2;  Optimized  Set.  Since  the  4s  wave  function  has  one 
very  prominent  hump,  a  good  starting  expansion  might  be  obtained  by  first 
optimizing  the  fit  between  a  three-term  expansion  and  the  single  hump. 

The  result  of  this  optimization  is  given  in  Table  XIII.  Although  there 
are  no  nodes  in  this  set,  the  PF  procedure  inserted  one. 

Set  3;  Burns  Set.  This  set  has  its  exponents  determined  from 
41 

the  Burns  screening  constants.  The  A^  parameters  are  arbitrarily  fixed. 

Set  4:  Mann’s  Set.  This  set  has  its  exponents  determined  from 
38 

Mann's  binding  energies.  The  A^  parameters  are  arbitrarily  fixed. 

Three  of  these  three-term  expansions  resulted  in  two  node  PF 
wave  functions.  This  seems  to  indicate  that  the  PF  wave  function  is 
using  one  Slater  orbital  to  reproduce  the  main  hump  of  the  4s  wave  function 
(and  thereby  the  F°  and  normalization  integrals)  and  two  Slater  orbitals 
to  accomodate  the  orthogonality  constraints. 
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Set  5:  Mann *8  Set.  This  five-term  expansion  is  determined 
from  Mann's  binding  energies  as  was  the  case  for  Set  A.  It  is  Interesting 
that  although  the  initial  standard  deviation  is  worse  than  that  of  Set  A, 
the  final  standard  deviations  are  almost  identical.  The  five  Slater 
orbitals  for  this  and  later  sets  are  chosen  so  that  there  is  a  Is,  2s, 
and  3s  Slater  orbital  for  the  orthogonality  conditions,  a  As  orbital 
for  the  normalization,  and  a  3d  orbital  to  help  form  the  main  hump  of 
the  wave  function. 

Set  6:  a  Peaked.  The  location  of  the  maximum  value  of  a 
_ m _ 

Slater  orbital  is  given  by 


x  ■  p  /a 
max  rm  m 


(28) 


where  p  is  the  power  of  x  (or  r)  in  the  Slater  orbital  (see  Eq.  (20) >. 
m 

The  x  of  the  Hartree-Fock  orbitals  which  these  Slater  orbitals  are 
max 

approximating  are  given  in  Table  XV.  The  are  determined  by  setting 

p_  to  the  conventional  value  of  £  +  1.  The  A  in  all  the  five-term 

m  th 

expansions  are  arbitrarily  set  to  the  x  value  of  the  n  orbital 

max 

divided  by  the  square  root  of  the  sum  of  the  squares  of  all  the  x 

max 

values.  The  signs  are  selected  to  give  three  nodes.  This  procedure 
yields  the  same  accuracy  as  the  last  two  sets. 

Set  7:  p  Peaked.  Instead  of  letting  p  be  l  +  1,  the  a 
_ _m _  m  m 

could  be  set  to  those  values  of  Set  5  and  p  could  be  found  from  Eq.  (28). 

m 

The  result  is  given  in  Table  XV.  These  "exact"  p  can  then  be  rounded 

m 

to  the  nearest  Integer  value,  and  the  a  solved  for.  Thus  s  Slater 

m 

orbitals  approximate  the  Is  and  3d  Hartree-Fock  orbitals  whereas  p 
Slater  orbitals  approximate  the  2s,  3s,  and  As  Hartree-Fock  orbitals. 

This  set  gives  a  standard  deviation  at  least  three  times  lower  than  the 
previous  sets.  This  illustrates  the  importance  of  having  each  Slater 
orbital  accurately  represent  the  main  hump  of  the  orbital  it  is  approximating. 
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Set  8;  Final  Set^  This  set  is  simply  a  continuation  of  Set  7 
and  illustrates  the  capability  of  PF  to  obtain  accurate  wave  functions. 

In  several  minutes  of  computing  time  the  standard  deviation  is  reduced 
from  0.0100  to  0.0002.  In  detail  the  errors  in  the  five  integrals  are 
given  by 
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(29) 


The  orthogonality  constraints  are  the  best  satisfied. 

The  final  result  of  Set  5  has  two  nodes  whereas  the  results 
of  the  last  three  sets  have  three  nodes.  The  most  accurate  sets  possessed 
three  nodes  although  this  is  clearly  no  guarantee  of  a  good  set. 

Thus  far  the  various  expansion  sets  have  been  compared  as  to 
how  well  they  reproduce  the  five  radial  integrals  of  (29).  However,  the 
real  question  is  how  well  they  reproduce  the  4s  Hartree-Fock  wave  function. 
The  answer  is  given  in  Table  XVI  where,  in  particular,  the  percent  standard 
deviation  between  the  PF  and  Hartree-Fock  wave  functions  is  given  in  the 
third  column.  It  is  quite  evident  that  a  good  fit  of  the  integrals  in 
(29)  does  not  guarantee  a  good  fit  of  the  wave  function.  Actually,  of 
the  first  six  sets,  the  Burns  Set  and  Optimized  Set  give  the  best  re¬ 
sults.  Of  course.  Set  7  (and  therefore  Set  8)  yields  the  best  wave 
function.  The  <r° >  integrals  describe  where  the  PF  expansion  gives  a 
poor  fit.  The  fits  are  usually  poorest  at  very  small  x,  very  good  in 
the  main  portion  of  the  wave  function,  and  gradually  become  poorer  in  the 

tail.  The  wave  function  derived  by  PF  Set  8  is  compared  to  the  Hartree- 
Fock  wave  function  in  Fig.  (41) .  This  figure  gives  ample  proof  that  the 
PF  method  works  and  works  well. 
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12.  CROUP  PROPERTIES 


Group  theory  has  successfully  elucidated  and  simplified  the 
angular  portion  of  atomic  wave  functions.  It  is  hardly  an  exaggeration  to 
state  that  much,  if  not  most,  of  our  knowledge  of  the  actinide  and  rare 
earth  atoms  has  come  from  their  group  properties.  However,  the  present 
status  of  the  radial  portion  of  atomic  wave  functions  is  quite  primitive  with 
respect  to  the  angular  portion  of  the  wave  function.  The  best  radial  wave 
functions  presently  obtainable  result  from  some  type  of  Hartree-Fock  calcula¬ 
tion.  This  type  of  calculation  is  very  difficult  and  time  consuming,  even 
with  the  use  of  modern  computers.  In  addition,  there  are  many  unsolved 
theoretical  problems  in  the  Hartree-Fock  procedure. 

The  greatest  weakness  of  Hartree-Fock  procedures  is  that  the  re¬ 
sulting  wave  function  has  no  known  transformation  properties  with  respect 
to  any  symmetry  groups  of  the  atom.  Because  the  angular  wave  functions  do 
transform  in  known  ways  with  respect  to  certain  groups,  many  selection  rules 
are  known,  and  even  more  important,  many  matrix  elements  are  known  which 
are  proportional  to  each  other.  This  latter  result,  a  consequence  of  the 
Wlgner-Eckart  theorem,  enables  one  to  calculate  the  values  of  one  set  of 
matrix  elements  and,  after  finding  a  single  constant  of  proportionality,  to 
write  down  directly  the  values  for  another  set  having  the  same  transforma¬ 
tion  properties.  Unfortunately,  since  the  transformation  properties  of  the 
radial  wave  functions  are  unknown,  such  a  labor-saving  simplification  cannot 
be  used  to  relate  radial  matrix  elements. 

The  mathematical  techniques  to  be  used  in  this  work  are  not  well 
defined.  Although  a  great  deal  of  effort  has  gone  into  the  general  prob¬ 
lem^,  it  is  not  certain  that  the  techniques  exploited  by  mathematicians  are 
the  best  techniques  to  apply  to  physical  problems.  Because  of  the  uncer¬ 
tainty  in  the  techniques  and  the  associated  inherent  difficulties,  the 
present  Investigation  is  confined  to  the  hydrogen  atom,  since  the  radial 
wave  functions  are  well  known  and  well  studied  in  this  case. 
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Although  the  hydrogen  atom  has  been  extensively  studied,  most  of 
the  group  theoretical  work  which  involves  radial  wave  functions  has  consisted 
of  studies  of  hydrogen  using  the  orthogonal  groups  0(4),  0(4,1),  etc.  How¬ 
ever  in  this  case,  the  radial  wave  function  has  not  been  studied  directly 
since  the  basic  functions  for  these  groups  are  composed  of  product  of 
radial  and  angular  parts.  This  report  studies  the  group  properties  of  the 
hydrogen  radial  wave  function  themselves. 

Two  approaches  for  this  study  immediately  come  to  mind.  The 
first  approach  is  to  search  for  a  selection  rule  based  on  a  symmetry  group 
of  the  hydrogen  atom:  R(4),  the  proper  rotation  group  in  four  dimensions. 

The  second  approach  is  to  study  the  groups  for  which  Laguerre  polynomials 
form  bases  of  representation.  These  groups  may  be  pertinent  because  the 
hydrogen  radial  wave  function  is  proportional  to  a  Laguerre  polynomial. 

These  approaches  are  studied  in  detail  in  Appendix  5.  The  R(4) 
group  does  explain  some  selection  rules  and  has  promise  for  explaining 
some  vanishing  matrix  elements.  The  second  approach  leads  to  the  complex 
3-dlmenslonal  rotation  group.  In  order  to  search  for  the  desired  selection 
rules,  the  representation  to  which  the  operators  r  8  belong  must  be  known. 
The  latter  problem  is  not  trivial  and  the  solution  is  not  known. 

Another  approach  is  examined  in  Appendix  6.  Here  hydrogenlc 

radial  wave  functions  are  shown  to  form  bases  for  representations  of  the 

algebra  of  a  three  dimensional  non-compact  group.  Powers  of  r  are  found 

to  transform  as  tensors  with  respect  to  this  algebra.  Through  the  proper- 

-N 

ties  of  the  group,  matrix  elements  of  r  (N  positive)  can  be  expressed  in 
terms  of  the  familiar  Clebsch-Gordan  coefficients  for  R(3).  For  the  first 
time  these  coefficients  explain  the  known  selection  rules  of  the  radial 
wave  functions  of  hydrogen. 
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13.  SUMMARY 


The  calculation  of  atomic  spectra  requires  the  knowledge  of 
radial  wave  functions  which  determine  the  transition  probabilities 
between  the  atomic  energy  levels.  Transition  probabilities  are  very 
difficult  to  measure  and  they  are  therefore  calculated  whenever  possible. 
For  example,  in  arc  discharge  emission  calculations  the  angular  wave 
functions  are  calculated  in  intermediate  coupling  from  experimentally 

determined  Slater  integrals  and  the  radial  wave  functions  are  calculated 

4  18 

from  the  Thomas-Fermi  or  the  Bates  and  Damgaard  coulomb  approximation, 

but  these  latter  methods  are  limited  to  light  atoms  with  one  or  two 

valence  electrons.  However,  it  should  be  possible  to  determine  the  radial 

wave  function  for  any  atom  by  fitting  a  power  series  expansion  with 

parametric  coefficients  of  the  radial  wave  function  to  the  experimentally 

determined  Slater  integrals  (used  for  the  angular  wave  function)  of  the 

radial  wave  functions,  keeping  the  radial  wave  function  normalized  and 

orthogonal  to  closed  shell  electrons.  This  parameter  fitting  (PF) 

technique  has  the  advantage  of  being  independent  of  the  complex  electronic 

shell  structure.  Thus,  PF  should  succeed  where  all  other  methods  fail. 

An  alternative  to  the  PF  method  is  to  apply  group  theory  directly 
to  the  calculation  of  the  transition  probability  matrix  elements,  thus 
bypassing  the  need  for  an  exact  knowledge  of  the  rfedial  wave  function. 

This  technique  has  worked  extremely  well  for  angular  wave  functions  but 
thus  far  has  failed  for  radial  wave  functions. 

This  study  of  radial  wave  functions  has  progressed  with  some 
significant  progress.  The  PF  technique  has  been  developed  and  tested 
on  a  Hartree-Fock  Cu  4s  radial  wave  function  which  was  reproduced  by  PF 
to  within  +8%.  This  technique  is  relatively  slow  but  for  the  first 
time  a  method  is  available  for  the  determination  of  radial  wave  functions 
for  any  atom.  Moreover,  this  technique  is  more  consistent  than  other 
methods  since  the  same  Slater  integrals  are  used  for  both  the  angular 
and  radial  wave  functions. 
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The  knowledge  of  "experimentally  determined"  radial  wave  func¬ 
tions  via  PF  is  intriguing.  Transition  probabilities  may  be  calculated 
more  accurately  and  for  more  complex  atoms  than  presently  done  in  arc  dis¬ 
charge  calculations  such  as  described  in  the  first  six  chapters  of  this 
report.  Furthermore,  these  transition  probabilities  would  be  "experimentally 
determined",  i.e.,  not  from  ab  initio  theoretical  calculations. 

A  potentially  more  interesting  application  of  PF  is  to  determine 

the  radial  wave  functions  of  atoms  in  crystals.  Whenever  the  crystal 

field  parameters  are  known,  they  can  be  treated  as  additional  Slater  integrals 

in  PF.  Thus,  the  radial  wave  function  of  an  atom  can  be  studied  in  differ- 

3+ 

ent  crystals,  e.g.,  Nd  in  YAG,  LaCl^,  YjO^,  Y2®2^’  etc‘  correlati°n 

between  the  physical  extent  of  the  radial  wave  function  and  the  fluorescing 
or  lasing  properties  of  the  atom  can  then  be  studied. 

Significant  progress  has  also  been  achieved  in  the  application 
of  group  theory  to  the  calculation  of  matrix  elements  over  the  radial 
wave  functions  of  hydrogen.  In  particular^  hydrogenic  radial  wave 
functions  have  been  shown  to  form  bases  for  representations  of  the  algebra 
of  a  three  dimensional  non-compact  group.  Through  the  properties  of  this 
group,  matrix  elements  of  r  (N  positive)  were  expressed  in  terms  of 
the  familiar  Clebsch-Gordan  coefficients  for  the  proper  rotation  group 
in  three  dimensions,  R(3).  For  the  first  time  these  coefficients  explain 
the  known  selection  rules  of  the  radial  wave  functions  of  hydrogen. 

This  successful  group  theoretical  investigation  of  the  hydrogen 
radial  wave  function  is  of  particular  long-range  significance  and  should 
be  immediately  applicable  to  heavier  atoms  via  at  least  two  methods:  the 
quantum  defect  method  and  the  1/Z  expansion  method.  Both  methods  use 
hydrogenic  radial  wave  functions.  The  PF  method  could  also  be  used  by 
replacing  the  Slater  orbitals  in  the  PF  expansion  with  hydrogenic  wave 
functions. 
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This  successful  investigation  strongly  suggests  that  group 
theoretical  investigations  of  more  complicated  radial  wave  functions  may 
be  possible.  In  addition,  it  is  particularly  gratifying  to  demonstrate 
that  the  radial  and  angular  parts  of  the  atomic  wave  function  can  be 
placed  on  a  more  even  footing  and  theoretically  treated  with  the  same 
general  techniques. 
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A  Comparison  of  Lamps  for  Use  in  High 
Continuous  Power  Nd:YAG  Lasers 

I.  Liberman  and  R.  L.  Grassel 


I  lie  efficiency  of  krypton,  miioii,  anil  tungsten  -inilinc  lamps  for  |mni|ihig  Nil:YA(i  are  compared  hot  It 
.sped  rosi'opieally  and  liy  lasei  pumping.  The  rolat  ivo  pumping  efficiencies  of  xenon  and  tungsten  to  kryp¬ 
ton  are  0,8(1  and  1,1,  respectively.  However,  liecause  the  tungsten  filament  uses  a  higher  per¬ 
centage,  of  its  power  in  overcoming  the  laser  threshold,  the  krypton  lamp  with  its  larger  power  density 
yielded  the  liest  over-all  efficiency.  An  output  power  of  Itl'i  W  was  olitained  at  an  input,  power  of  !f(i  10  W 
for  an  over-all  efficiency  of  2.0%. 


Introduction 

The  tungsten  filament  lamp  was  the  first  and  is  today 
the  most  common  lamp  used  for  pumping  NdiYAG 
in  the  cw  mode  of  operation.  Because  of  the  develop¬ 
ment  of  the  tungsten-iodine  cycle  these  lamps  are  com¬ 
pact  and  can  operate  for  appreciable  duration  near  the 
melting  point  of  tungsten.  Even  at  this  temperature, 
the  radiance  of  tungsten  is  low  compared  to  some  gas 
discharges.  However,  if  the  radiance  from  continuous 
discharges  is  to  compete  with  the  radiance  from  a 
tungsten  filament,  lamps  must  be  constructed  having 
electrodes  and  envelopes  which  can  withstand  the 
average  power  dissipated  at  high  current  densities. 

The  diameter  of  Nd:YAG  rods  used  for  cw  applica¬ 
tions  is  limited  to  the  order  of  1  cm  because  of  the 
radial  them  ml  gtudLlits  ucUlTIfilig  at  li'igti  aWTagv. 
power  and  because  of  growth  considerations.  In 
nut  imttgfag  pump  CMcUf.  Ihsu  is  w’vj'tu#’  lo 
lamp  diameters  being  larger  than  rod  diameters  since 
the  imaged  radiance  can  never  be  greater  than  the 
source  radiance  and  that  it  is  desirable  to  collect  the 
radiation  over  4ir  sr.  At  current  densities  obtainable 
continuously  in  a  1-cm  diam  discharge,  the  plasma  is 
optically  thin.  That  is,  the  radiation  from  the  contin¬ 
uum  is  much  less  than  that  of  a  blackbody  at  the  same 
temperature  as  the  plasma.  The  radiance  from  the 
line  spectrum  may  be  orders  of  magnitude  greater  than 
that  from  the  continuum.  Therefore  if  the  line  spec¬ 
trum  of  a  gas  has  a  good  overlap  with  the  absorption 
spectrum  of  Nd:Y'AG  it  might  be  superior  gas  for  use  as 
a  laser  pump.  However,  although  intense,  the  lines 
are  usually  narrow  and  therefore  may  only  contain  a 
negligible  part  of  the  energy  dissipated  by  the  lamp. 
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It  has  Iteeii  observed1  3  that  the  line  spectrum  from 
krypton  is  a  lie t ter  match  to  NdiYAG  than  the  line 
spectrum  of  xenon.  For  pulsed  pumping  of  lasers, 
xenon  is  superior  to  krypton  because  the  current  densi¬ 
ties  are  sufficiently  high  that  the  line  spectrum  is 
negligible1  and  xenon,  because  of  its  large  atomic 
weight,  is  a  more  efficient  continuum  radiator  than 
krypton.5  For  continuous  discharge  lamps  the  com¬ 
parison  becomes  difficult  because  both  the  continuum 
and  the  line  radiation  may  he  nearly  equivalent.  A 
number  of  methods  can  he  used  to  compare  lamps  al¬ 
though  none  of  them  are  entirely  satisfactory. 

The  simplest  method  is  to  measure  the  fluorescence 
from  a  Nd:YAG  sample  when  irradiated  by  various 
SOW c».s  pluVCil  ; iV  lulgb  >1  isYuml'S  ii'UIil  \  AV )  suTtqm,-. 
This  experiment  is  difficult  to  perform  because  the 
frrtffwily  «  W,  Y'l.i,  VTtttl  tmlfertS  m  is.  > /triers  si  tnagwr 


tilde  greater  than  the  fluorescence,  and  the  radiation 
around  1  .00  is  comparable  to  the  fluorescence.  Read 1 
performed  an  experiment  of  this  type  and  concluded 
that  the  relative  fluorescence  resulting  from  Kr,  Xe, 
and  W-l2  lamps  was  I. DO:  1.00:0.63,  respectively. 

A  second  approach  is  based  on  spectral  measure¬ 
ments.  While  this  technique  yields  a  great  deal  more 
information,  the  interpretation  is  difficult,  tedious,  and 
requires  very  precise  wavelength  measurements.  The 
importance  of  wavelength  precision  is  a  result  of  the 
narrow  absorption  lines  of  Nd+3  in  YAG  (about  10  A) 
and,  as  will  he  shown,  the  even  narrower  line  radiation 
from  Xe  and  Kr.  Previous  spectral  measurements2-4 
of  continuous  noble  gas  discharges  have  been  made 
using  relatively  poor  resolution.  In  one  comparison  be¬ 
tween  Xe  and  Kr2  the  resolution  stated  was  f)0.0  A. 
However,  the  band  replotted  curves  show  full  widths 
at  half  height  of  over  100  A  which  is  over  an  order  of 
magnitude  greater  than  the  actual  line  width  and  thus 
distorts  the  actual  contribution  from  the  lines.  There- 
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Fig.  1.  Linn  spectrum  of  (a)  krypton  and  (b)  xenon  lamp  at.  200 
mV  full  scale.  Spectral  radiance  =  I pWfsr  -  A  —  mm1  — 


fore  before  spectral  measurements  can  be  used  to  com¬ 
pare  gases,  high  resolution  must  be  obtained. 

A  third  method  of  comparison  is  the  actual  laser  test. 
Here  comparison  is  difficult  unless  the  laser  head  is 
optimized  for  each  specific  lamp.  Since  a  tungsten 
filament  lamp  has  well  defined  radiating  dimensions 
while  a  gas  discharge  does  not,  this  comparison  can  be 
subject  to  criticism.  Read  found  that  Kr  had  a  slope 
efficiency  about  2.4  times  as  great  as  Xe  or  W-Is.  How 
ever,  since  no  absolute  numbers  arc  given,  the  degree 
of  optimization  cannot  be  determined. 

Experimental  Apparatus 

The  discharge  lamps  used  for  our  tests  are  very 
similar  to  the  lamp  described  by  Stcarn  and  Colliver.* 
That  is,  cooling  water  not  only  flowed  over  the  quartz 
envelope,  but  also  circulated  inside  the  hollow  elec¬ 
trodes.  The  electrode  spacing  was  5  cm  and  the  bore 
diameter  .r>..r)  mm.  The  electrode  diameter  was  5  mm 
so  that  the  gas  space  behind  the  electrode  tips  was  kept 
at  a  minimum  thus  building  up  the  operating  pressure 
as  much  as  possible  over  the  fill  pressure  of  2  atm. 

The  rod  used  was  5  mm  X  50  mm  and  had  antire¬ 
flection  coatings  placed  on  the  flat  ends.  The  rod  was 
centered  in  an  optical  cavity  of  41-cm  length  terminated 
by  two  dielectric  mirrors  of  1-m  radius  of  curvature. 
The  output  reflectivity  was  9 5%  which  for  Kr  and  Xe 
pumping  was  found  to  be  more  nearly  optimal  than 
98%  or  92%  reflectivity.  The  pump  cavity  consisted 
of  a  25-cm  diam  gold  coated  sphere.7  While  we  believe 
this  cavity  is  superior  to  an  elliptical  cylinder,  it  has  a 
further  advantage  when  comparing  lamps  in  that  it 
provides  a  one  to  one  image.  Therefore  good  geo¬ 
metric  coupling  is  always  assured  if  the  lamp  diameter 
is  equal  to  or  smaller  than  the  rod  diameter.  In  an 
elliptical  cylinder  having  an  appreciable  eccentricity, 
the  image  coming  from  over  J  the  solid  angle  of  the 
lamp  radiation  is  magnified. 

Spectral  measurements  were  performed  using  a  Jar- 
rell-Ash  £-m  Ebert  monochromator.  The  curved  slits 
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were  40  n  wide  corresponding  to  a  resolution  of  0.64  A. 
The  detector  used  was  a  silicon  p,i,n  photodiode.  The 
system  was  calibrated  with  a  calibrated  tungsten  spectral 
radiance  standard  and  automatically  corrected  for  non¬ 
linearity  between  5000-10,000  A  by  using  a  variable  at¬ 
tenuator. 

The  laser  power  was  measured  with  a  Spectra-Physics 
model  401  photovoltaic  cell  in  conjunction  with  two 
quartz  plate  beam  splitters.  It  was  calibrated  using  a 
TRG  model  107  ballistic  thermopile.  Each  quartz 
plate  deflects  the  radiation  7.5°  from  the  normal,  and 
in  mutually  perpendicular  planes.  This  allows  for 
relatively  easy  alignment  and  removes  any  calibration 
variations  with  polarization.  The  discharge  lamps 
were  dc  excited,  while  the  tungsten  lamp  ran  on  60-Hz 
ac. 

Spectral  Measurements 

The  spectra  of  the  krypton  and  xenon  lamps  operat¬ 
ing  at  35  A  are  shown  in  Figs.  1  and  2.  Since  the  peak 
intensity  of  the  continuum  is  orders  of  magnitude  lower 
than  that  of  the  lines,  the  continuum  radiation  is 
clearly  optically  thin.  Line  widths  were  obtained  by 


Fig,  2.  Continuum  spectrum  of  (a)  krypton  Inmp  at  2  (A  full 
scale  and  (l>)  xenon  lamp  at  5  nV  full  scale.  Spectral  radiance  — 
ir>.!t/uW(sr  —  A  —  mm1  —  /u  V ) -•. 


Fig.  3  Power  conversion  of  fi  mm  X  •r>0  mm  Nd:YA(5. 
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Tabic  I.  Comparative  Spectral  Efficiency  of  Lamps  for  Pumping  Nd:YAG 


Wavelength 

region 

M 

Useful 

portion 

(A) 

Weighting 

factor 

2.28-kW  krypton 

Avg. 

intensity  product. 

(mV)  (A-sV) 

2.73-kW 

Avg. 

intensity 

(mV) 

xomin 

product 

(A-„v> 

2.4-kW  tungsten 

Avg. 

intensity  product. 

(mV)  '  (A-/.V) 

0.50-0.60 

300 

0.50 

0.8 

120 

1.4 

210 

1(1 

150 

0.74-0  83 

600 

0.75 

0  6 

270 

1.0 

450 

1  9 

855 

0.86-0  90 

30 

0  85 

0.5 

13 

1.0 

25 

I  9 

48 

Lines 

2° 

0  75 

3(H)* 

450 

— 

— 

- - 

— 

Total 

853 

685 

1053 

Total /kW 

374 

250 

430 

Fluorescent  efficiency 

4.3% 

2  9';, 

5.0% 

°  Average  line  width  of  krypton. 
1  Total  absorption  of  all  lines. 


using  ft  much  slower  sctui  speed.  For  each  gits  the 
line  width  of  the  various  lines  varied  about  20%.  An 
average  measured  full  width  at  half  height  for  krypton 
was  2.0  A  and  for  xenon  it  was  4.9  A.  Therefore, 
despite  the  lower  peak  intensities  and  fewer  number  of 
lines,  the  xenon  lamp  is  a  superior  emitter  of  line  radi¬ 
ation  in  the  wavelength  region  studied.  In  the  con¬ 
tinuum  the  xenon  lamp  is  also  more  efficient  than  kryp¬ 
ton  by  about  30%. 

Because  of  the  spectral  overlap  of  the  krypton  lines 
with  the  Nd:YAG  absorption  spectrum  the  krypton 
lamp  can  be  more  efficient  than  a  xenon  lamp  for  pump¬ 
ing  Nd:YAG.  The  overlap  was  measured  by  placing 
a  spectroscopic  Nd:YAG  sample  in  front  of  the  mono¬ 
chromator  entrance  slit.  The  results  show  that  the 
overlap  for  krypton  lines  is  appreciable  while  for  xenon 
the  overlap  is  negligible.  As  shown  in  Table  I  an 
estimate  of  the  relative  pumping  efficiency  of  the  lamp 
is  given  based  on  the  spectral  measurments.  The  effect 
of  reflection  from  the  Nd:YAG  spectroscopic  sample 
was  taken  into  account  in  making  the  calculation.  It 
was  assumed  that  all  photons  absorbed  in  the  measured 
wavelength  region  contribute  equally  to  photons  emit¬ 
ted  at  1.064  n.  However,  the  energy  absorbed  at 
each  wavelength  was  multiplied  by  the  ratio  of  the 
wavelength  to  1.064  g.  Since  the  spectral  measure¬ 
ments  are  absolute  numbers,  a  calculation  of  laser 
efficiency  can  be  approximated.  This  was  done  as 
given  in  the  last  row  of  Table  I.  The  collected  energy 
was  assumed  to  be  0.75  (3ir  sr)  of  the  emitted  radiance 
based  on  a  computation  of  the  spherical  pump  cavity 
efficiency.  Due  to  the  difficulty  in  determining  the 
spectral  overlap  the  accuracy  of  the  over-all  efficiency 
is  no  better  than  30%.  The  tungsten  measurements 
were  obtained  using  a  2550  K  spectral  radiance  standard 
and  converted  to  3200  K  by  using  a  table  of  blackbody 
radiance.  The  2.4-kW  input  power  was  calculated  as 
the  radiation  from  a  5  mm  X  50  mm  cylinder  of  0.4 
emissivity  at  3200  Ii  having  an  80%  radiative  efficiency. 
These  calculations  show  that  tungsten  is  about  16% 
superior  than  krypton  in  spectral  efficiency  for  pumping 
Nd:YAG. 


Laser  Tests 

One  lamp  was  used  for  both  the  xenon  and  the 
krypton  tests  so  that  all  other  parameters  would  be 
unchanged.  The  lamp  was  first  filled  with  two-afin 
xenon  and  pinched  off.  After  the  test,  only  a  negli¬ 
gible  amount  of  blackening  was  observed  on  (he 
quartz  wall.  The  lamp  was  opened,  refilled  with  2-atm 
krypton  and  pinched  closed.  For  the  case  of  tungsten 
pumping,  no  lamp  was  available  for  a  proper  geo¬ 
metric  match  to  the  laser  rod.  A  2  mm  X  50  mm 
single  coil  filament  was  used  which  limited  the  input 
power  to  1070  W.  The  results  are  shown  in  Fig.  3. 
The  saturation  of  the  output  under  xenon  pumping  is 
not  peculiar  to  the  xenon  lamp.  Saturation  was  ob¬ 
served  with  the  krypton  lamp  when  the  resonator  spac¬ 
ing  was  increased  to  50  cm.  It  is  believed  that  the 
saturation  is  a  result  of  a  strong  lens  effect  in  the  rod 
due  to  the  radial  thermal  gradient  present.  In  pump¬ 
ing  the  rod  without  the  resonators,  and  placing  a  He:Nc 
beam  axially  through  the  rod,  a  positive  lens  was  ob¬ 
served  with  a  waist  minimum  occurring  loss  t  han  20  cm 
from  the  rod.  This  strong  lensing  should  increase  the 
losses  for  many  of  the  oscillating  modes.  It  is  further 
hypothesized  that  since  the  xenon  lamp  is  more  effi¬ 
cient  than  the  krypton  lamp  in  pumping  the.  visible 
absorption  bauds,  more  heat  is  dissipated  in  the  rod  for 
the  same  output,  power.  Therefore,  the  equivalent 
lensing  should  occur  at  lower  output  power  levels  for 
the  xenon  lamp.  With  tungsten  pumping,  because  of 
the  low  input  power,  a  98%  output  mirror  produced  an 
over-all  efficiency  of  1.4%  compared  with  1.23%  using 


Table  II.  Summary  of  Pumping  Efficiency  of  Lamps 
Calculated  by  Several  Ways 


Efficiency 

Kr 

Absolute  % 
Xe 

W 

Kr 

It  dative. 

Xe 

W 

Spectral 

4.3 

2.9 

5.0 

1  <1 

0  08 

1 . 16 

Slope 

3.74 

3.22 

— 

1  0 

0.7(1 

— 

At  80  W 

2.03 

2.13 

— 

1.0 

0  81 

— 

At  15  W 

1.16 

0.89 

1.22 

1  (1 

0  77 

1  (HI 
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the  95%  mirror.  However,  for  comparison  purposes 
the  95%  mirror  results  arc  shown. 

Because  of  the  threshold  power  required  to  obtain 
oscillations,  the  over-all  efficiency  of  the  system  is  a 
function  of  input  power  thus  making  comparisons  sub¬ 
ject  to  interpretation.  Assuming  the  spectral  effi¬ 
ciency  of  the  lamp  does  not  vary  with  power,  the  slope 
efficiency  should  be  proportional  to  the  spectral  effi¬ 
ciency  and  equal  to  it  for  a  lossless  system.  Row  2  of 
Table  II  gives  these  results.  Row  8  and  4  of  Table  II 
gives  flic  over-all  efficiencies  at  outputs  of  80  W  and 
15  W. 

Discussion  and  Conclusions 

The  comparison  of  efficiency  by  spectral  and  laser 
tests  agree  within  experimental  error.  We  believe 
that  because  of  the  difficulty  in  interpreting  the  spec¬ 
tral  data,  the;  slope  efficiency  is  the  most  accurate  mea¬ 
surement.  Therefore  our  best  estimate  of  relative  effi¬ 
ciency  between  Kr,  Xe,  and  W  is  1:0.80:1.1,  re¬ 
spectively.  This  is  in  considerable  disagreement  with 
Read.1  Our  spectral  data  look  quite  different  from 
those  previously  published2  4  because  their  spectral 
resolution  is  much  poorer  than  the  true  line  width  of 
xenon  and  krypton  operating  under  the  conditions 
described. 

Using  the  krypton  lamp  an  output  power  of  105  W  at 
an  efficiency  of  2.92%  has  been  obtained.  A  potas¬ 
sium-mercury  spectral  iulditive  lamp  has  been  de¬ 
veloped”  which  has  a  spectral  efficiency  over  three 
times  that  of  tungsten  which  in  turn  has  a  better 


spectral  efficiency  than  krypton.  Therefore  we  pre¬ 
dict  that  with  development  a  ew  Nd.  YAG  laser  will 
have  an  efficiency  appreciably  greater  thau  the  results 
reported  here. 

We  wish  to  thank  C.  Spontak,  E.  Kowal,  and  A.  Zot- 
tola  for  constructing  the  electrodes  used  in  the  dis¬ 
charge  lamps.  We  also  thank  C.  H.  Church  for  de¬ 
veloping  the  foundat  ion  of  this  work. 

Xolr  wtilat  in  proof:  Additional  (lata  on  rare  gas  pumping  ef¬ 
ficiencies  has  been  recently  published  by  J.  It.  Oliver  and  F.  S. 
Barnes,  I 1SK1S  I.  Quant, tun  Flection.  QK5, 225,  2:’,2  ( HUH*)-  As  in 
lief.  I,  laser  test,  results  are  in  arbitrary  units  which  could  lead 
to  erroneous  conclusions  if  all  parameters  are  not  full v  opti¬ 
mized.  Also,  as  in  Itef.  2,  the  spectral  resolution  was  ItMt  A  which 
is  estimated  to  he  an  order  of  magnitude  greater  than  the  actuul 
line  width  under  their  operating  conditions.  However,  their 
conclusions  are  similar  to  ours. 
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APPENDIX  2 


ROD  TO  SLAB  TRANSFORMATION 

Al' hough  the  simple  slab  laser  geometry  has  been  used  for 
the  computer  calculations  reported  in  the  text,  the  most  common  laser 
geometry  is  that  of  a  cylindrical  rod  in  a  cavity.  In  order  to  fully 
utilize  the  computed  results  for  the  slab  laser  geometry,  the  cylindrical 
rod  in  a  cavity  geometry  must  be  transformed  to  an  equivalent  slab  laser 
geometry.  This  transformation  is  given  in  this  appendix. 

The  emittance  (power  per  unit  surface  area)  absorbed  by  an 
arbitrary  body  as  a  function  of  wavelength  is  given  by 


(2.1) 


Pfl  =»  Pq  /  (1  -  exp  (-ac)]  du 


where  a  is  the  spectral  absorptivity,  c  is  the  chord  length  traversed 
by  the  light  in  the  body,  Pq  is  the  incident  light  emittance,  and  w 
is  the  solid  angle.  The  integration  is  over  the  entire  absorbing  body. 
An  effective  chord  length  c  can  be  defined  as  that  length  which  absorbs 
from  a  single  light  ray  the  same  power  density  as  the  actual  body,  i.e., 

(2.2)  P  =  P  U  -  exp  (-ctc)J 

3  o 

Thus  an  effective  chord  length  can  be  found  for  a  given  wavelength  for 
any  absorbing  body  from  an  evaluation  of  Eq.  (2.1). 

The  effective  chord  length  for  the  slab  geometry  is  easily 
found  when  the  light  is  perpendicularly  incident  on  the  slab  because 
then 


(2.3) 


P  =  P  [1  -  exp  (-a  t) ] 
a  o  r  s 
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where  t  is  the  slab  thickness  and  ag  is  the  absorptivity  of  the  slab. 
Therefore 

(2.4)  c  ■  t 

s 

which  is  rather  obvious. 

The  effective  chord  length  for  a  laser  rod  geometry  can  be 
found  from  Eq.  (2.1)  with  one  modification.  The  rod  in  a  spherical 
laser  cavity  subtends  a  very  small  solid  angle  when  viewed  from  the 
lamp  and  therefore  almost  all  the  light  absorbed  by  the  laser  rod  is 
first  reflected  by  the  laser  cavity  walls.  The  rod  is  thus  bathed  in 
this  reflected  lamp  light.  However,  the  light  intensity  emitted  by 
the  lamp  has  an  approximate  Lambertian  distribution  with  respect  to 
the  angle  0  between  the  emitted  light  ray  and  the  normal  to  the  outer 
lamp  wall.  Specular  reflecting  laser  cavities  maintain  this  Lambertian 
distribution  at  the  rod  surface.  A  surface  roughened  rod  will  slightly 
modify  the  Lambertian  distribution  but  probably  not  by  a  factor 
larger  than  the  accuracy  of  the  general  pump  efficiency  calculations. 

In  any  event,  we  assume  the  light  distribution  within  the  rod  to  be 
Lambertian  and  that  Eq.  (2.1)  can  be  written  as 


(2.5) 


■  (Pq/tt)  /  [1  -  exp  (-a^c)]  cos  0  dw 


where  the  integration  is  over  the  laser  rod  and  is  the  rod  absorptivity. 
The  7T  in  the  denominator  keeps  the  same  light  emittance  incident  on  the 
laser  rod  and  slab.  The  integral  in  Eq.  (2.5)  has  not  been  analytically 
solved  but  Pukhov  has  numerically  evaluated  it  for  an  infinite  cylinder 
of  diameter  d  and  expressed  his  results  in  the  simple  form: 


(2.6) 


c 

r 


d+6.44 

r 


0.27) 


d 
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The  rod  to  slab  transformation  is  completed  by  equating  the 
emittance  absorbed  by  the  slab  to  the  emittance  absorbed  by  the  rod, 
provided  the  incident  emittance  is  the  same.  In  other  words  Eq.  (2.2) 
is  the  same  value  irrespective  of  which  effective  chord  length  is  used, 
i. e. , 


(2.7) 

which  implies  that 


a  c  =  a  c 
s  s  r  r 


(2.8)  t  -  D  cr 

where  D  is  the  activator  doping  of  the  rod  relative  to  the  activator 
doping  of  the  slab. 

The  transformation  is  strictly  valid  for  only  one  wavelength 
but  for  the  5  x  50  mm  Nd:YAG  rod  used  for  the  laser  system  efficiency 
measurements  Eq.  (2.8)  is  fairly  constant  for  all  wavelengths,  as  can 
be  seen  by  using  the  experimental  parameters  to  obtain 


(2.9) 


D  -  1.3/0. 7 


and 


(2.10)  0.63  cm  <  t  <  0.93  cm 

where, using  Eqs.  2.6  and  2.8, the  maximum  absorptivity  (10  cm  ^)  within 
the  main  absorption  bands  of  the  rod  determines  the  lower  limit  and 
zero  absorptivity  determines  the  upper  limit.  An  average  value  for 
the  absorptivity  is  roughly  2  cm  ^  which  gives  an  average  value  of 

(2.11)  t  =  0.84  cm 
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This  value  is  better  than  one  might  expect  at  a  causal  glance.  The 
average  value  of  2  cm  1  for  the  absorptivity  is  admittedly  crude  but 
only  appears  in  the  denominator  of  Eq.  (2.6)  and  a  +  1  cm  error  in  this 
average  value  only  gives  an  error  of  +  0.04  cm  for  the  effective  slab 
thickness.  Thus,  to  a  good  approximation,  the  laser  rod  used  in  the 
experimental  measurements  reported  in  Appendix  1  is  equivalent  to  a 
slab  of  thickness  0.84  cm. 
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APPENDIX  3 


NEWTON’S  METHOD  FOR  PARAMETRIC  FITTING 
OF  RADIAL  WAVE  FUNCTIONS  IN  INTEGRALS 

A.  J.  Federowicz 

The  radial  wave  fitting  problem  has  already  been  formulated 
as  one  of  solving  sets  of  simultaneous  non-linear  equations  in  which  the 
number  of  fitting  parameters  (variables  in  the  fitting  problem)  equal 
the  numbers  of  equations.  These  simultaneous  equations  result  from 
applying  normality  and  orthogonality  conditions  which  the  radial  wave 
function  must  satisfy  and  from  attempting  to  obtain  agreement  with 
certain  experimental  observations. 

If  the  number  of  equations  were  larger  than  the  number  of 
parameters,  this  proposed  procedure  would  be  modified,  but  not  drastically. 
What  would  then  be  done  would  be  to  fit  such  overdetermined  sets  of 
equations  using  least  squares.  As  will  be  seen,  a  least  squares 
criterion  is  suggested  here  for  a  different  reason  and  to  apply  this 
criterion  to  overdetermined  sets  of  equations  would  be  a  minor 
modification. 

The  proposed  procedure  is  to  use  a  modified  form  of  Newton's 
method  for  solving  systems  of  simultaneous  non-linear  equations.  The 
modifications  are  aimed  at  guaranteeing  convergence  of  the  iteration 
and  at  avoiding  the  use  of  analytic  partial  derivatives. 

Newton' s  method  for  solving  a  set  of  simultaneous  non-linear 

equations  is  as  follows:  as  the  problem  is  generally  phrased,  for  a 

set  of  m  functions  of  m  variables,  one  is  to  find  a  point  at  which  the 

m  functions  are  simultaneously  zero,  that  is,  find 

x  =  (x  *  x0w  ...  x  *)  such  that 
12  m 

(3.1)  fx(x*)  =  f2(x*)  =  ...  =  fm(x*)  =  0 

or  in  compact  form 

(3.2)  f(x*)  =  0. 
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Newton's  method  consists  of  making  simultaneous  linear  approximations 

-*(  n) 

to  the  m  functions  at  some  point,  x  ,  and  solving  this  set  of  linear 
approximations  for  the  point  where  all  will  be  zero.  That  is 


(3.3)  f^)  *  fl(x(n))  ♦  f  %  +  f  +  ...  ♦  ,  Ax 

2  m 

f  ( x )  —  f(x^n^)+f  Av+f  Ay  +  +f 

2W  r2vx  )  +  2  ^  +  t  ax2  *■  ...  +  Aj( 

a.  mm 


f  (x)  *  f  (XU))  +  f  AX  +  f  AX  +  ...  +  f  AX 

m  m  mx,  1  mx_  2  mx  n 


which  upon  being  set  equal  to  zero  yields 


or  in  compact  form 

(,.5)  5i(n)  =  -  t  fj(?(n))]-1 

The  iteration  then  consists  of  using  the  obvious  recurrence  relation 

(3.6)  ;("n)  =  2(n)  +  3i(n) 

If  there  is  only  one  equation  in  one  unknown,  this  formulation  reduces 
to  the  familiar  form  of  Newton’s  method  in  which 
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V 


(3.7) 


In  the  radial  wave  function  fitting  problem  some  of  the 
equations  will  not  be  in  analytic  form  but  only  as  numerical  approxima¬ 
tions  to  integrals.  Difference  approximations  to  the  partial  derivatives 
should  thus  be  used  in  constructing  the  [f^]  matrix.  This  modification 
will  considerably  simplify  the  necessary  computer  programming  since 

only  one  subroutine,  which  simply  evaluates  f,  will  be  necessary  in 
r  ^  -i 

setting  up  f  and  L  for  each  case.  It  is  necessary  to  use  differences 

x_i 

which  will  estimate  accurately,  however.  Attention  must  thus  be 
paid  to  the  question  of  how  many  significant  digits  result  from  the 
differencing  process. 

In  order  to  guarantee  convergence  of  the  iteration, 

computed  using  ( 1 . 6) ,  will  be  accepted  only  if  it  represents  a  better 

( n ) 

solution  to  the  set  of  simultaneous  equations  than  x  does.  For 
reasons  of  simplicity  "better"  will  be  taken  to  mean  that 


(3.8)  ||f(i(nrn))|!2  <  ||f(x(n))||2 

with  weighting  of  the  components  of  f  being  an  option.  Equation  (l,8) 
says  that  the  sum  of  the  squared  errors  must  be  reduced  from  one  step 
to  the  next. 

If  x(n+l)  results  in  an  increase  in  ||f(x)||^  then  the  chord 

joining  x^  and  x^1"1"^  will  be  searched  for  a  better  value.  This  can 

^  -»(n)  ] 

be  done  assuming  linearity  of  f  between  x  and  x  as  follows: 

(3.9)  f(t(n)+  s£(n))  i  f(x(n) )  +  s(f(x(n+l))  -  f(x(n))) 


=  f(n)  +  sSf(n) 
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The  point  at  which  the  norm  of  this  linear  approximation  is  minimized 
can  be  used  to  determine  s.  The  minimizing  s  is  given  by 


(3.10) 


s* 


f(n).  Af(n) 

||3f||2 


f  2 

If  this  interpolating  procedure  does  not  produce  a  smaller  llfll  at 
_k(n)  -*^(n)  -*(n) 

x  +  si*  Ax'  ,  then  it  may  be  repeated  several  times  using  x  and 

x(n)  +  9*  Ax^n^  as  new  end  points.  If  this  strategy  fails  then  a 

different  one  dimensional  search  strategy,  such  as  interval  halving, 

may  be  employed. 

—A 

In  the  case  that  f  =  0  has  no  solution,  as  is  likely  to  be 

the  case  in  this  problem,  the  determinant  of  [f^]  will  tend  to  zero  as 

| | f | |  tends  to  a  positive  value.  The  analogy  with  the  one  dimensional 

case  is  that  f'(x^)  is  zero  at  any  point  where  | |f^||  is 

minimized,  if  that  minimum  is  non-zero,  which  is  the  case  when  f  =  0 

has  no  solution.  It  is  instructive  to  examine  what  Newton's  method 

does  when  used  to  look  for  a  zero  of  f  which  does  not  exist  in  the  one 

dimensional  case.  The  method  will  usually  oscillate,  neither  approaching 

a  limit  nor  diverging,  in  such  a  case.  For  example,  if  the  problem  is 

2 

to  find  a  zero  of  f  =  1  +  x  then  one  has 


(3.11)  Ax 


(n)  _ 


1+x 


(nr 


2x 


TV 


and 

(3.12) 


,(n+l)  _ 


.(n) 


2x 


(n) 


For  large  |x^|  the  term  l/2x^  can  be  ignored  and  this  formula  says 

that  x(n)  will  be  halved  at  each  step.  As  Ix^l  becomes  small, 

/  \ 

however,  this  term  dominates  and  causes  |x'  to  again  become  large, 
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while  also  causing  x^  to  change  sign.  A  similar  behavior  occurs  in 

-1 

multi-dimensional  problems  in  which  f  =  0  has  no  solution.  By  insisting 

r\ 

that  jlf||  be  reduced  on  each  step,  this  possibility  of  oscillation  is 
avoided. 

-I 

A  problem  which  may  also  occur  is  that  f  =  0  may  have 
several  solutions  and/or  that  ||f||  may  have  several  local  minima. 

Since  this  problem  can  be  rigorously  resolved  only  in  the  case  that  one 
is  dealing  with  a  convex  function,  different  starting  points  should  be 
tried  for  the  iteration.  Then  if  the  iteration  consistently  converges 
to  the  same  answer  in  a  given  case,  one  can  be  reasonably  certain  that 
the  answer  is  unique;  if  convergence  to  different  answers  occurs,  then 
the  best  answer  should  be  sleeted  and  one  must  be  careful  about  finding 
the  best  answer  in  other  cases. 

As  a  final  point,  it  should  be  mentioned  that  a  number  of 
other  optimum  search  procedures  are  potentially  applicable  to  this 
problem.  These  include  Direct  Search,  GOOP  and  a  non-linear  least 
squares  program  written  by  R.  D.  Fardo,  three  methods  which  have  been 
used  extensively  at  this  location.  An  essential  feature,  which  must 
be  maintained  in  any  computer  program  of  this  proposed  procedure  in 
order  to  make  such  other  techniques  easily  usable,  is  the  separation 
of  the  functional  evaluation  subroutine  from  the  logic  of  the 
iterative  method. 
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APPENDIX  4 


COMPARISON  OF  DIRECT  SEARCH  AND  NEWTON'S  METHOD 


A.  J.  Federowicz 


This  appendix  is  a  post-mortem  on  why  Newton's  method,  suggested 
in  Appendix  3,  did  not  work  well  on  the  radial  wave  function  parametric 
fitting  problem.  Conversely  light  is  shed  on  the  question  why  Direct 
Search  did  a  relatively  much  better  job. 

The  characteristic  exhibited  by  the  modified  Newton's  method 
was  that  it  became  hung  up  in  local  areas  during  the  search.  Very  little 
improvement  in  the  least  squares  criterion  would  occur  for  many  iterations. 
At  the  time  that  this  occurred,  a  computer  program  bug  was  suspected.  A 
bug  was  not  the  cause,  however,  but  rather  the  problem  is  an  ill-conditioned 
one  for  Newton's  method.  As  an  example  of  slow  convergence,  the  following 
sums  of  squares  occurred  in  one  case: 

Iteration  #  Error  Sum  of  Squares 


0 

1 

2 

3 

4 


6.319567 

3.107108 

3.103914 

3.102053 

3.101231 


By  comparison,  Direct  Search  produced  the  following  sequence  starting  near 
the  parameter  point  of  iteration  #1: 
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Iteration  // 


Error  Sum  of  Squares 

0  3.108122 

1  2.538027 

2  .603644 

3  .338447 

4  .217680 

5  .116751 

• 

12  .012570 

There  are  a  number  of  equivalent  ways  of  explaining  the  observed 
phenomena.  One  explanation  is  that  the  radial  wave  function  parameters 
are  highly  correlated  in  their  effect  upon  the  sum  of  square  criterion. 

The  result  of  this  correlation  is  that  Newton's  method  selects  a  very  poor 
direction  for  minimizing  the  sum  of  squares  locally.  The  selected  direction 
would  be  very  good  if  the  linear  approximation  used  by  Newton's  method 
were  valid  far  from  the  current  point.  However,  the  problem  is  highly 
non-linear  and  the  use  of  a  linear  approximation  to  select  a  point  far 
from  the  current  one  is  thus  weak. 

In  order  to  illustrate  what  can  occur,  assume  for  the  moment 
that  there  are  two  functions,  f^  and  for  which  zeroes  are  being  sought 
and  two  parameters  x^  and  x^  which  are  being  searched  for  these  zeroes. 

Then  the  search  procedure  of  Appendix  3  consists  of  making  linear  approxi¬ 
mations  to  both  and  solving  for  where  these  linear  approximations  are 
simultaneously  zero.  A  case  in  which  the  two  parameters  are  correlated  in 
their  effect  upon  the  sum  of  squares  is  illustrated  in  the  Fig.  42.  The 
correlation  is  identified  by  the  nearly  parallel  gradients  Vf^  and  '^^2’ 

Based  upon  these  gradients,  a  good  local  direction  in  which  to  search  is 
the  X2  direction.  Based  upon  the  contours,  the  step  length  to  take  would 
be  approximately  the  length  of  the  pictured  Vf^.  However,  the  modified 
Newton's  method  selects  the  intersection  of  the  zeroes  of  the  linear 
approximations.  The  resultant  step  length  overtaxes  the  linear  approxi¬ 
mations  and  only  a  small  improvement  in  the  sum  of  squares  can  be  achieved 
in  this  direction. 
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A  measure  of  the  correlation  between  parameters  in  the  illustrated 
case  is  the  sine  of  the  angle  between  the  two  gradient  vectors .  This 
is  given  by 

i  Det  [Vf  Vf  ]| 
sin(Vf  _  ,Vf_)  =  - - - - - — 

1  2  IVfjl  l»fjl 


As  this  measure  increases  from  0  to  1,  the  gradient  vectors  go  from 
parallel  to  being  orthogonal. 

The  similar  statistic  in  3  parameters  measures  how  nearly  co- 
planar  the  3  gradient  vectors  would  be.  The  identity 


Det  [Vf x  Vf2  Vf 3] 


sin  (VJf1,Vf2)  sin  (VfjVf^Vf  ) 


puts  this  co-planar  relationship  into  terms  similar  to  the  2  dimensional 
case.  In  3  dimensions  the  measure  equals  the  product  of  the  sine  between 
two  of  the  gradients  multiplied  by  the  sine  of  the  angle  between  the 
last  gradient  and  the  plane  of  the  first  two.  This  final  product  of  sines 
relationship  extends  to  an  arbitrary  number  of  dimensions,  one  sine  term 
being  added  for  each  additional  dimension. 

For  the  5  parameter  case  noted  earlier  in  which  the  modified 
Newton's  method  converged  very  slowly,  this  measure  was  found  equal  to 
.00019,  a  value  indicating  that  the  5  gradients  are  very  nearly  co-planar 
and  thus  explaining  the  slow  convergence. 

Direct  Search,  by  contrast,  works  directly  with  the  sum  of  squares 
criterion,  computing  it  at  neighboring  points  to  the  current  point. 

Referring  to  Fig.  42,  Direct  Search  would  choose  a  direction  more  nearly 
in  line  with  the  gradients  and  would  thus  achieve  a  bigger  local  improve¬ 
ment  in  the  criterion.  It  does  not  extrapolate  a  long  distance  for  zero 
and  come  up  with  nothing  as  a  result. 
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APPENDIX  5 


GROUP  THEORETICAL  PROPERTIES  OF  RADIAL  WAVE  FUNCTIONS 
Lloyd  Armstrong,  Jr. 

Group  theory  has  been  applied  with  great  success  to  studies  of 
the  angular  portion  of  atomic  wave  functions.  It  is  hardly  an  exaggeration 
to  state  that  much,  If  not  most,  of  the  present  day  knowledge  of  the  actinide 
and  rare  earth  atoms  could  scarcely  have  been  obtained  without  use  of  this 
powerful  tool.  The  state-of-the-art  concerning  the  radial  portion  of 
atomic  wave  functions  is,  however,  quite  primitive  with  respect  to  that 
of  the  angular  portion  of  the  wave  function.  The  best  radial  wave  functions 
obtainable  result  from  some  type  of  Hartree-Fock  calculation.  This  type  of 
calculation  is  very  difficult  and  time  consuming  to  carry  out,  even  using 
the  best  of  modern  computers.  In  addition,  there  are  many  unsolved  theo¬ 
retical  problems  relating  to  the  Hartree-Fock  procedure. 

Rather  than  discuss  any  of  these  difficulties,  however,  let  me 
point  out  what  I  consider  to  be  the  greatest  weakness  of  such  a  procedure: 
the  resulting  wave  function  has  no  known  transformation  properties  with 
respect  to  any  symmetry  groups  of  the  atom.  Because  the  angular  wave 
functions  do  transform  in  known  ways  with  respect  to  certain  groups,  one 
is  able  to  immediately  write  down  a  number  of  selection  rules  on  angular 
matrix  elements  and,  even  more  important,  to  pick  out  types  of  matrix 
elements  that  must  differ  in  value  by  only  a  proportionality  factor.  This 
latter  result,  a  consequence  of  the  Wigner-Eckart  theorem,  enables  one  to 
calculate  the  values  for  one  set  of  matrix  elements  and,  after  finding 
a  single  constant  of  proportionality,  to  write  down  directly  the  values 
for  another  set  having  the  same  transformation  properties.  Unfortunately, 
since  no  transformation  properties  of  the  radial  wave  functions  are  known, 
such  a  labor-saving  simplification  cannot  be  used  to  relate  radial  matrix 
elements. 
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The  purpose  of  this  report  is  to  describe  ways  in  which  one  might 

hope  to  determine  the  transformation  properties  of  radial  wave  functions. 

We  wish  to  start  by  studying  the  simplest  possible  system,  hydrogen. 

52 

Pasternack  and  Stemhelmer  have  reported  a  very  interesting  selection 
rule  on  hydrogenic  radial  wave  functions: 

R  R 

(5.1)  ;  -gi-Hi!  dr  =  0  s  -  2,3. . . ,  | l-l'  |  +  1, 

S 

r 


where  the  complete  hydrogenic  wave  function  is  written 


n&m 


Rn*(r> 


Y 

SL  m 


(9,<f>)  . 


This  selection  rule,  which  has  successfully  defied  all  attempts  to  give 

it  a  simple  explanation,  looks  suspiciously  like  the  type  of  selection  rule 

that  results  from  a  group  property.  Thus  the  primary  goal  of  this  type 

of  research  must  be  to  explain  the  above  result  using  group  theory. 

Two  avenues  of  investigation  immediately  present  themselves. 

The  first  involves  searching  for  a  selection  rule  based  on  a  symmetry 

group  of  the  hydrogen  atom:  R(4),  the  proper  rotation  group  in  four 

dimensions.  The  second  involves,  a  study  of  the  groups  for  which  Laguerre 

polynomials  form  bases  of  representation.  This  latter  is  of  interest 

2i+l 

because  is  proportional  to  a  Laguerre  polynomial,  .  The  former 

may  perhaps  be  the  most  susceptable  to  study,  so  we  shall  consider  it 
first. 
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a)  The  Symmetry  Group  of  the  Hydrogen  Atom. 

The  R(4)  symmetry  group  of  the  hydrogen  atom  has  an  algebra 
consisting  of  the  three  components  of  the  angular  momentum  vector  ^  and 
the  three  components  Lenz  vector 

£ 

M  »(l/2y)(P  xL-LxP)-k- 
2 

where  k  ■  ze  .  Only  one  of  the  commutators  of  the  algebra  is  of  real 
Interest  at  the  moment 

[L2,M.  ]  3  -  2M  L ,  +  2M  L  +  2M  . 

i  Z  T  t  Z  I 
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% 


If  we  consider  the  operation  of  this  commutator  on  the  state  ^  nnf  we  find 


L  «+w  • <1+1)  <*«>  <v„u> 


which  implies  that 


Vnu  -  A*ul+M+1 


where  A  is  some  constant.  A  detailed  calculation  reveals  that 


(5.2) 


VnU 


Y£+1,&+1  (6,lj>)  £+1 

/ 2£+3  p  r 

Rnr 


We  shall  return  to  this  below. 

One  can,  in  fact,  find  a  simple  commutator  which  explains  Eq.  (5.1) 
if  one  uses  the  raising  operator  for  hydrogenic  radial  wave  functions  ^,55, 


(5.3a) 


R* ,  ±_  _  i±i  +  hL. 

K+  dr  r  £+1 


(5.3b) 


R~  R  =  B  R 
+  n£,  n£+l. 


We  consider  first  the  commutator 


*2T2  . 

[H  -  ^r,  CR]  -  0 
2mr2  q 


-&2  1  d2r  A2  zef 

whereH"  2^T  7  ^2  ^  "  r 


and  is  a  spherical  tensor  of  rank  k.  Multiplying  the  commutator  on 
the  left  by  (ll>  £J,  on  the  right  by |^n, ?  l 5 , ) ,  and  obtain 


rnH'£' 
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^En£  En'£'^  ^n£m  ^Cql  ^n'£'m'^ 

A2  ck 

m'ta  (Klm  ^  *n'1V)  [Z(l+1)  ~  ^'C^’+l)] 

For  n  =  n',  EnA  =  En,£l.  The  angular  portion  of  the  matrix  element  on 
the  right  above  can  easily  be  seen  to  be  non-vanishing  for  at  least  one 
value  of  k,  q,  leaving  the  result 


R  *  R  .  | 

6  rs  /  —2 - _B -  dr  (£(£+1)  _  £'(£•+!)) 


which  implied 
(5.4) 


r  Rn£  Rn£* 

J  - o -  dr  »  0 


if  £  +  V 


Consider  next  the  commutator 

fH  ^2L2  .kid  .  A 
[H - c  r  r]  -  0. 

2mr  *»  r  dr 

Proceeding  as  before,  and  letting  n  -  n' 


t(M-l) 

(5.5) 

- 

|ck± 
q  r 

d 

dr 

or 

k  1_  d_ 
r3  dl 


IMt+l)  -  t'U’+l)]  »nlm  l-fljrl 


+  U'W+"  1^1  -  0 


Again,  we  can  carry  out  the  angular  integration  which  is  the  same  for  both 

matrix  elements  and  which  has  at  least  one  non-vanishing  value.  We  then 
have 
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UU+l)  -  *'  U'+l)]  /  R  ,  dr 

T 

R  R 

+  2a'U'+l)  S  -n*  ■■  n£ -  dr  =  0 


Assume  that  £’>£,  and  use  Eqs.  (5.3a)  and  (5.3b)  to  treat  the  —  R  „ . 

dr  nl1 

term  above.  One  obtains 


i*(a+i)  -  t'a'+D]  /  t1—  r  Jl, 

r 


+  A  \v  *  B  R„,  dr 


-2S'(l'+l)  /  -7”-'  dr. 


Because  of  Eq.  (5.4),  the  above  simplifies  greatly,  leaving 


{[£(£+1)  -  £'(£'+l)]  (£’+l)  +  21'  (£'+l) } 

,  Rn£  Rn£.'  , 

/  - - -  dr  -  0 


or 


Rn£  Rni* 


(U'+l)  [i(t+l)  -  £'(£'-1)]  /  ■■  ■  3-  ■  dr  =»  0 

r 

which  implies 

;  nl  n*'  dr  -  0  £’>£+1. 

r 

All  the  remaining  selection  rules  can  be  proven  by  considering  commutators 
of  the  type 

<u_lt}L  ci<  i  £jl i.o 

lH  2m  2’  Cq  r  ,  n  °* 

r  H  dr 


and  proceeding  as  above. 

The  next  step  is  clearly  to  try  to  express  the  simple  commutators 
above  in  terms  of  the  generators  of  R(4).  In  particular,  if  we  consider 
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equations:  such  as  Eq.  (5.5)  with  m=£,  we  should  be  able  to  replace, 

in  some  fashion,  with  M+  (Eq.  (5.2)).  This  approach  provides  a  very 
encouraging  method  of  studying  the  vanishing  matrix  elements;  one  can 
hope  for  early  success  in  the  use  of  R(4)  symmetry  to  obtain  selection 
rules  on  the  radial  matrix  elements. 

b)  Symmetry  Groups  of  Laguerre  Polynomials. 

This  subject  has  been  studied,  in  some  detail  by  Miller.32 
We  shall  only  mention  here  some  of  the  major  results  of  his  work.  He 
considers  the  general  4  dimensional  complex  Lie  Algebra  defined  by 

[J+,J-]  =  2a2J3  -  bE;  [J3,J+]  »  J+; 

[J3,J_]  >  -  J“;  [J+,E]  *  [J~,E]= 

[J3,E]  =  0. 
for  all  complex  a,b. 

Clearly,  if  a=l,  b=0,  one  has  the  algebra  for  the  complex  three  dimensional 
rotation  group.  Miller  has  shown  that  the  basis  functions  for  a  realization 
in  two  variables  of  this  algebra  are  confluent  hypergeometric  functions. 

(The  Laguerre  polynomials  are  simply  a  special  case  of  the  confluent  hyper¬ 
geometric  functions.)  In  addition,  he  finds  that  when  b=l,  a=0,  one  also 
may  construct  basis  functions  from  confluent  hypergeometric  functions. 

The  algebra  in  this  case  is  the  complex  form  of  the  algebra  of  S^.  It 
would  appear,  at  this  time,  that  this  latter  group  is  of  little  value  and 
we  shall  confine  our  attention  to' the  case  a*l,  b*0. 

The  realization  in  two  variables  of  the  algebra  with  a*l,  b«0 
(G(1,0))  can  be  written 

J3  -  J*  -  (±f-  -  k(x)  ~  +  j  (x)) 

3y’  3x  3y 

E  =  u, 

where  u  is  a  complex  constant;  k(x)  and  j (x)  must  be  determined  so  as  to 
make  the  above  operators  satisfy  the  proper  commutation  relations.  Miller 
finds  two  possible  types  of  solutions  for  k(x)  and  j (x) : 
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Type  A  k(x)  *  cot  (x+p),  j (x) 


_ 9 _ 

sin (x+p) 


Type  B  k(x)  =  i,  j (x) 

58 


where  p,q  are  complex  constants.  As  an  example,  we  consider  Type  B 
representations  of  the  variety  labelled  D(u,m.).  We  first  change  co- 

_  4  y 

ordinate  systems  by  letting  t  ■  y  -  Itt / 2 ,  z  -  -  ie  .  Then 


9_ 

3t* 


XT 


3 

‘3z 


qz) 


If  the  solution  is  written  in  the  form  f  ■  Z  eraT ,  we  find 

m  m 

(^  +  m  -  qz)  Zn  -  -  <m  -  u) 

'  "  +  1*)  Z„  *  -  <■  +  “>  Zm-1 

which  has  solutions 

Zm  ■  (2qz)U+^  e  ^Z  (u-m+1 ;  2u+2;  2qz) 

and 

Z'  ■  (2qz)  u  e  -F.  (-u  -m;  -2u;  2qz) 
m  11 

where  ^F^  is  a  confluent  hypergeometric  function.  In  terms  of  Laguerre 
polynomials,  we  have 


_  .u+1  -qz  T 2u+l  . 

Z  -  (2qz)  e  L  .  (2qz) 
m  m-u-1 

Z;  -  <2«z>"u  e"qz  Cf1  <2<z>- 

This  type  of  representation  has  u  +  m  ^  integer. 

Of  greater  interest  to  us  are  the  basis  states  f  for  represen- 

m 

tation  labelled  +  ;  here 


f 

m 


T(2u)n!  T 2u-l  ,_x  ^u+n 
r (n+2u)  Ln  (z)  C 


where  u,n,  are  Integers.  Acting  on  this  state  with  raising  operators, 
one  finds 

T+  t  2u-l  _ 2u-l 
J  L  -*•  L  . 
n  n+1 

Basically,  this  result  implies  that  a  basis  for  the  representation  4- 
can  be  formed  by  the  collection  of  states  of  the  hydrogen  atom  obtained 
by  holding  £  fixed  and  letting  the  principal  quantum  number  run  over 
all  possible  values.  One  can  work  out  Clebsch-Gordan  Coefficients  for  the 
Kror.ecker  product  41  x  4-  ,  without  too  much  difficulty.^ 

One  therefore  has  the  transformation  properties  of  the  hydrogenic 
radial  wave  functions  with  respect  to  the  complex  3-dimensional  rotation 
group.  In  order  to  search  for  the  desired  selection  rules,  one  must 
determine  the  representation  to  which  the  operators  l/rS  belong.  This  is 
not  trivial  and  at  the  present  time  the  solution  is  unknown. 

This  latter  approach  would  appear  to  be  more  general  and  powerful 
than  that  of  (a) ,  but  the  problems  involved  are  correspondingly  more  dif¬ 
ficult.  None  the  less,  both  approaches  seem  at  the  present  to  be  rather 
hopeful. 
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Most  contemporary  studies  of  atomic  structure  depend 
heavily  on  the  use  of  group  theory.  In  the  central  field  model 
of  the  atom,  wavefunctions  can  be  written  as  a  product  of  an 
angular  function  and  a  radial  function.  The  radial  function  is 
usually  calculated  through  use  of  some  variational  technique 
(e.  g. ,  Hartree-Fock)  or  simply  treated  in  terms  of  a  variable 
parameter;  the  angular  portion  is  expressed  in  terms  of  pro¬ 
ducts  of  spherical  harmonics.  It  is  to  this  angular  portion  of 

59 

the  wavefunction  that  group  theory  has  traditionally  been  applied. 

Quite  impressive  simplifications  have  resulted  from  this 
application  of  groups  to  the  angular  wavefunction.  For  example, 
selection  rules  based  on  Kronecker  products  permit  one  to 
identify  many  matrix  elements  as  being  identically  zero  without 
explicit  (and  often  laborious)  calculation.  Use  of  the  Wigner- 
Eckart  theorem  often  enables  one  to  ascertain  sets  of  angular 
matrix  elements  which  are  proportional  to  one  another,  thus  saving 
a  great  deal  of  computational  effort. 

The  use  of  group  theory  has  been  invaluable  in  advancing 
the  understanding  of  angular  wavefunctions  and  atomic  shell 
structure  to  its  present  state.  Unfortunately,  although  very  im¬ 
pressive  techniques  have  been  developed  in  the  last  few  years 
for  the  calculation  of  radial  wavefunctions,  very  little  has  been 
done  on  the  study  of  the  symmetry  properties  of  radial  wave- 
functions.  Most  of  the  group  theoretical  work  which  involves 
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radial  wavefunctions  has  consisted  of  studies  of  hydrogen  using 

0(4),  0(4,1),  etc.  In  this  case,  however,  the  radial  function 

has  not  been  studied  directly,  since  the  basis  functions  for  these 

groups  are  composed  of  a  product  of  radial  and  angular  parts. 

We  wish  to  report  on  a  study  of  the  group  properties  of  the 

hydrogenic  radial  functions  themselves.  Two  previous  results 

suggested  strongly  that  these  functions  should  be  amenable  to 

a  group -theoretic  study.  First,  a  well-known  result  of  Paster - 

52 

nack  and  Sternheimer  shows  that,  for  hydrogenic  radial  wave- 
functions , 


J  RntRnt' 
0 


dr 


0 


t>  V 

s  s  2,  3, .  1 


(we  choose  our  radial  wave  function  to  be  R  n^(r)/r).  This  result 
is  quite  suggestive  of  a  group -theoretical  selection  rule. 

Second,  the  expectation  values  of  rn  (n  both  positive  and  nega¬ 
tive)  evaluated  with  hydrogenic  functions6® can  be  seen  to  be 
proportional  to  Clebsch-Gordan  coefficients61  in  which  i>  plays 
the  part  of  the  angular  momentum,  and  n,  the  projection  of  the 
angular  momentum. 

The  group  algebra  pertinent  to  this  study  is  formed  by  the 
operators 


J 

J 


3 

± 


(1) 
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which  satisfy  the  commutation  relations 


[J+.J_]  =  2J3 


(2) 


(Operators  of  this  same  general  type  have  been  considered  by 
Miller.  )  We  wish  to  construct  functions  f^  such  that 


J3  {ln  =  n  fln 


J±‘tn 


[Hl+D  -  n(n±l)]^n±1  . 


(3) 


A  function  satisfying  Eqs.  (3)  is  given  by 

,  ((i)n-M  e-,«  ‘(Y, 

n  2ir  (n  +  f )  !  J 

a  62 

where  L  ^  is  the  Laguerre  polynomial  of  Morse  and  Feshback. 

For  fixed  n,  this  function  is,  of  course,  proportional  to  the  radial 

62 

wavefunction  for  hydrogen.  We  define  a  Hilbert  space  as  the 
space  of  functions  f^  with  inner  product 

<  f<t,Jf<,«n'  >  =  Jf4nfi«n,d° 

-2 

where  df)  =  x  dTdz;  the  parameter  t  can  range  from  0  to  2tt, 
the  parameter  z  from  0  to  ®  .  One  can  easily  show  that 


In  addition,  one  has 
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<  ftnlJ3lft'n' >  =  <J3£<.nlft'n'> 


Consideration  of  the  operators  shows  that 


J -*1,1+1 


0  . 


providing  a  lower  bound  to  the  set  of  functions  f  ;  clearly, 
there  is  no  upper  bound.  The  representation  thus  formed  i.B 
infinite  dimensonal.  It  can  also  be  Bhown  to  be  irreducible, 
implying  that  the  group  described  by  the  algebra  (1)  and  (2) 
is  non-compact. 

Ik) 

We  define  a  spherical  tensor  T'  '  by  the  relationships 


tJ±’T<k’ql  *  (k(k+l)-q(q±l)]^  T(k)q±  , 

[Jj.T(k,q)  -  qT(k,q  .  (5) 


This  definition  assures  that  the  operators  T'  ^  transform  in 

the  same  way  as  a  state  f.  .  Note,  however,  that  the  operators 

Kq 

T^  (or  the  states  f,  )  with  |q|  s  k  form  a  finite  dimensional 
q  Kq 

representation  of  the  algebra.  We  shall  be  concerned  here  only 
with  operators  of  this  type.  Because  of  the  transformation 
properties  of  the  T^  we  can  use  the  Wigner-Eckart  theorem 
to  show  that 


<WI'°Vt'n'>  *  A(kq.t'n'|/n)  U||T(k)|U')  (6) 
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where  (/  ||T^||/  ')  is  a  reduced  matrix  element  which  is  inde¬ 
pendent  of  q,  n',  and  n;  and  A(kq;  t'n'I^n)  is  a  coefficient  which 

Ik) 

depends  on  the  transformation  properties  of  f^,  f^,  rI  ,  and  Tv  ' 

but  not  on  their  explicit  form.  One  can  determine  the  coefficients 
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A  following  the  technique  of  Racah.  By  considering  the  matrix 
elements  of  [J±  ,  between  the  states  <  f^n|  and  |f^,n,>  , 

and  the  coupling  of  states  f^,^  and  to  a  state  f^,  one 

can  easily  obtain  recursion  relations  which  demonstrate  that  the 
dependence  of  A(kq,  fn'  |-tn)  on  k,  q,  V,  n',  L  and  n  is  identical 
to  the  dependence  of  the  algebraic  form  of  the  Clebsch-Gordan 
coefficient  (kq,  i'n'  |{,n)  on  the  same  quantum  numbers.  The  A 
coefficient  is  not  strictly  identical  to  the  corresponding  Clebsch- 
Gordan  coefficient,  however,  since  in  the  present  case  n  >  t, 
n'  >  <?',  values  for  which  the  Clebsch-Gordan  coefficient  vanishes 
identically.  Nonetheless,  one  can  obtain  an  algebraic  expression 
for  A(kq,  l  'n'  |{,n)  for  a  fixed  k  and  q,  and  general  n', -C, ,  n 
by  determining  the  algebraic  expression  for  the  equivalent  Clebsch- 
Gordan  coefficient  and  relaxing  the  restriction  that  | n*  |  <  -f  •  , 

|  n  |  <  l . 

Changing  for  convenience  from  Clebsch-Gordan  to  3-j 
symbols,  we  can  write 

<ftolT"0,IW>  *  ('  k  \)  «»T“V>  <7) 

\  -n  q  n*  / 
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where  the  subscripted  3-j  symbol  refers  to  the  algebraic  form 
of  the  usual  3-j  symbol  with  the  restrictions  on  the  range  of 

n  and  n'  relaxed.  Algebraic  forms  of  certain  3-j  symbols 

61  ,,  64 

have  been  given  by  c.  g.  ,  Edmonds  and  Falkoff  ctjal.  . 

We  now  consider  the  tensorial  ranks  of  the  operators 

eimT/rN  (Ns  0).  By  using  Eqs.  (5),  one  can  easily  show  that 

such  an  operator  corresponds  to  the  mth  component  of  a  tensor 

of  rank  N,  .  By  relating  the  function  (4)  to  the  radial  wave- 

function  of  hydrogen,  we  can  therefore  make  the  identification  (N>  1) 


dr 


1  ,1*?  (-!)*""( ♦i)*'" 

[(2<.+  l)(2*'+l)J 


FT 


(~> 


if1 

1*  \  -n 


N-2 


(4||t(N'2)|U‘) 


One  notes  immediately  that  the  selection  rules  of  Pasternack 
52 

and  Sternheimer  are  contained  in  the  triangular  selection  rules 
of  the  3-j  symbol  contained  in  Eq.  (8). 
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_N 

The  proportionality  of  expectation  values  of  r  to  3-j  symbols 
is  also  made  obvious  by  Eq.  (8). 


In  order  to  fully  utilize  Eq.  (8)  one  must  be  able 
to  evaluate  the  pertinent  reduced  matrix  elements.  This  is 
easily  done:  assume  that  £> /’,  and  consider  matrix  elements 
in  which  the  bra  is  the  state  <f  .  |.  Then,  for  example 

VM'  '  ^ 


<  f 


/.  l+l 


t+1 


>  = 


1  )  Kl|T(k||k') 

0  1*1 1 


-z  L+V  -k 
z 


Li.i<  <*> dz 


.m 


The  integral  above  can  be  simply  evaluated,  as  can  the  3-j 
symbol.  One  then  obtains 


ul|T(k)|U«)  =  vlmV 


[ 


(2-t+l )(2^»+l )k!  2(m'-k)!  A 
(k+^+^'V  1 )!  (I  -j'+k)!  {l+k-l*)'.  J 


(10) 
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♦ 


Use  of  these  expressions  in  conjunction  with  algebraic  forms 

of  3-j  symbols  leads  immediately  to  the  known  expectation 
-N 

values  r  (N  positive). 

This  successful  application  of  the  techniques  o.  group 
theory  to  a  study  of  radial  wavefunctions  suggests  strongly 
that  such  studies  of  more  complicated  radial  functions  may  be 
possible.  In  addition,  it  is  particularly  gratifying  to  demon¬ 
strate  that  the  radial  and  angular  parts  of  the  atomic  wave- 
function  can  be  placed  on  a  more  even  footing  and  treated 
theoretically  using  the  same  general  techniques. 


The  author  gratefully  acknowledges  invaluable  discussions 
with  Professors  B.  R.  Judd  and  G.  Domokos. 
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TABLE  I.  The  Measured  Energy  Levels  ^  and  SLJ  Labels  for  the  Lowest 
Five  Configurations  of  Xe  I. 

5  2 

The  core  configuration  is  always  5p  (  P)  and  therefore  is 
not  listed. 


n'l' 

SLJ 

E(cm“1) 

n'l' 

SLJ 

E(cm  ^ 

5p 

lso 

0 

5d 

3po 

79,772 

6s 

3p2 

67,068 

3°i 

79,987 

3pl 

68,046 

3f4 

80,197 

3po 

76,197 

3f2 

80,323 

lpi 

77,186 

3f3 

80,971 

6p 

3Di 

77,270 

3D2 

81,926 

3d2 

78,120 

3d3 

82,431 

3d3 

78,404 

3pl 

83,890 

3pl 

78,957 

3p2 

91,153 

3p2 

79,213 

1d2 

91,448 

3po 

80,119 

lp3 

91,747 

3S1 

88,380 

lpl 

93,619 

!d2 

89,163 

7s 

3p2 

85,189 

lpl 

89,279 

3pl 

85,441 

lso 

89,861 

3po 

95,721 

lpl 

95,801 
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TABLE  II.  The  Measured  Energy  Levels  and  SLJ  Labels  for  Fourteen 
Configurations  of  Kr  I. 

5  2 

The  core  configuration  is  always  4p  (  P)  and  therefore  is 
not  listed.  The  values  in  parentheses  are  approximate 
values. 


I'l' 

SLJ 

E(cm  ^  ) 

n'l' 

SLJ 

E(cm 

4p 

\ 

0 

6s 

3p2 

99,628 

5s 

3p2 

79,973 

99,895 

3P1 

80,918 

3po 

105,092 

3po 

85,192 

lpi 

105,147 

S 

85,848 

6p 

3p! 

102,888 

5p 

3°1 

91,169 

3p>3 

103,116 

S 

92,295 

3°2 

103,122 

3°2 

92,308 

3pi 

103,314 

3P1 

92,965 

3p2 

103,363 

3p2 

93,124 

3po 

103,762 

3po 

94,094 

3p! 

108,439 

3S1 

97,597 

lp! 

108,515 

1P1 

97,920 

1d2 

108,569 

\ 

97,946 

\ 

108,822 

lso 

98,856 

5d 

3®1 

103,803 

4d 

3po 

96,772 

3po 

104,074 

\ 

97,086 

3p4 

104,631 

3p2 

97,689 

3p3 

104,917 

3p4 

97,798 

3p2 

105,008 

3p3 

98,227 

3°2 

105,164 
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TABLE  II 


(Continued) 


n’l* 

SLJ 

E(cm 

n'l' 

SLJ 

E(cm  ■*■) 

Ad 

S 

98,868 

5d 

S 

105,209 

3d3 

99,080 

\ 

105, 6A9 

99.6A7 

S 

110, 10A 

lD2 

103,267 

3p2 

110,123 

103, AA3 

lp3 

110,238 

lp3 

103,702 

lf>l 

110.73A 

\ 

10A.888 

7s 

3p2 

105, 6A8 

9s 

3p2 

109,753 

3pi 

105,772 

^1 

109,780 

3po 

(111,00A) 

3po 

(115,128) 

lpi 

lll.OOA 

1P2 

115,128 

Af 

\ 

105,965 

10s 

3p2 

110,609 

3p2 

105,966 

3P! 

110,620 

3(35 

105,990 

3po 

(115,961) 

\ 

(105,990) 

lpi 

115,961 

3(I3 

106,022 

11s 

3p2 

111,155 

3°2 

106,022 

3pi 

111,172 

3F3 

106, 0A8 

3po 

(116,527) 

3fA 

106, 0A8 

\ 

116,527 

1(5a 

111,379 

12s 

3p2 

111,529 

3°3 

111,379 

3pi 

111,537 
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Table  II.  (Continued) 


1 1* 

SLJ 

E(cm  "*■) 

n'l’ 

SLJ 

E(cm  ) 

4f 

lp3 

111,381 

12s 

3p0 

(116,904) 

1d2 

111,382 

lpl 

116,904 

8s 

3p2 

108,326 

\ 

108,374 

3po 

(113,711) 

lp, 

113,711 
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TABLE  III.  The  Intermediate  Coupling  Wave  Functions  for  the  Five  Lowest  Configurations 
of  Xe  I. 
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TABLE  III  (Cont) 
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TABLE  IV.  The  intermediate  coupling  wave  functions  for  fourteen 

configurations  of  Krl.  The  core  configuration  is  always 
4p5(2P) . 


1*1* 

SL 

Eigenvectors 

5p 

3P 

0.738 

-0.675 

xs 

0.675 

0.738 

6p 

3P 

0.786 

-0.618 

xs 

0.618 

0.786 

5s 

3P 

0.717 

-0.697 

1P 

0.697 

0.717 

5p 

3d 

0.192 

0.908 

0.370 

0.027 

3P 

0.725 

0.059 

-0.487 

-0.484 

3s 

0.482 

-0.068 

-0.147 

0.861 

3P 

0.453 

-0.408 

0.777 

-0.153 

4d 

3d 

0.144 

0.961 

0.237 

3P 

-0.712 

-0.065 

0.699 

XP 

0.687 

-0.269 

0.675 

6s 

3P 

0.798 

-0.603 

XP 

0.603 

0.798 

6p 

3d 

0.225 

0.843 

0.486 

0.045 

3P 

-0.767 

-0.091 

0.474 

0.424 

3S 

0.474 

-0.151 

-0.038 

0.867 

1p 

0.370 

-0.508 

0.733 

-0.259 

5d 

3d 

0.173 

0.945 

0.277 

3P 

0.766 

0.048 

-0.641 

Ill 


TABLE  IV.  (Continued) 


n'l' 

SL 

Eigenvectors 

XP 

0.620 

-0.323 

0.716 

7s 

3P 

0.809 

-0.588 

XP 

0.588 

0.809 

8s 

3P 

0.814 

-0.581 

*P 

0.581 

0.814 

9s 

3P 

0.815 

-0.580 

XP 

0.580 

0.815 

10s 

3P 

0.816 

-0.578 

XP 

0.578 

0.816 

11s 

3P 

0.815 

-0.579 

*P 

0.579 

0.815 

12s 

3P 

0.816 

-0.578 

XP 

0.578 

0.816 

5p 

3D 

0.653 

-0.208 

-0.728 

3P 

0.477 

0.360 

0.181 

XD 

0.589 

-0.466 

0.661 

4d 

3f 

0.357 

0.844 

0.401 

-0.015 

3d 

-0.615 

0.024 

0.518 

0.594 

3P 

0.373 

0.081 

-0.472 

0.795 

1D 

0.596 

-0.530 

0.590 

0.124 

6p 

3d 

0.653 

-0.265 

-0.709 

3P 

0.483 

0.867 

0.121 
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TABLE  IV.  (Continued) 


Eigenvectors 


1 

D 

0.583 

-0.422 

0.694 

3f 

0.372 

0.795 

0.479 

0.009 

3d 

-0.658 

-0.014 

0.524 

0.541 

3P 

0.447 

-0.034 

-0.307 

0.840 

XD 

0.479 

-0.605 

0.634 

-0.047 

3f 

0.617 

0.786 

0.021 

3d 

-0.539 

0.404 

0.739 

XD 

0.573 

-0.467 

0.674 

3f 

-0.493 

0.292 

0.819 

3d 

0.610 

0.788 

0.086 

XF 

0.620 

-0.542 

0.567 

3f 

-0.620 

0.200 

0.759 

3d 

0.520 

0.829 

0.206 

XF 

0.588 

-0.522 

0.617 

3g 

0.485 

0.674 

0.558 

-0.001 

3f 

-0.652 

0.001 

0.566 

0.505 

3d 

0.336 

0.073 

-0.380 

0.859 

XF 

-0.477 

0.735 

-0.474 

-0.086 

3g 

-0.750 

-0.011 

0.661 

3f 

0.327 

0.863 

0.386 

:lg 

0.575 

-0.505 

0.643 

113 


TABLE 


Atom 


Xe 


Kr 


.  The  values  of  Slater  integrals  and  spin-orbit  coupling  constants  which  was 
determined  by  fitting  the  theoretical  energy  levels  to  the  measured  levels. 
The  parameters  and  the  rms  deviation,  o,  of  the  fit  are  given  in  cm 


n'l' 

Fo 

F2 

Go 

G1 

G2 

G3 

G4 

6nl 

6n’l’ 

6 

6s 

71040 

951 

6239 

233 

6p 

81866 

173 

629 

33 

7005 

382 

105 

5d 

84161 

186 

298 

28 

6892 

-175 

191 

7s 

88875 

175 

700  L 

32 

5s 

82512 

800 

3482 

4 

5p 

94118 

176 

650 

42 

3538 

153 

46 

4d 

99747 

175 

177 

14 

3501 

-  35 

138 

6s 

101626 

176 

3608 

54 

6p 

104922 

45 

204 

13 

3556 

60 

17 

5d 

106520 

105 

88 

6 

3539 

11 

245 

7s 

107507 

73 

3545 

40 

4f 

107800 

4 

0 

0 

3580 

0 

1 

8s 

110150 

29 

3580 

16 

9s 

111561 

16 

3578 

8 

10s 

112399 

6 

3566 

3 

11s 

112956 

10 

?578 

5 

12s 

113326 

5 

3582 

3 
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TABLE 


Atom 


Xe 


VI.  The  transition  probability  rates  A  for  emission  between 

mn 

the  energy  levels  of  the  five  lowest  configurations  of  Xe  I 

and  between  the  energy  levels  of  fourteen  configurations  of  Kr  I. 

The  core  configuration  is  always  5p  (  P)  and  4p  (  P)  for 

xenon  and  krypton,  respectively.  Only  values  of  A  larger  than 
8  X 

0.1  x  10  sec  was  listed. 

State  n  State  m  X  A 

o  Q®®  —  1 

n*!* SLJ  n'T1  SLJ  A  10  sec 


ls0 

7s 

lpi 

1044 

0.9 

5d 

lpi 

1068 

27.8 

7s 

3pi 

1170 

1.1 

5d 

3pi 

1192 

13.0 

3°1 

1250 

1.1 

6s 

lpi 

1296 

5.2 

3pi 

1470 

3.9 

3po 

6p 

1P1 

7644 

0.3 

lpl 

lso 

7890 

0.4 

3po 

3S1 

8209 

0.1 

3p2 

3p2 

8234 

0.3 

lpi 

lpi 

8269 

0.1 

3pi 

3po 

8282 

0.4 

lpi 

lD2 

8349 

0.4 

3p2 

3°3 

8822 

0.3 

lpi 

3S1 

8933 

0.2 

3pi 

3p2 

8955 

0.1 

3p2 

3°2 

9048 

0.1 
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TABLE  VI.  (Continued) 


Atom 


Kr 


State 

.’l* 

n 

SLJ 

State 

n"l" 

m 

SLJ 

X 

0 

A 

A 

gtnn 

10  sec 

3pi 

3pi 

9165 

0.3 

S 

3°1 

9802 

0.2 

3pi 

S 

9926 

0.2 

4p 

lso 

6s 

lpi 

951 

0.8 

4d 

lpi 

953 

11.9 

6s 

3pi 

1001 

1.1 

4d 

3pi 

1004 

9.9 

\ 

1030 

1.2 

5s 

1P1 

1165 

4.7 

3pi 

1236 

3.9 

5s 

3pi 

5p 

3po 

7589 

0.5 

3p2 

3p2 

7604 

0.3 

1P1 

lso 

7687 

0.4 

3po 

lpi 

7857 

0.2 

5p 

\ 

5d 

\ 

7916 

0.1 

5s 

3po 

5p 

3§1 

8062 

0.2 

3p2 

3°2 

8107 

0.1 

\ 

8115 

0.4 

3pi 

3p2 

8192 

0.1 

lpi 

S 

8266 

0.4 

lpi 

8283 

0.2 
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TABLE  VI.  (Continued) 


State  n 

State  m 

A 

A 

Atom 

n'l*  SLJ 

n"l"  SLJ 

0 

A 

in8mn  -1 
10  sec 

3pi 

8300 

0.3 

1P1 

3si 

8511 

0.2 

3pi 

8779 

0.2 

3p2 

8931 

0.2 
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Table  VII.  Comparisons  of  the  four  lamps  pumping  a  0.7%  NdrYAG  slab. 

P.  is  the  electrical  power  density  into  the  lamp  and  P 
in  out 

is  the  power  density  out  of  the  laser  slab. 
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TABLE  VIII.  A  comparison  of  accuracy  and  speed  of  different  integration 
methods  for  the  evaluation  of  orthonormality  and  Slater 
integrals  for  the  Cu  4s  electron. 

Two  Laguerre  integrations  with  n  =  15  (Lj^)  and  n  =  48  (L^g) 
are  compared  with  Simpson's  rule  of  integration.  The 
per  cent  errors  in  the  Gk  values  for  the  L15  integration 
were  invalidated  by  a  programming  error  and  are  not  given. 
However,  the  quoted  execution  times  are  for  the  evaluation 
of  all  integrations.  Only  errors  larger  than  0.1%  are  listed. 


n'l' 

Integration 

Method 

Absolute 
Error  in 
Orthogonality 
Integral 

] 

F° 

Per  Cent 
G° 

Errors 

G1  G2 

Execution 

Time 

In 

Seconds 

4s* 

l15 

20.3 

l48 

0.1 

S 

0.1 

3d 

l15 

1.4 

0.7 

l48 

1.6 

6.2 

S 

0.6 

32.8 

3P 

l15 

0.5 

0.6 

l48 

2.0 

5.7 

S 

0.5 

30.3 

-5 

3s 

L15 

2  x  10 

0.6 

0.5 

l48 

<10"6 

1.0 

5.1 

S 

<10-6 

0.2 

27.0 

2s 

l15 

1  x  10"3 

2.9 

0.5 

l48 

<10"6 

3.1 

4.9 

S 

7  x  10"5 

0.1 

27.0 

Is 

l15 

4  x  10-3 

12.4 

1.3 

00 

<10-6 

8.8 

4.9 

s 

<10-6 

0.3 

26.7 

*  These  execution  times  are  included  in  Table  IX. 
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TABLE  IX.  A  comparison  of  accuracy  and  speed  of  different  integration 
methods  for  the  evaluation  of  <rn>  radial  integrals  for  the 
Cu  4s  electron.  Two  Laguerre  integrations  with  n  =  15  (L15) 
and  n  =  48  (L^g)  are  compared  to  the  trapezoidal  rule  (T) 
and  Simpson's  rule  integration  (S).  The  accuracy  is  given 
in  per  cent.  Errors  less  than  0.1%  are  not  listed. 


Integration 

Method  -2 

-1 

n  in  <rn> 

0  1 

2  4  6 

Execution 

Time 

in 

Seconds* 

L15 

10.1 

2.0 

0.1 

0.7 

r 

is 

00 

3.3 

T 

6.4 

0.1 

0.1  0.1 

0.1  0.1  0.1 

13.9 

S 

4.5 

15.0 

*  Includes  the  calculation  of  F° . 


120 


TABLE  X.  The  screening  constants  proposed  by  Burns 


i 

iY 

n  -  1 

n 

n  + 

s 

s 

0.90 

0.40 

0.10 

P 

0.75 

0.35 

0.10 

d 

0.50 

0.35 

0.10 

f 

0.  AO 

0.20 

p 

s 

1.00 

0.50 

P 

0.90 

0.35 

d 

0.70 

0.35 

f 

0.30 

0.20 

d 

s 

1.00 

0.75 

0.15 

P 

1.00 

0.50 

0.05 

d 

0.90 

0.35 

0.05 

f 

1.00 

0.20 

f 

s 

1.00 

1.00 

0.50 

P 

1.00 

0.75 

0.15 

d 

0.85 

0.50 

f 

1.00 

0.35 

» 


TABLE  XI.  Average  screening  constants  obtained  from  Mann's  binding 
energies . 


5, 

£’\n' 

n  -  2 

n  -  1 

n 

n  +  1 

n  +  2 

n  +  3 

n  +  4 

n  +  5 

s 

s 

1.000 

0.875 

0.644 

0.133 

0.157 

0.141 

0.103 

0.062 

P 

1.000 

0.776 

0.505 

0.286 

0.200 

0.130 

0.070 

0.028 

d 

1.000 

0.918 

0.560 

0.321 

0.183 

0.115 

0.062 

0.025 

f 

0.982 

0.884 

0.610 

0.354 

0.186 

0.072 

0.029 

p 

s 

1.000 

1.000 

0.397 

0.079 

0.150 

0.162 

0.135 

0.085 

p 

1.000 

1.000 

0.706 

0.308 

0.237 

0.168 

0.105 

d 

1.000 

0.936 

0.576 

0.322 

0.220 

0.150 

0.095 

f 

1.000 

0.918 

0.628 

0.358 

0.187 

0.106 

d 

s 

1.000 

1.000 

1.000 

-0.134 

0.038 

0.148 

0.155 

P 

1.000 

1.000 

1.000 

0.231 

0.235 

0.198 

d 

1.000 

1.000 

0.747 

0.313 

0.237 

0.175 

f 

1.000 

0.996 

0.661 

0.362 

0.224 

f 

s 

1.000 

1.000 

1.000 

1.000 

0.020 

0.020 

P 

1.000 

1.000 

1.000 

0.801 

0.113 

d 

1.000 

1.000 

1.000 

0.218 

0.120 

f 

1.000 

1.000 

0.866 

0.352 

n  +  6 


0.025 
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TABLE  XII.  Hartree-Fock  values  of  radial  integrals  involving  the  Cu  4s 
wave  function. 38 ,48 


Slater  Integrals  between  4s 
nY  F° 

and  n'S,' 
k 

Wave  Functions 

Gk 

4s 

n 

Radial  Integrals 
<rn> 

4s 

.5450896 

-2 

.8077270 

3d 

.7122117 

2 

.0636135 

>1 

.3822107 

3p 

.7235111 

1 

.0232793 

0 

1.000000 

3s 

.7248400 

0 

.0164793 

1 

3.331133 

2 

13.08433 

4 

310.3408 

6 

12034.72 

123 


TABLE  XIII.  Initial  and  final  Slater  parameters  for  four  different 

three-term  expansions.  The  final  parameters  are  the  values 
after  100  seconds  of  computing  time.  The  standard  deviation 
is  between  the  five  true  radial  Integrals  and  the  five 
calculated  radial  integrals. 


Expansion 

Set 

am 

^m 

Standard 

Deviation 

Description 

No. 

m 

n£ 

Initial 

Final 

Initial 

Final 

Initial 

Final 

Switched  Input 

1 

1 

4s 

0.53 

0.26 

0.70 

2.15 

2 

3d 

0.50 

0.84 

-0.50 

-1.49 

3 

3p 

1.25 

1.03 

-0.31 

8.983 

0.069 

Optimized  Set 

2 

1 

4s 

0.22 

0.23 

-0.64 

-0.71 

2 

3d 

0.21 

0.35 

2.27 

0.35 

3 

3p 

0.23 

-0.79 

1.21 

0.871 

0.052 

Burns 

3 

1 

4s 

0.53 

0.16 

0.70 

1.11 

2 

3d 

1.03 

3.46 

-0.50 

0.23 

3 

3p 

1.25 

0.50 

-0.59 

0.559 

0.041 

Mann 

4 

1 

4s 

0.20 

0.11 

0.95 

-1.38 

2 

3d 

0.29 

0.28 

0.23 

2.21 

3 

3p 

0.72 

-0.23 

0.14 

0.187 

0.030 
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TABLE  XIV.  Initial  and  final  Slatar  parameter*  for  three  different 

five-tern  expansione.  The  final  paranetera  for  the  flrat 
three  aeta  are  the  rcault  of  100  aeconda  of  computing  time. 
Set  no.  8  la  the  continuation  of  aet  no.  7.  The  atandard 
deviation  la  between  the  five  true  radial  integrals  and  the 
five  calculated  radial  integrals. 


Expansion 


Description 

Set 

No. 

m 

nt 

*m 

- tm. 

Initial 

Final 

Standard 

Initial 

Deviation 

Final 

Mann 

5 

1 

3d 

0.29 

0.90 

0.70 

2 

4s 

0.20 

0.43 

-0.44 

3 

3a 

0.91 

-0.09 

-0.10 

A 

2a 

2.60 

0.03 

0.03 

5 

la 

7.39 

-0.01 

-0.01 

0.554 

0.032 

Sjj  Peaked 

6 

1 

3d 

0.10 

0.23 

-0.41 

2 

4a 

0.11 

0.94 

-0.60 

3 

3a 

0.47 

-0.23 

-0.02 

4 

2a 

1.43 

0.07 

0.07 

5 

la 

8.33 

-0.01 

-0.01 

1.321 

0.030 

pm  Peaked 

7 

1 

3a 

0.47 

0.23 

0.42 

2 

4p 

0.23 

0.94 

0.97 

3 

3p 

0.93 

-0.23 

-0.67 

4 

2p 

2.86 

0.07 

0.06 

5 

la 

8.33 

-0.01 

-0.03 

0.130 

0.010 

Final  Set 

8 

1 

3a 

0.47 

0.42 

-0.01 

2 

4p 

0.23 

0.97 

1.06 

3 

3p 

0.93 

-0.67 

-0.31 

4 

2p 

2.86 

0.06 

0.10 

5 

la 

8.33 

-0.03 

-0.02 

0.010 

0.0002 
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TABLE  XV.  The  determination  of  the  Initial  set  of  Slater  parameters 
for  expansion  sets  5,  6,  and  7. 


m 

n£ 

xmax 

Set  S 

Set  6 

Set  7 

am 

am 

P  m 

exact 

P  m 

rounded 

n'i' 

am 

1 

3d 

2.1A 

0.29 

1  .A0 

0.61 

1 

3s 

0.A7 

2 

As 

8.86 

0.20 

0.11 

1.76 

2 

Ap 

0.23 

3 

3s 

2.1A 

0.91 

0.A7 

1.95 

2 

3p 

0.93 

A 

2s 

0.70 

2.60 

1.A3 

1.82 

2 

2p 

2.86 

5 

Is 

0.12 

7.39 

8.33 

0.89 

1 

Is 

8.33 
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TABLE  XVI.  Accuracy  of  the  PF  wave  functions  for  the  eight  final  PF 
expansions  given  in  Tables  XIII  and  XIV. 


Expansion 

Set 

No. 

Standard 

Deviation 

of 

Integrals 

Percentage 
Standard 
Deviation 
of  Wave 
Function 

Percentage  Errors 

of  <rn> 

tM 

1 

a 

c 

-1 

0 

1 

2 

A 

6 

1 

0.069 

69.1 

7A.6 

20. A 

-A. A 

-16.7 

-31.2 

-59.2 

-78.3 

2 

0.052 

35.2 

-53.1 

O* 

• 

CM 

1 

0.9 

5.8 

8.1 

5.1 

-A.  1 

3 

0.0A1 

32.2 

-53.  A 

12.6 

1.2 

-5.9 

-8.2 

-1.2 

18.1 

A 

0.030 

8A.7 

-76.6 

-13.7 

1.6 

11. A 

25.8 

1A0.3 

57A.1 

5 

0.032 

A6.A 

-80.8 

-17.7 

2. A 

1A.7 

22.2 

22.7 

7.7 

6 

0.030 

50.5 

-17.0 

2.6 

2.5 

18.3 

52.3 

208.9 

57A.7 

7 

0.010 

13.1 

-2A.8 

3.6 

-1.5 

-6.0 

-10.5 

-19.1 

-27.5 

8 

0.0002 

8.2 

-A8.1 

-3.2 

-0.0 

0.3 

-0.8 

-6.0 

-1A.0 
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Fig.  5-Average  electron-neutral  momentum-transfer  cross  sections  for  krypton 
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Fig.  7  -^-Factor  for  xenon  —  SchlUter-, - Biberman  et  al> 

— — Yankov, . Liebermann;  o  measurements  by  Berge  et  at 
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Fig.  8— ^-factor  for  krypton  — Schluter-,  -  -  Biberman  et  al; 

Ameasurements  by  Carls 


Curve  592300-8 


[ 

r 


Fig.  9-Percentage  of  the  total  radiation  from  a  2.5  atmospheres  arc  discharge 
which  is  continuum  radiation.  The  wave  length  interval  is  from  2000 

to  20,000  angstroms 
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Curve  592298- B 


Fig.  10-Total  radiation  from  a  2. 5  atmospheres  arc  discharge  in  the  wave 
length  interval  from  2000  to  20,000  angstroms 


Curve  592296-A 


Fig.  11— The  electrical  characteristics  of  a  2.5  atmospheres  arc 
discharge.  The  central  core  temperatures  are  given  in  103  °K.  The 
points  on  the  curves  indicate  the  temperatures  at  which  the 
characteristics  were  calculated 
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jre  profiles  of  a  40  amp-2.5  atm  arc  discharge. 
□  and  •  denote  calculated  values 


Fig.  13-Spectrai  radiance  at  the  surface  of  a  40  amp-2. 5  atm 
xenon  arc  discharge 
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Fig.  14-Spectral  radiance  of  the  8Z32A  line  at  the  surface  of  a 
40  amp-2. 5  atm  xenon  arc  discharge 
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Fig.  16-Radiant  emittance  of  an  emission  line  at  9200°K  in  a 
§  40  amp-2. 5  atm  xenon  arc  discharge 


WAVELENGTH  (HNGSTRGV5/ 


Fig.  17-Measured  spectral  radiance  of  a  40  amp-2.5  atm  xenon  arc  discharge. 

p 

Full  scale  is  4.96  x  10  watts/(cm  ster  A) 


Fig.  18~Spectrai  radiance  at  the  surface  of  a  40  amp-2. 5  atm  xenon  arc  discharge. 
The  absolute  values  of  the  continuum  radiance  are  in  error  ( see  subsection  2. 2. 1) 
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Fig.  19-Spectral  radiance  at  the  surface  of  a  40  amp-2.5  atm 
xenon  arc  discharge 
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Fig.  20— Spectral  radiance  at  the  surface  of  a  40  amp-2. 5  atm 
krypton  arc  discharge 
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Fig.  21— Spectral  radiance  of  the  7606 A  line  at  the  surface  of  a 
40  amp-2. 5  atm  krypton  arc  discharge 
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Fig.  25-Spectral  overlap  integrand  for  the  W-Nd:  YAG  system  for  various  slab  thickness 
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Fig.  26-Spectral  overlap  integrand  for  the  K-Hg-Nd:  YAG  system  for  a  1.2cm  slab  thickness 
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Fig.  29-Slab  efficiency  for  0.7%  Nd.YAG  being  pumped  by  four  lamps 


X  Value  for  the  Air-Cooled  K*Hg  Lamp  Estimated  from  Incomplete  Data 
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Fig.  30-Pump  efficiencies  for  four  lamps  pumping  a  0. 7%  Nd:YAG  slab 


Curve  592332-A 


Fig.  31-Measured  power  conversion  of  a  5x  50  mm  Nd:YAG  rod 
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Fig.  32- Activator  efficiency  for  four  lamps  pumping  a  0. 7%  NdrYAG  slab.  The  input  power  densities  have  been 

normalized  to  200  watts/cm2 


FIG.  33  Linear  plot  of  the  experimental  and  theoretical  values  for  the  spectral  radiance  of  a  117  cm 
thick  xenon  plasma  of  homogeneous  temperature  corresponding  to  the  6400.1  -  44X0  A  cm*  peak 
current  density  series.  Flashtubc  :  1-27  cm  inside  dia..  30  cm  long,  filled  to  150  torr  xenon. 


Spectral  R 


34  Time  history  of  8231  6  A  line  of  xenon,  780J  input.  127  mm  dia,  tube.  150  torr  initial  pressure. 
Peak  current  1140  amps.  Peak  electric  field.  29  V/cm. 
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Fig.  35  Spectrophoto metric  absorption  spectrum  of  ruby.  The 
absorption  cross  section  <r  (right  scale)  is  given  by  a*=a/No,  where 
a  is  the  linear  absorption  coefficient  and  No  the  number  of  ab¬ 
sorbing  centers  per  cm8. 


Fig.  36  Fluorescent  quantum  and  power  efficiencies  for  ruby: 
ratio  of  output  fluorescence  to  input  absorption  as  a  function  of 
wavelength. 
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Fig.  37-Slab  efficiency  for  a  xenon  flash  lamp  pumping  a  ruby  slab 


Fig.  38~Pump  efficiency 
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Fig.  39-Effective  lifetimes  of  excited  chromium  ions  in  a  ruby  laser  rod  pumped  by  a  xenon  flash  lamp 


Curve  592731 -A 


Fig.  40”Energy  level  diagram  of  ruby  with  lifetimes  at  300°K 
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Fig.  42  Newton' s  method  applied  to  two  parameters  which  are  almost  linearly  dependent 
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